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FOREWORD 

The  work  described  herein  was  performed  at  the  Aerojet  Liquid  Rocket 
Company  under  NASA  Contract  NAS  8-33553,  with  Mr.  Fred  Braam,  NASA-Marshall 
Space  Flight  Center,  as  Project  Manager.  The  ALRC  Program  Manager  was 
Mr.  Jeffery  W.  Salmon,  and  the  Project  Engineer  was  Mr.  Gregory  M.  Meagher. 

The  technical  period  of  performance  for  this  study  was  from  9  July  1979 
to  31  October  1980.  The  overall  objective  of  this  study  was  to  upgrade  per¬ 
formance  and  heat  transfer  analysis  techniques  for  the  dual-throat  and  dual¬ 
expander  combustion  chamber  concepts. 

The  author  wishes  to  acknowledge  the  efforts  of  the  following  personnel 
for  their  significant  contributions: 


Dick  Ewen,  ALRC 
Jerry  Pieper,  ALRC 
Jim  Ouey,  ALRC 


Gary  Nickerson,  SEA  Inc. 

Irwin  Alber,  Consultant 
Joe  Hoffman,  Purdue  University 


This  report  presents  the  Summary,  Section  I  -  Introduction,  Section 
II  -  Conclusions  and  Recommendations,  and  Sections  III,  IV,  V  and  VI  - 
Technical  Discussion. 
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SUMMARY 


A.  STUDY  OBJECTIVES  AND  SCOPE 

The  objective  of  this  study  was  to  expand  and  extend  the  analysis 
models  and  parametric  studies  previously  performed  for  both  the  dual-throat 
and  dual-expander  engine  concepts.  This  was  accomplished  within  five  basic 
tasks,  including  (1)  improvements  to  the  existing  dual-throat  aerodynamic  and 
performance  prediction  computer  model,  (2)  preliminary  geometric  analysis  of 
the  dual-expander  concept,  (3)  preliminary  flow  field  analyses  of  the  dual¬ 
expander  concept,  (4)  further  preliminary  heat  transfer  analysis  of  both  con¬ 
cepts,  and  (5)  engineering  analysis  of  data  from  the  NASA/MSFC  hot-fire 
testing  of  a  dual-throat  thruster  model  thrust  chamber  assembly.  As  a  result 
of  this  effort,  a  quantitative  comparison  of  the  aerodynamic,  performance,  and 
thermal  design  characteristics  between  the  two  dual -nozzle  concepts  can  be 
derived. 


B.  STUDY  RESULTS 

During  this  study,  analytical  models  to  predict  performance  and 
operating  characteristics  of  dual-nozzle  concepts  have  been  developed  and 
improved.  Aerodynamic  models  are  now  available  to  define  flow  characteristics 
and  bleed  requirements  for  both  dual-throat  and  dual-expander  concepts. 
Advanced  analytical  techniques  were  utlized  to  provide  quantitative  estimates 
of  the  bleed  flow,  boundary  layer,  and  shock  effects  within  dual-nozzle 
engines.  Thermal  analyses  were  performed  to  define  cooling  requirements  for 
baseline  conf igurations,  and  special  studies  of  unique  dual-nozzle  cooling 
problems  (i.e. ,  lip  and  secondary  throat  regions)  defined  feasible  means  of 
achieving  adequate  cooling.  Similar  analysis  techniques  for  the  dual¬ 
expander  concept  are  now  available.  However,  these  techniques  are, 
in  general,  less  developed  than  those  for  the  dual -throat  concept 
and  have  not  been  substantiated  with  empirical  aerodynamic  correla- 
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Summary  (cont.) 


tions.  A  summary  of  these  results  is  contained  in  the  following  paragraphs, 
and  the  technical  details  are  provided  in  the  technical  discussion  of  this 
report. 


The  dual -throat  aerodynamic  model  has  been  expanded  to  include 
effects  not  considered  prior  to  the  subject  contract.  Previously,  the  aero¬ 
dynamic  model  was  limited  in  its  ability  to  predict  the  base  pressure  that 
developed  within  the  exterior  thrust  chamber.  As  presented  in  Section  III,B, 
the  aerodynamic  model  computer  program  has  been  modified  to  include  a  proce¬ 
dure  for  predicting  base  pressure.  This  model  assumes  that  the  base  pressure 
is  controlled  by  the  loss  in  total  pressure  which  occurs  when  the  shear  layer 
flow  passes  through  an  oblique  shock  during  attachment  at  the  exterior 
thruster  wall.  An  empirical  constant,  called  the  Nash  factor,  is  used  to 
calibrate  the  model  with  experimental  data.  A  further  development  of  the 
computer  program  has  been  to  provide  for  dissimilar  exhaust  gases  to  be  used 
for  the  interior  and  exterior  (bleed)  flow.  Each  stream  is  allowed  to  possess 
a  given  constant  value  for  specific  heat  at  constant  pressure,  molecular 
weight,  and  total  temperature. 

As  a  result,  the  aerodynamic  model  has  been  generalized  so  that  it 
can  be  applied  to  studies  of  various  dual-throat  concepts,  including  tripro¬ 
pellant  earth-to-orbit  engines  and  bipropellant  engine  applications  such  as 
the  Orbit  Transfer  Vehicle  Engine  (OTV-E)  and  liquid  oxygen/hydrocarbon 
boosters  for  uprating  the  Space  Shuttle.  This  model  was  also  used  to  define 
test  conditions  and  predictions  in  support  of  the  MSFC  hot-fire  testing  of  a 
dual -throat  thruster  model.  Subsequent  analysis  of  these  test  results,  as 
discussed  in  Section  VI,  has  shown  excellent  agreement  between  model  predic¬ 
tions  and  test  measurements. 


Summary  (cont.) 


The  capability  to  optimize  the  dual -throat  engine  performance  was 
enhanced  through  the  development  of  a  computer  program  that  allows  the  design 
of  an  optimum  nozzle  wall  contour  for  a  given  approach  flow.  This  was 
achieved  by  modifying  an  existing  Rao  nozzle  computer  program,  as  described  in 
Section  III,E.  As  a  result,  it  is  now  possible  to  design  the  contour  down¬ 
stream  of  the  exterior  nozzle  throat  so  that  it  will  provide  maximum  perform¬ 
ance  during  Mode  II  operation.  In  order  to  accurately  predict  the  performance 
of  nozzles  designed  using  these  procedures,  revisions  were  also  made  to  the 
TDK  computer  program  and  to  the  aerodynamic  model  computer  program.  This  pro¬ 
cedure  was  used  to  define  two  secondary  nozzle  contours  for  the  baseline  con¬ 
figuration:  one  optimized  for  booster  phase  operation  (Mode  I)  and  one  opti¬ 
mized  for  sustainer  operation  (Mode  II).  In  this  case,  the  performance  for 
the  operating  mode  in  which  the  nozzle  contour  was  optimized  was  about  one  (1) 
second  greater  than  its  performance  with  its  non-optimum  nozzle.  Although  a 
minor  difference  in  this  case,  it  shows  that  the  capability  now  exists  to 
optimize  the  performance  for  either  mode  of  operation  if  this  proves  to  be 
advantageous  from  a  total  mission  standpoint. 

Another  analysis  task  conducted  during  this  study  assessed  the 
effects  of  shock  waves  on  nozzle  performance.  Shock  waves  are  introduced  in 
the  supersonic  expansion  by  the  shape  of  the  jet  boundary  and  by  its  attach¬ 
ment  to  the  exterior  thruster  wall.  The  resulting  shock  structure  is  inher¬ 
ently  two-dimensional  and,  thus,  a  two-dimensional  numerical  procedure  must  be 
used  to  predict  the  flow  field.  The  VNAP  computer  program  was  used  for  this 
purpose.  The  results  of  this  study  task,  presented  in  Section  III,D,  indicate 
that  shocks  will  have  only  a  minor  influence  on  the  performance  and  thermal 
operating  characteristics  with  properly  designed  bleed  flows.  This  analysis 
was  also  used  to  define  boundary  conditions  for  the  thermal  analysis  of  the 
adverse  condition  of  Mode  II  operation  without  any  bleed  flow. 

The  boundary  layer  analysis  and  performance  loss  of  dual-throat 
engines  was  investigated  further  by  using  the  JANNAF  Reference  computer  code 
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BLIMP,  in  conjunction  with  the  aerodynamic  shear  model  (described  in  Section 
III.C).  This  procedure  provided  a  rigorous  method  of  assessing  boundary  layer 
influences,  including  the  effects  of  bleed  flow  addition,  and  guided  an 
improvement  to  the  simplified  boundary  layer  loss  procedure. 

Finally,  all  of  the  aforementioned  improvements  were  incorporated 
into  the  existing  dual-throat  simplified  performance  model,  and  results  were 
compared  to  previous  performance  predictions,  as  noted  in  Section  1 1 1 ,  F .  In 

general,  the  revised  model  predicts  lower  specific  impulse  (*wl% )  for  the  base¬ 
line  configuration  in  Mode  II,  but  as  a  result  of  this  investigation,  the 
uncertainty  of  this  prediction  has  been  significantly  reduced.  Parametric 
trends  with  operating  conditions  proved  similar;  however,  the  variation  of 
specific  impulse  with  thrust  split  is  altered  and  should  be  considered  in 
future  mission  analysis  studies. 

As  a  second  phase  of  this  study,  discussed  in  Section  IV,  perform¬ 
ance  models  for  dual -expander  thrusters  were  developed.  The  dual  expander  is 
similar  in  concept  to  the  dual-throat  thruster  and  consists  of  an  exterior, 
annular  thrust  chamber  that  surrounds  a  conventional  thrust  chamber.  Two  modes 
of  operation  are  treated:  Mode  I,  the  booster  mode  of  operation,  during  which 
both  engines  are  on,  and  Mode  II,  the  sustainer  inode  of  operation,  during 
which  the  annular  thruster  operates  alone.  During  Mode  II  operation,  a  free 
shear  layer  develops  along  the  jet  boundary  and  is  subsequently  turned  paral¬ 
lel  to  the  axis  of  symmetry  and  transformed  into  a  turbulent  wake.  Bleed  flow 
can  be  used  to  control  the  shape  and  properties  of  this  viscous  region. 

An  aerodynamic  model  similar  in  concept  to  the  dual-throat  model 
was  developed  to  determine  the  shear  layer  and  base  recirculation  flow  fields 
during  Mode  II  operation.  Use  of  this  model  showed  that  the  bleed  flow 
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Summary  (cont.) 


requirements  for  a  dual-expander  nozzle  were  significantly  influenced  by  the 
truncation  point  (Mach  number)  of  the  inner  nozzle.  The  overall  bleed  flow 
required  was,  in  general,  0  to  3%,  which  is  significantly  less  than  that 
required  for  the  dual  throat.  In  fact,  it  was  concluded  that  the  dual¬ 
expander  engine  could  be  designed  and  operated  without  bleed  flow  while  still 
obtaining  controllable  aerodynamic  flow  characteristics. 

A  model  for  predicting  shock  losses  and  flow  effects  during  Mode 
II  operation  was  also  applied  to  predict  shock  location  and  performance 
effects  with  varying  inner  nozzle  truncation  and  bleed  flow.  The  results 
indicated  little  performance  differences  with  nozzle  truncation  and  no  signi¬ 
ficant  shock  influences  over  the  range  of  bleed  flows  considered.  However, 
the  VNAP  results  were  not  sufficiently  accurate  in  its  present  form  to  be  used 
for  performance  prediction  or  for  optimum  nozzle  design.  To  overcome  this 
problem,  development  of  a  method  of  characteristics  analysis  with  shock  capa¬ 
bility  is  recommended.  Most  of  the  point  procedures  required  for  this 
approach  are  available  and  could  be  used  in  assembling  the  new  computer  pro¬ 
gram. 


Procedures  were  also  developed  to  specify  preliminary  geometric 
design  parameters  for  a  dual-expander  nozzle.  These  procedures  are  based  on 
the  criterion  that  the  flow  angle  and  static  pressure  must  be  matched  at  the 
primary  nozzle  lip  during  Mode  I  operation.  During  this  study,  the  VNAP 
analysis  was  used  to  calculate  the  flow  characteristics  in  the  annular  nozzle 
region.  The  subsequent  results  verified  that  the  design  criteria  had  been 
met. 


Thermal  design  and  parametric  analyses,  described  in  Section  V, 
were  performed  for  both  dual-nozzle  concepts.  A  satisfactory  dual-throat 
baseline  design  was  obtained  by  using  three  parallel  cooling  circuits. 


S 


Summary  (cont.) 


The  primary  circuit  coolant  flows  from  an  area  ratio  of  1.58:1  in  the  primary 
nozzle  to  the  primary  injector.  This  circuit  has  the  highest  pressure  drop 
requirement  of  the  three  parallel  circuits,  so  its  coolant  flow  fraction  is 
selected  to  provide  the  minimum  pressure  drop.  A  separate  circuit  cools  the 
inner  annulus  surface,  the  lip  region,  and  the  rest  of  the  primary  nozzle. 

The  majority  of  the  channel  pressure  drop  in  this  circuit  is  required  for  the 
lip  region,  for  which  a  three-dimensional  SINDA  model  was  developed.  A 
counterflow  circuit  cools  he  outer  surface  of  the  secondary  chamber;  this 
circuit  is  a  composite  design  which  considers  both  Mode  I  and  Mode  II  oper¬ 
ation.  The  channel  design  over  a  short  region  upstream  of  the  secondary 
throat  is  defined  by  the  primary  plume  impingement  associated  with  a  bleed 
flow  malfunction  in  Mode  II.  An  exact  solution  for  the  inviscid  plume 
impingement  problem  was  obtained  with  the  VNAP  program.  This  solution 
provided  edge  conditions  for  an  integral  boundary  layer  analysis  which  defined 
the  wall  heat  transfer  edge  conditions.  The  initial  momentum  thickness  for 
the  boundary  layer  analysis  was  varied  in  order  to  bound  the  problem,  with  the 
channel  design  based  on  the  worst  case.  The  maximum  heat  flux  ranged  from 
essentially  the  Mode  I  value  to  almost  50%  higher. 

Individual  cooling  circuit  characteristics  for  the  dual-throat  baseline 
design  have  been  defined.  Although  the  fuel  system  pressure  drop  of  1160  psi 
is  set  by  the  primary  circuit,  the  coolant  Mach  number  is  slightly  higher  in 
the  secondary  circuit  when  the  latter  is  designed  to  accommodate  the  maximum 
postulated  Mode  II  heat  flux  for  the  zero  bleed  flow  case.  Parametric  study 
results  indicate  that  it  is  the  coolant  Mach  number,  rather  than  pressure 
drop,  which  limits  increases  in  chamber  pressure  and  Mode  1  thrust  split. 
Increasing  the  chamber  pressures  to  4500/3150  psia  or  the  thrust  split  to 
80/20  results  in  a  primary  circuit  Mach  number  of  0.37  (an  unacceptably  high 
design  value). 


Summary  (cont.) 


The  dual-expander  designs,  similar  to  those  of  the  dual  throat,  are 
cooled  with  three  parallel  circuits.  However,  the  inner  annulus  circuit  is 
now  a  major  circuit  in  terms  of  length  and  heat  load  and  also  includes  a 
throat  region.  Regenerative  cooling  of  the  baselne  design  point  with  hydrogen 
only  is  difficult.  Coolant  flows  in  the  secondary  and  inner  annulus  circuits 
must  be  selected  to  avoid  minimum  channel  depth  and  Mach  number  limits, 
respectively,  and  the  remaining  flow  available  for  the  primary  is  less  than 
optimum.  As  a  result,  the  system  pressure  drop  of  1720  psi  is  130  psi  above 
the  minimum  primary  circuit  value.  Using  oxygen  to  cool  the  secondary  circuit 
and  reducing  the  interface  area  ratio  between  the  primary  circuit  inlet  and 
the  inner  annulus  circuit  outlet  allows  the  hydrogen  pressure  drop  to  be 
reduced  to  1550  psi.  Transpirat ion  cooling  the  throat  region  of  the  primary 
chamber  did  not  prove  to  be  effective  since  the  hydrogen  flow  fractions 
required  for  transpiration  cooling,  coupled  with  the  flow  requirements  of  the 
other  circuits,  severely  limit  the  flow  available  for  regeneratively  cooling 
the  remainder  of  the  primary  circuit.  This  results  in  relatively  high  pressure 
drops. 


Parametric  study  designs  of  the  dual  expander  could  not  be  obtained 
with  increased  chamber  pressures  or  increased  Mode  I  thrust  split,  or  with  the 
secondary  chamber  pressure  reduced  to  1500  psia.  The  latter  result?  from  the 
high  primary  chamber  pressure  drop  associated  with  reduced  coolant  inlet  pres¬ 
sure.  These  results,  along  with  the  baseline  design  details,  indicate  that 
the  baseline  design  point  is  near  the  optimum  chamber  pressure  ratio  and  is 
very  close  to  the  chamber  pressure  versus  thrust  split  characteristic  which 
represents  a  practical  cooling  limit. 


I. 


INTRODUCTION 


Propulsion  systems  for  future  vehicles,  such  as  the  Si ngle-Stage-to- 
Orbit  (SSTO)  and  Heavy  Lift  Launch  Vehicle  (HLLV),  may  embrace  such  capabili¬ 
ties  as  dual-mode  operation  and  in-flight  changes  in  area  ratio  for  altitude 
compensation.  These  vehicles  benefit  from  dual-mode  operation  through  reduced 
vehicle  volume  by  taking  advantage  of  relatively  high  bulk  density  propellants 
in  one  mode  and  lower  density,  but  higher  performing  propellants  in  the  other 
mode.  Area  ratio  change  during  flight  provides  an  increase  in  performance  as 
ambient  pressure  decreases  with  altitude. 

Dual-nozzle  engines  combine  both  operating  capabilities  in  a  single 
design.  Their  dual  combustors  allow  use  of  either  one  or  two  propellant  com¬ 
binations  and,  with  their  two  separate  nozzle  throats  and  a  fixed  nozzle  exit 
area,  allow  for  a  shift  in  area  ratio  without  resorting  to  nozzle  translating 
mechanisms. 

The  dual-nozzle  concept  can  be  applied  to  both  bipropellant  and  tripro¬ 
pellant  engines  with  resultant  advantages.  The  baseline  engine  application, 
analzyed  in  this  and  earlier  studies,  is  a  tri propel lant  si ngl e-stage-to-orbit 
(SSTO)  engine.  Dual  nozzles  could  also  be  applied  to  bipropellant  dual-thrust 
Orbit  Transfer  Vehicle  (OTV)  engine  and  Liquid  Rocket  Booster  (LRB)  engine 
designs.  For  OTV-E,  a  high-thrust  priority  payload  and  a  low-thrust  cargo 
transfer  engine  could  be  combined  in  a  single  thrust  chamber.  For  a  LOX/ 
Hydrocarbon  LRB,  mission  thrust  tailoring  and  an  area  ratio  increase  at  alti¬ 
tude  could  be  achieved.  This  concept  has  been  analyzed,  with  favorable 
results,  during  the  recent  Advanced  Oxygen-Hydrocarbon  Rocket  Engine  Study 
(Ref.  1). 

Two  types  of  dual-nozzle  designs  have  been  conceived  --  the  dual  throat 
and  the  dual  expander.  An  engine  system  preliminary  analysis  using  the  dual¬ 
throat  concept  was  performed  during  Contract  NAS  8-32967  to  examine  potential 


I,  Introduction  (cont.) 


power  cycles  and  generate  parametric  data  for  a  tripropellant  SSTO  vehicle 
engine.  A  preliminary  performance  prediction  methodology  based  on  a  new  aero¬ 
dynamic  bleed  flow  computer  model  was  developed.  This  model  was  formulated  by 
using  the  results  of  cold-flow  tests  conducted  with  a  subscale  dual-throat 
thruster  configuration.  The  dual-expander  engine  concept,  prior  to  the  sub¬ 
ject  contract,  has  received  less  formal  analysis.  Some  preliminary  engine 
system  parametric  data  have  been  generated,  and  a  conceptual  baseline  engine 
system  has  been  established. 

A  dual-throat  thruster  consists  of  a  large  thrust  chamber  assembly  that 
contains  within  it  a  smaller  thrust  chamber.  Because  of  this  novel  construc¬ 
tion,  dual-throat  thrusters  provide  the  means  to  obtain  a  large  area  ratio 
adjustment  within  a  single  thrust  chamber  assembly  without  the  need  for 
extendible  nozzles.  This  area  ratio  adjustment  is  accomplished  by  changing 
the  operating  mode  of  the  thruster.  For  example,  in  the  high-thrust  mode,  the 
exterior  and  interior  thrust  chambers  are  both  operating.  When  the  exterior 
thrust  chamber  is  extinguished,  the  low-thrust  interior  thrust  chamber  remains 
in  operation,  but  now  utilizes  the  exterior  thrust  chamber  exhaust  nozzle  as 
its  exhaust  expansion  contour.  In  this  manner,  the  interior  thruster  is 
caused  to  operate  at  a  high  expansion  ratio.  The  attachment  of  the  interior 
thruster  exhaust  to  the  exterior  thruster  wall  contour  can  be  controlled  by 
injecting  mass  flow  into  the  exterior  chamber,  thereby  creating  a  free  shear 
layer  that  develops  along  the  jet  boundary.  Mass  injection  will  raise  the 
base  pressure  in  the  chamber  of  the  exterior  thruster  and  move  the  jet 
attachment  point  downstream.  Proper  location  of  the  free  jet  can  serve  to 
lower  the  losses  in  interior  chamber  performance  that  are  caused  by 
aerodynamic  effects. 

The  sustainer  mode  of  operation,  as  described  above,  is  called  Mode  II, 
ai.u  the  booster  mode  of  operation  is  called  Mode  I.  The  dual-throat  thrust 
chamber  configuration  and  flow  fields  for  Mode  I  and  II  of  the  tripropellant 
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I,  Introduction  (cont.) 
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SSTO  engine  case  are  illustrated  in  Figure  1-1.  During  Contract  NAS  8-32967, 
the  baseline  engine  design  gas  generator  flow  was  employed  as  the  Mode  II 
bleed  flow.  During  Mode  I,  the  gas  generator  flow  would  be  dumped  in  the 
nozzle  exhaust. 

The  dual  expander  is  similar  in  concept  to  the  dual-throat  thruster. 

The  dual  expander  consists  of  an  exterior  annul uar  thrust  chamber  that  sur¬ 
rounds  a  conventional  thrust  chamber.  Figure  1-2  illustrates  the  two  modes  of 
operation  for  the  SSTO  engine  application  that  are  treated  in  this  study. 

Mode  I  is  the  booster  mode  of  operation  during  which  both  engines  are  on.  The 
annular  thruster  uses  a  LOX/H2  propellant  combination,  and  the  interior 
thruster  uses  LOX/RP-I.  Thus,  in  the  dual  expander,  the  exterior  chamber  is 
the  sustainer  thruster  whereas,  in  the  dual-throat  design,  the  interior 
chamber  is  the  sustainer  thruster.  Mode  II  is  the  sustainer  mode  of  operation 
during  which  the  annular  thruster  operates  alone.  During  Mode  II  operation,  a 
free  shear  layer  develops  along  the  jet  boundary  and  is  subsequently  turned 
parallel  to  the  axis  of  symmetry  and  transformed  into  a  turbulent  wake.  Bleed 
flow  can  be  used  to  control  the  shape  and  properties  of  this  viscous  region. 

Technological  investigation  of  the  dual -throat  nozzle  concept  was  ini¬ 
tiated  in  1975  by  the  NASA  George  C.  Marshall  Space  Flight  Center.  In  July  of 
1976,  the  NASA  and  the  Aerojet  Liquid  Rocket  Company  (ALRC)  collaborated  on  an 
analytical/empirical  program  to  evaluate  the  concept.  A  bailment  agreement 
was  established  wherein  NASA  hardware  was  supplied  for  use  on  an  ALRC- 
sponsored  cold-flow  program.  Aerodynamic  analyses  were  conducted,  and  the 
cold-flow  test  program  was  completed.  After  the  promising  aerodynamic  fea¬ 
tures  had  been  confirmed,  NASA  contracted  ALRC  (NAS  8-32666,  Ref.  1)  to  expand 
the  data  base  and  to  initiate  development  of  a  dual-throat  performance  predic¬ 
tion  model . 
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I,  Introduction  (cont.) 

The  Dual-Fuel,  Dual-Throat  Engine  Preliminary  Analysis  Program  (NAS 
8-32967,  Ref.  2)  was  subsequently  conducted  to  determine  dual-throat  engine 
system  parameters  and  vehicle/engine  integration  requirements  for  the  tri- 
propellant  SSTO  engine  application.  This  study  concluded  that  while  the 
dual-throat,  dual-fuel  engine  was  a  viable  SSTO  candidate,  a  number  of  uncer¬ 
tainties  remaining  in  the  heat  transfer  and  performance  areas  needed  to  be 
resolved.  Also,  as  the  Dual-Expander  Engine  had  been  conceived  during  this 
period  (Ref.  3),  there  was  a  need  to  make  accurate  trades  between  the  two 
dual-nozzle  concepts.  As  a  result,  the  effort  described  in  this  report  was 
conducted. 

To  accomplish  this,  the  Dual  Nozzle  Aerodynamic  and  Cooling  Analysis 
study  was  conducted.  The  results  are  discussed  in  this  report.  The  major 
objectives  of  the  program  were  as  follows: 

(1)  Improve  the  existing  dual-throat  aerodynamic  bleed  flow  model  and 
performance  prediction  methodology; 

(2)  Develop  basic  geometric  and  aerodynamic  data  defining  performance 
for  a  dual -expander  chamber; 

(3)  Conduct  preliminary  thermal  analyses  of  both  the  dual  throat  and 
dual  expander; 

(4)  Provide  engineering  analysis  of  data  from  government  hot-fire 
testing  of  a  dual-throat  thruster  model  chamber. 
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I,  Introduction  (cont.) 

The  objectives  were  accomplished  through  a  series  of  five  tasks: 

Task  I:  Dual-throat  aerodynamic  model  and  performance  prediction 

improvements.  This  task  resulted  in  improving  the  bleed 
flow  model  to  analytically  predict  base  pressure  for 
dissimilar  compositions  of  the  primary  and  bleed  flows. 
Secondly,  a  computer  program  was  developed  to  optimally 
design  the  secondary  nozzle  contour  for  Mode  II  operation. 
Third,  simplified  procedures  used  to  predict  boundary  layer 
and  divergence  losses  were  verified  by  a  rigorous  analysis. 
This  analysis  also  investigated  the  influence  of  shocks  on 
engine  operation. 

Task  II:  Dual -expander  preliminary  geometric  analysis.  A  design  cri¬ 
terion  was  established,  and  the  interrelationships  between 
nozzle  geometry  and  engine  operating  characteristics  were 
defined. 

Task  III:  Dual-expander  preliminary  flow  field  analysis.  The  flow 
field  characteristics  of  the  dual-expander  operation  were 
investigated,  verifying  the  design  approach  from  Task  II  and 
identifying  the  insensitivity  of  nozzle  performance  to  inner 
chamber  truncation  and  the  low  bleed  flow  requirements  of 
the  dual  expander. 

Task  IV:  Preliminary  thermal  analysis  of  both  dual  throat  and  dual 
expander.  This  task  showed  that  the  lip  regions  of  the 
inner  nozzle  in  both  concepts  could  be  cooled.  Completely 
regenerati vely  cooled  designs  were  established  for  both 
engines. 

Task  V:  Test  data  analysis.  Hot-fire  data  analysis  verified  the 

improved  aerodynamic  bleed  flow  model. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  major  conclusions  formed  during  conduct  of  the  program  are 
listed  below. 

1.  The  dual-throat  aerodynamic  bleed  flow  model  has  been  veri¬ 
fied  on  the  basis  of  cold-flow  and  hot-fire  data  correla¬ 
tions  and  is  ready  for  use  as  a  general  design  tool. 

2.  The  predictive  accuracy  of  the  aerodynamic  bleed  flow  model 
could  be  bettered  by  improving  the  current  plume  shape  calcu 
lational  procedure. 

3.  Dual-throat  engine  performance  is  relatively  insensitive  as 
to  whether  or  not  the  secondary  nozzle  contour  is  optimized 
for  the  Mode  I  or  Mode  II  expansion. 

4.  The  previous  simplified  technique  for  predicting  the  dual¬ 
throat  Mode  II  boundary  loss  was  inaccurate.  The  rigorous 
procedure  employed  during  this  program  has  resulted  in 
significant  improvement  of  the  simplified  prediction  proce¬ 
dure. 

5.  Dual-throat  nozzle  performance  in  Mode  I  is  comparable  to 
that  of  a  conventional  nozzle.  Mode  II  performance  is 
approximately  1-2%  lower  than  that  of  a  conventional  nozzle 
optimized  for  altitude  ( i . e. ,  sustainer)  operation. 


II,  A,  Conclusions  (cont.) 


6.  Preliminary  performance  analysis  results  for  the  dual¬ 
expander  nozzle  indicate  that  high  Mode  II  expansion  effi¬ 
ciencies  can  be  achieved.  The  results  indicate  low  bleed 
flow  requirements  and  insensitivity  to  the  primary  nozzle 
truncation  area  ratio. 

7.  There  is  a  degree  of  uncertainty  associated  with  the  dual¬ 
expander  performance  results.  The  cal culational  procedures 
employed  during  this  program  (i.e.,  the  VNAP  computer  model) 
should  be  replaced  by  a  method  of  characteristics  solution 
that  treats  shock  waves. 

8.  Both  the  dual-throat  and  dual-expander  baseline  tripropellant 
SSTO  engine  designs  can  be  cooled  with  three-circuit  regener¬ 
ative  designs.  The  primary  nozzle  lip  and  secondary  nozzle 
throat  gas-side  heat  transfer  coefficient  assumptions  should 
be  verified  through  hot-fire  testing. 

B.  RECOMMENDATIONS 

The  following  recommendations  for  future  work  are  made  on  the 
basis  of  the  significant  results  obtained  on  this  program.  They  are  listed  in 
the  order  of  highest  to  lowest  priority. 

1.  A  demonstration  program  should  be  conducted  to  verify  the 
major  thermal  design  assumptions  utilized  during  the  study. 
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II,  B,  Recommendations  (cont.) 


2.  The  dual-expander  flow  field  analysis  results  should  be  veri¬ 
fied  through  conductance  of  a  model  cold-flow  test  program. 

3.  The  dual-throat  SSTO  engine  parametrics  developed  on  Contract 
NAS  8-32967  should  be  updated  on  the  basis  of  the  improved 
performance  prediction  results  obtained  on  this  program. 

4.  A  method  of  characteristics  solution  that  treats  flow  field 
shocks  should  be  developed  for  Mode  II  dual  expander 
analysis. 

5.  The  plume  shape  calculational  procedure  contained  within  the 
aerodynamic  bleed  flow  model  should  be  improved  and  corre¬ 
lated  with  existing  test  data. 

6.  A  procedure  for  optimization  of  dual-expander  secondary 
nozzle  contour  should  be  developed. 
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III.  DUAL  THROAT  AERODYNAMIC  MODEL  AND  PERFORMANCE  PREDICTION 

Improvement 


The  basic  objectives  of  Task  I  were  to  improve  the  dual  throat  aero¬ 
dynamic  model's  capability  and  verify  the  approaches  used  in  the  dual  throat 
performance  prediction  methodology  developed  by  ALRC  for  NASA/MSFC  during 
Contract  NAS  8-32666. 

The  specific  objectives  of  Task  I  were  as  follows: 

°  Improve  the  base  flow  pressure  and  shear  layer  portion  of  the 
aerodynamic  bleed  flow  model. 

°  Assess  the  approach  used  to  calculate  boundary  layer  loss  for  Mode 
II  operation. 

0  Assess  the  approach  used  to  calculate  nozzle  divergence  efficiency 
for  Mode  II  operation  and  determine  how  shocks  affect  it. 

°  Develop  a  Rao  nozzle  optimization  program  for  the  Mode  II  opera¬ 
tion  and  determine  the  effect  of  Mode  II  contour  optimization  on 
the  delivered  performance  of  both  modes. 

0  Assess  changes  in  predicted  performance  based  on  the  use  of  the 
improved  aerodynamic  model  and  update  the  performance  model  as 
required. 

°  Review  the  existing  cold-flow  data  bank  and  define  areas  that 
would  benefit  from  further  cold-flow  test  data. 
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III,  Dual  Throat  Aerodynamic  Model  and  Performance  Prediction 
Improvement  (cont.) 

A.  LITERATURE  SURVEY 

In  support  of  the  improvement  of  the  aerodynamic  model,  a  compu¬ 
terized  literature  survey  was  conducted  for  the  purpose  of  obtaining  refer¬ 
ences  that  might  be  useful  in  formulating  the  analytical  model  as  well  as  data 
that  could  be  used  for  calibration  of  the  model.  Conducting  a  computerized 
literature  search  has  the  advantage  of  significantly  reducing  the  amount  of 
man-hours  otherwise  required  for  this  task,  thus  significantly  reducing  cost 
and  time  as  well. 

The  literature  search  was  carried  out  via  a  remote  data  terminal 
to  gain  access  to  a  computerized  library.  The  method  used  involves  the  formu¬ 
lation  of  a  search  strategy  based  on  key  words  or  phrases.  These  may  be  com¬ 
bined  using  boolian  operators  to  obtain  all  literature  from  a  given  data  base. 
The  search  is  conducted  interactively  on  the  main  frame  computer.  For 
example,  reference  titles  that  contain  a  given  key  word,  but  not  another  key 
word,  may  be  listed  on  the  terminal.  If  a  particular  title  appears  to  be  of 
interest,  a  request  by  terminal  is  made  and  the  abstract  of  the  article  is 
printed.  If  desired,  the  article  can  then  be  ordered  over  the  terminal. 

In  order  to  prepare  a  search  strategy,  the  initial  computer  ses¬ 
sion  was  spent  in  checking  various  key  author  names  and  key  words  for  fre¬ 
quency  of  citation  in  the  NTIS  index.  Titles  obtained  in  this  way  were  used 
to  construct  the  boolian  operators  to  be  used  in  the  general  search.  A 
strategy  was  established  to  be  accomplished  in  two  steps. 

The  first  step  was  to  identify  those  papers  which  make  reference 
to  the  principal  works  in  the  field.  This  was  accomplished  with  the  Science 
Citation  Index  capability  of  SCISEARCH.  Papers  which  referred  to  the  works  of 
Korst,  Bauer,  and  Alber  were  sought  in  this  way. 
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Ill,  A,  Literature  Survey  (cont.) 


The  second  part  of  the  literature  search  utilized  various  key 
words  and  employed  six  data  bases.  Information  obtained  in  this  manner  is 
summarized  in  Table  III-I,  which  gives  the  number  of  references  found  by 
category  by  source.  Abstracts  were  obtained  for  these  references.  After 
reading  the  abstracts,  approximately  30  papers  were  ordered. 

B.  AERODYNAMIC  BLEED  FLOW  MODEL 

The  aerodynamic  bleed  flow  model's  capability  of  predicting  the 
aerodynamics  unique  to  Mode  II  operation  has  been  expanded.  The  Mode  II  (sus- 
tainer  mode)  operation  is  illustrated  in  Figure  1 1 1 -1 .  The  interior  thrust 
chamber  is  called  the  primary  thruster,  and  its  exhaust  is  called  the  primary 
flow  stream.  A  secondary  flow  stream  (bleed  flow)  is  injected  into  the 
exterior  thrust  chamber.  The  thruster  expands  the  primary  flow  to  a  super¬ 
sonic  condition  at  the  nozzle  exit.  Further  expansion  then  occurs  in  the  form 
of  a  Prandtl-Meyer  fan  at  the  nozzle  lip.  The  result  is  an  exhaust  plume  with 
a  constant  pressure  boundary.  Flow  is  introduced  into  the  exterior  chamber  so 
as  to  control  the  location  of  the  plume  boundary  in  order  to  minimize 
reattachment  shocks  which  occur  when  the  flow  impinges  on  the  exterior 
thruster  wall.  The  plume  boundary  is  an  exhaust  streamline  path  and  acts  very 
much  like  a  nozzle  wall.  However,  a  shear  layer  develops  along  this  boundary 
due  to  viscous  interaction  of  the  exhaust  jet  and  the  gases  recirculating  in 
the  exterior  chamber.  Analysis  of  this  shear  layer  mixing  region  and  of  the 
phenomena  described  above  are  an  important  feature  of  the  aerodynamic  model. 

The  principal  assumptions  employed  in  the  aerodynamic  model  are 
listed  below: 


The  flow  leaving  the  primary  nozzle  is  one-dimensional  and 
supersonic. 


SUMMARY  OF  LITERATURE  SEARCH-AEROJET  DUAL  NOZZLE  ANALYSIS 


Ill,  8,  Aerodynamic  Bleed  Flow  Model  (tout.) 

°  The  primary  nozzle  exhaust  flow  is  represented  by  an  ideal 
gas  with  constant  specific  heat  and  molecular  weight. 

°  The  bleed  flow  is  represented  by  an  ideal  gas  with  constant 
specific  heat  and  molecular  weight.  These  values  can  differ 
from  those  of  the  primary  exhaust  stream.  The  total  temper¬ 
ature  of  the  bleed  flow  can  also  differ  from  that  of  the  pri¬ 
mary  flow  stream. 

0  The  shear  layer  is  represented  by  a  model  of  the  type 
developed  by  Chapman,  Korst,  Bauer,  Chow  and  Addey,  and 
others  (see  Refs,  b,  6,  7,  8,  and  9).  The  model  used  here, 
however,  is  nor- i soenerget ic ,  as  defined  by  the  assumptions 
given  above. 

°  The  shear  layer  is  treated  as  two-dimensional.  This  assump¬ 
tion  has  been  shown  to  be  accurate  (see  Bauer,  Ref.  7),  pro¬ 
vided  that  the  projected  thickness  of  the  mixing  zone  on  the 
radius  of  symmetry  is  less  than  .3  of  the  radius  of  symmetry. 
This  condition  |>.  always  satisfied  in  dual  throat  thrusters. 

The  boundary  1  aver  at  the  primary  nozzle  lip  is  thin. 

The  shear  layer  >s  a  constant  pressure  surface,  both  in  cross 
section  and  along  its  on’ 're  length.  At  the  inner  edge  of 
the  shear  layer,  the  g,r.,  velocity  is  negligible.  At  the 
outer  edge  of  the  shear  layer,  the  velocity  is  constant  and 
equal  to  the  velocity  of  the  plume  boundary  streamline. 

The  shear  layer  is  superimposed  on  the  inviscid  jet  plume 
bounda  ry. 
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Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 


°  The  plume  boundary  can  be  represented  by  the  methods  devel¬ 
oped  by  Herron  (Ref.  10). 

1.  Shear  Layer  Modeling 

Some  features  of  the  free  shear  layer  that  develops  along  the 
jet  boundary  are  shown  in  Figure  1 1 1-2 -  The  shear  layer  is  assumed  to  begin 
at  the  exit  of  the  primary  nozzle  where  its  thickness  is  negligible.  It  then 
develops  as  a  turbulent  shear  layer  along  the  boundary  of  the  primary  nozzle 
exhaust  plume.  Viscous  interaction  of  the  exhaust  plume  with  gases  recircu¬ 
lating  in  the  chamber  of  the  secondary  nozzle  form  the  shear  layer.  In 
steady-state  operation,  the  system  will  maintain  a  stable  equilibrium  base 
pressure.  This  base  pressure  determines  the  plume  boundary.  The  theory 
assumes  that  a  streamline  divides  the  shear  layer  such  that  gases  on  one  side 
of  the  streamline  are  recirculated  into  the  base  region.  This  is  the 
so-called  "dividing  streamline"  concept.  At  steady-state  operation ,  all  of 
the  gases  exhausting  from  the  primary  nozzle  will  also  exhaust  from  the  secon¬ 
dary  nozzle.  Gases  in  the  recirculating  region  (in  this  case  the  chamber  of 
the  secondary  nozzle)  have  been  trapped  because  they  are  degraded  in  total 
pressure  to  the  extent  that  they  cannot  penetrate  the  static  pressure  rise 
produced  by  turning  off  the  flow  at  the  point  where  the  shear  layer  impinges 
on  the  secondary  nozzle  wall.  It  is  assumed,  in  this  model,  that  the  dividing 
streamline  becomes  a  stagnation  point  at  the  secondary  wall.  This  stagnation 
occurs  abruptly  either  by  a  compression  that  is  essentially  isentropic,  or  by 
a  non-i sentropic  compression  through  a  shock  structure.  Details  of  this  pro¬ 
cess  govern  the  value  that  will  be  obtained  for  a  base  pressure,  i.e.,  the 
pressure  along  the  plume  boundary. 

The  introduction  of  mass  flow  into  the  recirculation  region 
will  tend  to  raise  the  base  pressure  and  shift  the  plume  boundary  downstream. 
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Figure  III-3.  The  "d"  and  "j"  Streamlines  within  the  Shear  Layer 


Ill,  8,  Aerodynamic  Bleed  Flow  Model  (cont.) 


Mass  flow  introduced  into  the  recirculation  region  at  a  given  rate  must  exit 
the  nozzle  at  the  same  rate.  If  this  flowrate  is  sufficiently  small,  all  of 
the  injected  flow  will  become  entrained  within  the  shear  layer.  For  this 
reason,  the  analysis  used  here  is  called  a  "small  bleed  flow"  model. 

Figure  1 1 1 -3  shows  the  location  of  two  important  streamlines 
within  the  shear  layer.  The  “d"  streamline  is  the  so-called  "dividing  stream¬ 
line"  that  separates  the  recirculation  gases  from  the  gases  that  exit  the 
nozzle.  The  recirculation  gases  are  assumed  to  be  a  mixture  or  primary  and 
bleed  flow  gases.  The  "j"  streamline,  also  shown  in  Figure  1 1 1 -3 ,  is  a 
limiting  streamline  for  the  bleed  flow  gas.  On  the  primary  side  of  the  "j" 
streamline,  there  is  "pure"  primary  flow,  whereas,  on  the  secondary  side, 
there  is  "mixed"  primary  and  secondary  flow.  All  of  the  bleed  flow  introduced 
during  steady-state  operation  must  exit  the  nozzle  and  is  contained  between 
the  "d"  and  "j"  streamlines.  The  "d"  and  "j"  streamlines  coincide  if  the 
bleed  flow  is  reduced  to  zero. 

The  velocity  profile  at  any  position  across  the  shear  layer 
will  appear  as  shown  in  Figure  1 1 1-2.  The  shear  layer  will  be  turbulent,  but 
for  either  laminar  or  turbulent  flow  the  velocity  profile  can  be  represented 
by 


0  =  (1  +  erf  n)  /  2 


where 


0  =  u/u^  ,  0<  0  <1 

and  erf  is  the  error  function,  i.e., 

n 

erf  n  =  2  n  S  exp  (-n^)  dn 

o 


•  \ 
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Ill,  B,  Aerodynamic  Bleed  flow  Model  (con!..) 
n  is  the  nondimensional  mixing  length 


In  the  above  equation,  y  is  t^o  ^  ;  t  nr  <  ■(>  inured  normal  to  the  centerline  of 
the  shear  layer,  and  >'  is  the  !erg‘h  the  shear  layer  measured  along  the 
centerline  from  its  origin.  Thr  variable  ’s  known  as  the  "shape  factor"  or 
"similarity  parameter."  The  difference  between  laminar  and  turbulent  mixing 
is  reflected  in  the  rel at i onshi r  used  f - ■'  It  is  necessary  to  use  an  empir¬ 
ical  expression  for  the  shape  factor.  Numerous  expressions  for  r  are  avail¬ 
able  in  the  literature  (see  Refs.  b,  7,  9,  11),  but  the  empirical  expression 
given  by  Korst  (Ref.  6)  has  been  used  here,  i.e., 

n  --  12  +  2.758  M. 

b 

where  Mb  is  the  Mach  number  along  The  inviscid  primary  jet  boundary. 

Since  the  bleed  flow  gas  properties  can  be  different  from 
those  in  the  primary  stream,  if  necessary  to  assume  a  variation  in  pro¬ 
perties  across  the  mixing  /one,  .<■*.,  across  the  shear  layer  from  the  base 
reci rculation  region  to  the  "g"  sfreaml ’ ne.  It  is  assumed  that  the  distribu¬ 
tion  of  mass  fraction  across  tin.  > 1  ng  /one  is  linear  with  respect  to  velo¬ 
city  ratio,  thus: 


n 


in 


!’-*)+  0C 

PP 


where  Cpr  is  the  specific  heat  m  the  recirculation  region,  i.e., 


r  : -v  )  +  C  Y 

pr  ps  pr '  pp  pr 

Y  is  the  amount  of  mass  from  fhe  primary  flow  entrained  in  the  recirculation 
pr 

region . 
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Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 
The  gas  constant  in  the  mixing  region  is 


R  =  Rr  (1-0)  +  0Rp 

where  Rr  is  the  gas  constant  in  the  recirculation  region,  i.e., 


R„  =  R. 


with 


(1-Y  )  +  R  Y 

pr;  p  pr 


Rs  =  R/Mws 

and 

Rp  =  R/MWp 

where  R  is  the  universal  gas  constant.  The  specific  heat  ratio  in  the  mixing 
zone  is 


Using  the  relations 


it  follows  that 


Y 


r=  y<cp-R), 


Y  *  C  /  C 

p  V 


Yr  (1-0)  +  Yp 


0 


(1-0)  + 


0 
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Ill,  8,  Aerodynamic  Bleed  Mow 


'•',;del  (cunt.) 


In  like  manner,  the  specific  hea* 
to  be 


v(. 


ratio  along  the  dividing  streamline  is  found 


M-rfJ  + 


'  r 
'f 


*A 


Next,  the  velocity  ratio  at  the  dividing  streamline,  ^  >  is  sought.  Tbr 
energy  equation  is 


C  T  *  '? !?  -  C  T 
p  no 


so  that 


y2  -  2CP  (T,  -  T) 


and 


M  <Tod  ‘  V/cpp  !top  '  V 

where  is  the  temperature  alone  fhe  plume  boundary.  Defining  the  Crocco 
number  as 

r:-  I  -  T/T 


gives 


»/  •  "V-  ■  lCd2  cpd  Tod)/(Cb2  CPP  V 

Expressing  the  Crocco  number  in  terms  of  pressure  gives 

M  -  1  -  <ypop)<v,)/1'p 
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Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 
and 


C 


d 


2 


<VPod> 


(Yj“1 )/y. 


where  Pod  is  the  total  pressure  on  the  "d"  streamline.  However,  at  the 
shear  layer  attachment  point,  the  l:d“  streamline  is  stagnated  after  passing 
through  an  oblique  shock.  The  maximum  strength  of  this  shock  corresponds  to 
the  deflection  angle  experienced  by  the  plume  boundary  when  turned  abruptly  by 
the  nozzle  wal 1 ,  i .e. , 


Pod/Pb 


where  P/Pb  is  obtained  from  the  oblique  shock  relations.  On  the  other  hand, 
if  the  shock  imparts  no  loss  in  total  pressure  to  the  "d"  streamline,  then  the 
flow  is  compressed  i sentropical ly  and 


pod  =  pb 

It  is  assumed  that  the  true  recompression  criterion  falls  between  these  two 
extremes,  so  that 


pod  =  $  +  (1-N)  Pb 


and 


0  <  N  <  1. 


The  quantity,  N,  is  known  as  the  Nash  factor  (see  Ref.  11,  12).  In  the 
analysis  presented  here,  N  is  used  to  calibrate  the  shear  layer  model  with 
experimental  data. 


Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 


Total  enthalpy  across  the  mixing  layer  is  also  assumed  to 
vary  linearly  with  respect  to  velocity  ratio,  i.e.,  at  the  "d"  streamline 


(’■  i) 


so  that 


C  T 
PP  op 


C  T 
pr  r 

C  T 
PP  op 


+  /  1  - 


op  J 


From  the  definition  of  the  Crocco  number,  and  since  Tr  =  Tocj ,  it  follows 


C  T  /  C  T  \ 

*('-  "d 

PP  op  \  PP  op  J 


Solving  the  above  equations  gives: 


^d  2 


A  .  cg r  .V  \ 

Cb  V  cpp  Top  ) 


Cpr  Tr 
Cpp  Top 


The  expressions  given  above  are  sufficient  to  determine  Td,  cd2,  and 
0d2  by  means  of  numerical  iteration,  provided  that  conditions  in  the  base 
recirculation  region  are  known.  The  similarity  coordinate  of  the  "d"  stream¬ 
line  is  then  given  by 

n,  =  erf"1  (2  0.  -1 ) 


#.  A*  *  »  "■ 


A 


III,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 


The  similarity  coordinate  for  the  "j"  streamline,  n.  can  be  found  from  the 
mass  and  momentum  conservation  equations.  Using  these  equations,  it  can  be 
shown  that  for  a  two-dimensional  shear  layer* 


(p/p  )  0dn  = 


S  (0 


/p  )  0  (1-0)  dn 


The  term  (p/p^),  which  is  the  ratio  of  the  gas  density  across  the  shear  layer 
to  the  density  at  the  primary  edge,  Pb>  can  be  evaluated  as  follows: 

(p/p^)  =  (p/pb) 


From  the  perfect  gas  law,  and  since  the  pressure  is  constant  across  the  shear 
layer. 


(c/0  )  =  Tb  _  Rp  Tb/Top  'op 

p/pb'  R  T  R  T/T  T 


S  i  nee 


T/Tn  =  1  -  c2 


Rn  !-Ck 

(p/ob}  =  / 


*For  example,  see  Bauer  (Ref.  7)  and  neglect  the  axisymmetric  terms. 


32 


Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 

As  shown  previously,  it  follows  directly  from  the  energy  equation  that 
2 

(6  = 


2 

L  C  T 


and  is  assumed  to  vary  linearly  with  velocity  across  the  shear  layer  as  follows 
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Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 


where 


C  =  C  (1-0)  +  0  C 
P  pr  PP 


R  =  Rr  0-0)  +  0  Rp 


The  similarity  coordinate,  n,  is  then  found  from  the  following  relationship: 


nj 


C 

cpp 


R 


r  2  h2 
Cb  * 


-3 


R^  0  0-0) 

c  R 
PP 


A  -  Cb  2  02 


da 


The  far  edges  of  the  shear  layer  are  taken  as  +3  and  -3.  The  value  for  nj 
is  obtained  by  integrating  the  left-hand  side  of  the  above  expression  until  it 
is  equal  in  value  to  the  right-hand  side.  For  computational  convenience,  the 
followng  functions  are  defined; 


A 


i,  («  -  j 


I,  (»)  * 


so  that 

’’j 


c  R 

X  -£ 

1  t  R 

3  or 


-  Cu 


2  .2 


dn  = 


c 

-3 

pp 

A 

c 

\ 

c 

) 

PP 

-3 

3 

C  R 

"5 

J>_  -  p 

C  R 

-3 

PP 

R 

R 


4. 

R 


0 

A  -  Cb2  02 


0 

A  -  cb  0 


dn 


dn 


d*>- 


■'  '  Cb  ?  ^ 


>  :  R 

S  -P  _E 

)  C  R 

-3  op 


.  'T7 

a  -  nb  «t 


dn 


J1  (nj)  =  Ij  (3)  -  I3  (3) 


Next,  a  method  is  presented  for  determining  the  species  concentrations  in  the 
base  recirculation  region.  The  shear  layer  to  the  left  of  the  "d"  streamline 
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III,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 

is  turned  into  the  recirculation  region  at  a  height  rw  above  the  primary 
flow  centerline.  By  applying  the  mass  and  energy  balance  equations  to  a 
control  volume  surrounding  the  recirculation  region,  it  can  be  shown  that 
(Ref.  13)  equation  [18]) 

■^d 

“s  <Wd  -  V>  ■  £  2'  rw 0  u(v  -  y  dy 


The  bleed  flow  is  entirely  made  up  of  the  secondary  flow  species,  so  that 


YPBleed  °p/pbleed  0 


Thus, 


0  *  V  ‘  ^  rd 

«s  ) 


PU  <Ypr  -  Yp>  dy 


For  unity,  Schmidt  and  Lewis  No. 


Y  -  Y  =  (1  -  Y  ) 

p  pr  v  pr' 


So  that 


2  nr. 


0  =  Y  + 
pr 


$ 


>U0dy  ( 1  -Ypr ) 


Changing  to  nondimensional  coordinates,  and  proceeding  as  before,  gives 


0  =  Y 


pr 


+  2”rw  ^pbUb  ^ 


{1-cb )  [i3  ^  n-vpr)] 
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III,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 

It  can  be  shown  that  the  mass  flow  in  the  secondary  stream  is* 

“s  ■ 

Combining  the  above  two  equations  gives 


where 


"S  o-cb2)  [Ii  ("j’  -  !i  (nd>] 


a  =  Y  /Y 

pr  sr 


The  total  enthalpy  in  the  base  recirculation  region  is 


therefore 


Cps  Tos 

r  T  . 

PP  op 


Cpr  Tr 
r — T — 
PP  op 


CPr  Tr 
CPP  Top 


-1 


(2) 


The  total  mass  flow  from  the  primary  nozzle  is 


-p..%  v 


i~vvop" 


(ir) 


V1 


Thus,  the  ratio  of  the  secondary  mass  flowrate  to  the  primary  mass  flowrate 


VWn 
s  p 


2 

o 


er 


w 


YP  +  1 

2^) 


[I1  C^j)  -  I-j  (nd)]  (3) 

•for  example,  see  Bauer  (Ref.  7)  and  neglect  the  axisymmetric  terms. 
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HI,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 

Equations  (1),  (2),  and  (3)  above  form  a  system  of  three  nonlinear  simultan¬ 
eous  equations  in  three  unknowns:  Vpr,  Pb/P0p»  and  Tr/T0p.  The 
aerodynamic  bleed  flow  computer  program  numerically  solves  for  these  three 
unknowns,  using  input  initial  estimates  of  their  values. 

2.  Aerodynamic  Bleed  Flow  Model  Computer  Program 

The  small  bleed  flow  model  for  the  shear  layer  derived  at  in 
the  previous  subsection  is  used  in  the  aerodynamic  computer  bleed  flow  pro¬ 
gram.  A  flow  chart  for  this  computer  program  is  presented  in  Figure  III -4 . 

The  procedures  shown  in  Figure  III -4  are  discussed  below. 

At  the  beginning  of  the  program,  initial  and  default  values 
are  set  and  the  input  data  is  read,  using  FORTRAN  NAMELIST.  (See  Appendix  A 
for  a  description  of  the  program  input  and  a  sample  problem.)  The  primary  and 
secondary  wall  geometries  are  then  constructed.  Both  the  primary  and  the 
secondary  nozzles  are  assumed  to  have  geometries  that  can  be  described  ana¬ 
lytically  by  the  following  parameters: 

Chamber  contraction  ratio 

Circular  wall  radius  connecting  the  chamber  and  the  nozzle 

i  nl  et 

Conical  nozzle  inlet 

Circular  wall  radius  upstream  of  the  nozzle  throat 

Circular  wall  radius  downstream  of  the  nozzle  throat 

The  nozzle  wall  geometries  are  defined  by  simple  geometrical  inputs  such  as 
wall  slope,  throat  radius  of  curvature,  etc.  Complete  shapes  are  specified  so 
that  the  wall  contours  can  be  plotted  as  well  as  used  in  the  analysis. 
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Figure  I II -4.  Flow  Chart  for  the  Aerodynamic  Bleed  Flow  Model  Computer  Program 


Ill,  B,  Aerodynamic  Bleed  Mow  Model  (cont.) 


In  order  to  be  used  easily  by  the  computer  programs,  the  wall  geometries  are 
converted  to  finely  spaced  tables  of  points. 


Next,  an  iterative  procedure  is  initiated  which  consists  of 
the  steps  enclosed  in  the  dashed  line  rectangle  shown  in  Figure  1 1 1-4.  The 
steps  taken  in  this  procedure  are  individually  discussed  below. 


(1) 


Conditions  at  the  exit  of  the  primary  nozzle  are  deter¬ 
mined  in  the  MAIN  subroutine  by  using  the  one-dimensional 
relation  between  Mach  number  and  area  ratio,  i.e.. 


=  M 
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(2)  The  exhaust  flow  is  turned  at  the  nozzle  lip  and  expands 
to  equal  the  boundary  Mach  number  M^.  The  angle  of  this 
turn,  is  determined  in  the  PLUME  subroutine  by  using 
the  Prandtl -Meyer  relations  for  supersonic  expansion 
at  a  point ,  as  f o 1  lows  : 
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1 1 1 , B ,  Aerodynamic  Bleed  Flow  Model  (cont.) 

(3)  The  plume  boundary  is  determined  by  using  the  modified 
method  of  Herron  (Ref.  10).  Since  the  plume  boundary  is 
a  circular  arc,  all  of  its  geometric  properties  (such  as 
its  position,  arc  length  to  any  point,  and  slope  at  any 
point)  are  known  analytically.  This  step  is 
accomplished  by  using  subroutine  PLUME. 

(4)  The  intersection  of  the  plume  with  the  secondary  nozzle 
wall  is  determined  geometrically  by  using  subroutine 
INTERS. 

(5)  The  shear  layer  model  allows  the  shear  layer  profile  to 
be  determined  exactly  at  the  point  of  intersection  of 
the  plume  boundary  and  the  wall.  All  properties  of  the 
profile  can  be  evaluated  from  various  integrals  of  the 
error  function,  as  described  in  Section  1 1 1  ,B ,  1 . 

These  integrals  are  evaluated  in  subroutines  IONEV, 

ITHRV ,  I1INVR,  and  IZEROS. 

(6)  The  plume  boundary  is  a  streamline.  The  streamline  will 
be  turned  abruptly  (i.e.,  deflected)  to  follow  the 
secondary  nozzle  wall,  resulting  in  a  shock  structure. 

The  cases  of  most  interest  correspond  to  solutions  of 
the  shock  relations  of  the  weak  family.  Cases  where  a 
lambda  shock  structure  develops,  so  that  the  flow  is 
shocked  subsonically  in  part  of  the  region,  are  of  less 
interest  since  they  correspond  to  aerodynamically  poor 
designs.  The  obliquie  shock  relations  are  used  to  deter¬ 
mine  the  rise  in  static  pressure  across  the  shock  and 
the  loss  of  total  pressure  across  the  shock.  Since  the 
development  of  shock  structure  in  the  nozzle  is  funda¬ 
mentally  a  two-dimensional  phenomenon,  no  attempt  is 
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ir.adt?  to  treat  the  1 1  uw  in  the  nozzle  downstream  of  the 
origin  ot  the  shoot.  The  oblique  shock  relations  are 
f  a1-  d  v.  hr  '  i  ne  WAVANG. 

The  cal  cut  at  iw'.  i-s>  pitied  above  require  a  first  estimate  for 
the  variables  IV,  T,-,  a:  d  •  p  •  I ‘  ;s  necessary  to  find  converged  values 
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flow  given  m  an  input  table,  converged  values  for  a  solution  are  automatic¬ 
ally  used  as  first  estimates  for  the  next  case. 
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III,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 

multiplier  for  the  similarity  parameter,  M^/y,  in  the  Herron  plume  correla¬ 
tion.  (M5  is  the  Mach  number  along  the  plume  boundary,  and  y  is  the  ratio 
of  specific  heats.)  Altering  the  plume  scaling  influences  the  plume  impinge¬ 
ment  angle  and,  hence,  the  resulting  shock  strength. 

A  plume  scaling  factor  of  1.2  resulted  in  model  predictions 
that  agreed  well  with  the  experimental  data  for  small  nozzle  spacing  as  shown 
in  Figure  II 1-5.  However,  for  tests  with  larger  nozzle  spacing,  the  model 
predictions  became  less  accurate  when  a  constant  scaling  factor,  as  shown  in 
Figure  1 1 1 -6 ,  was  used.  Decreasing  the  plume  scaling  factor  to  0.9  improved 
the  accuracy  for  large  nozzle  spacing  (Figure  1 1 1 -7 ) .  A  correlation  between 
plume  scaling  factor  and  nozzle  spacing,  as  shown  in  Figure  III -8 ,  can  be 
inferred,  but  the  validity  of  applying  it  to  other  dual  throat  geometries 
needs  to  be  determined. 

The  second  parameter,  the  Nash  factor  (N),  is  used  in  the 
recompression  criteria  relating  the  pressure  behind  the  plume  attachment  shock 
to  the  recirculation  pressure.  This  parameter  adjusted  the  stagnation  pres¬ 
sure  loss  to  match  the  cold-flow  data.  While  the  plume  scaling  factor  shifts 
the  base  pressure  versus  bleed  flow  curve  up  or  down,  the  Nash  factor  changes 
the  actual  shape  of  the  curve  (Figure  1 1 1-9).  A  Nash  factor  of  0.4  was  found 
to  result  in  model  predictions  that  agreed  well  with  the  cold-flow  data. 

As  shown  in  Figure  1 1 1- 5 ,  the  bleed  flowrate  for  the  "blow- 
off"  condition  matches  the  cold-flow  data  to  within  1/2%,  but  the  base  pres¬ 
sure  prediction  becomes  asymptotic  after  3.5%  bleed.  At  3.5%,  the  plume 
impingement  results  in  a  3%  loss  in  stagnation  pressure  according  to  the 
oblique  shock  calculations.  The  Nash  "recompression  criteria"  indicate  the 
actual  stagnation  pressure  loss  to  be  approximately  0.4  of  the  3%.  So  the 
situation  that  exists  is  that  increasing  the  bleed  flowrate  from  3.5%  to  5.5% 
results  in  a  small  change  in  stagnation  pressure  loss  and  in  the  corresponding 
base  pressure. 
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Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 


The  third  empirical  factor  is  the  jet  spreading  parameter 
(sigma),  which  is  used  in  the  shear  layer  calculations  and  defines  the  rate  at 
which  the  shear  layer  jets  spread.  Presently,  the  bleed  flow  model  uses  a 
correlation  developed  by  Korst  which  states  that  the  spreading  parameter  is  a 
function  of  plume  boundary  Mach  number.  From  the  derivation  of  the  shear 
layer  calculations  presented  in  Section  III,B,1,  it  can  be  seen  that  the  jet 
spreading  parameter  will  influence  the  mass  fraction  Ypr  and  recirculation 
temperature  Tr/T0.  Calibration  to  the  cold-flow  data  was  not  possible 
since  the  recirculation  temperature  was  not  measured  during  cold  or  hot-fire 
testi ng. 


4.  Parametric  Analysis 

A  parametric  analysis  was  performed  to  investigate  the 
effects  of  varying  bleed  flow  thermodynamic  properties  and  primary  nozzle 
geometry  on  bleed  flow  requirements.  Changes  were  made  in  molecular  weight 
(MW)  and  specific  heat  (Cp)  for  both  primary  and  bleed  flows,  bleed  flow  to 
primary  flow  stagnation  temperature  ratio,  nozzle  axial  spacing,  primary 
nozzle  area  ratio,  primary  nozzle  size  and  contour. 

Table  1 1 1  - 1 1  lists  the  effects  of  bleed  flow  thermodynamic 
properties  on  bleed  flow  requirements.  The  nozzle  configuration  simulated  was 
that  used  for  the  dual  throat  cold-flow  tests  (NAS  8-32666).  Varying  the 
ratio  of  specific  heats  (yD)  f°r  the  bleed  flow  appears  to  have  little 

D 

effect  on  the  required  bleed  flowrate.  Variations  in  bleed  flow  molecular 
weight  and  stagnation  temperature  appear  to  have  significant  effects  on  bleed 
flow  requirements.  Decreasing  the  molecular  weight  of  the  bleed  gas  results 
in  a  lower  required  flowrate,  whereas  a  decrease  in  stagnation  temperature 
results  in  an  increase  in  the  required  bleed  flow.  An  increase  in  the  primary 
flow  gamma  (y  )  will  also  result  in  requiring  less  bleed  flow  to  turn  the 
plume  for  a  shock-free  attachment. 
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III,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 

The  effect  of  geometric  design  parameters  on  bleed  flow 
requirements  was  investigated  by  single  parameter  variations  about  a  baseline 
configuration.  A  baseline  conf iguration  was  generated  from  the  operating  spe¬ 
cifications  presented  in  the  final  report  of  the  Dual-Fuel,  Dual  Throat  Engine 
Preliminary  Analysis  (NAS  8-32967).  Table  1 1 1  - 1 1 1  lists  the  main  geometric 
design  parameters  for  this  design,  and  Figure  1 1 1 -10  shows  the  baseline 
geometry  and  plume  contours  for  various  bleed  flowrates.  For  the  baseline 
design,  the  bleed  flow  model  calculated  that  a  5.94%  bleed  flowrate  was  needed 
for  the  plume  to  attach  shock-free  to  the  secondary  nozzle,  i.e.,  the  blowoff 
condition.  Reducing  the  bleed  flow  from  the  blowoff  condition  to  approxi¬ 
mately  3%  showed  that  the  resulting  plume  boundary  pressure,  plume  contour, 
and  shock  strength  were  practically  the  same  in  the  blowoff  condition. 

The  results  of  varying  the  nozzle  spacing  are  presented  in 
Table  III-IV.  Figure  1 1 1 - 1 1  shows  the  configurations  used  and  the  resulting 
plume  contours  for  0%  bleed  flow.  For  this  case,  the  plume  attachment  point 
moves  further  upstream  as  the  nozzle  spacing  increases.  The  plume  boundary 
pressure  also  increases  with  increased  nozzle  spacing,  resulting  in  a 
decreased  slope  of  the  plume  at  the  attachment  point.  Further,  as  the  nozzle 
spacing  increases,  the  bleed  flow  requirement  for  flow  attachment  at  the 
secondary  throat  also  increases.  The  net  effect  of  these  trends  is  that  the 
total  flow  angle  decreases  with  increased  nozzle  spacing,  resulting  in  reduced 
shock  strength. 


Table  III-V  lists  the  bleed  flow  requirements  for  varying  the 
primary  nozzle  area  ratio.  As  the  primary  nozzle  area  ratio  increases,  the 
bleed  flow  requirement  decreases.  For  0%  bleed  flow,  increasing  the  primary 
nozzle  area  ratio  causes  the  plume  attachment  point  to  move  upstream  and  the 
plume  boundary  pressure  to  decrease.  The  decrease  in  plume  boundary  pressure 
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ABLE  1 1 1- II I 


DUAL  TIPOAT  BASELINE  GEOMETRY 

Mode  T  Sea  Level  Thrust  -  600,000  IbF 
Mode  I  to  Mode  II  Thrust  Ratio  -  3.2 

SECONDARY  NP7ZLE: 

Chamber  Pressure 
Mixture  Ratio 
Area  Ratio 
Contraction  Ratio 
Throat  Radius 
Percent  Mode  I  Thrust 
Inlet  Radius 
Upstream  Throat  Radius 
Downstream  Throat  Padius 
Nozzle  Length 

PRIMARY  N07ZLE : 


Chamber  Pressure 

4000  psia 

Mixture  Ratio 

7.0 

Area  Ratio 

2.7 

Contraction  Ratio 

2.5 

Throat  Padius 

2.90  in. 

Percent  Mode  I  Thrj.< 

30% 

Inl^t  Radius 

2.90  in. 

Upstream  T^roa*  Pad’ us 

2.90  in. 

Downstream  Throat  Radius 

1.80  in . 

Nozzle  Lenoth 

6.3  in. 

Nozzle  'paring 

2.75  in. 

2800  psia 
2.8 
45.6 
2.5 

6.41  in. 
70% 

6.41  in. 
6.41  in. 
3.25  in. 
123.7  in. 
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Throat  Baseline  Configuration 


TABLE  1 1 1 -IV 


variation 
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.07170 
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- 
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3.  5 

0 
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15.4 
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-22.4 

3.0 

.07793 

2.29 

9.8 

14.7 

-4.9 

8.18 

.08215 

- 

- 

- 

- 

4.0 

0 

.06239 

2.41 

13.4 

36.8 

-23.4 

3.0 

.08187 

2.27 

9.0 

17.0 

-8.0 

★ 

8.0 

.08765 

2.23 

7.9 

2.4 

5.5 

* 

9.76 
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2.23 
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Vpo  " 

pe/pc  Validity  of  Results  Questionabl 
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M 

rP 

2.25 
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O 

o 

o 

psia 

P  /P 
c  c 

.08389 

Rtp 

2.90 

’■p 

2.707 

F  Ratio 

=  3.2 

XLE  -  Spacing  Between  Primary  Nozzle  Exit  and  Secondary  Nozzle  Throat 

WS/WD  -  Bleed  flowrate 

Pg/Pp  -  Plume  Boundary  to  Chamber  Pressure  Ratio 

Mp  -  Mach  Number  Alonq  Plume  Boundary 

‘Plump  "  Plum?  Angle  at  Attachment  Point 

sri  -  Total  Flow  Turnina  Angle 

Flow 

1 |  -  Wall  Analp  At  Attachment  Point 

'.g  ^  -  Shear  Layer  Thickness 

f.  -  Primary  Nozzle  Area  Ratio 

P 

R  -  Primary  Nozzle  Throat  Radius 

E  Ratio  -  Mode  I  to  Mode  II  Thrust  Ratio 

Pe/PQ  -  Pressure  Ratio  for  Primary  Nozzle 


6S.L. 

.10763 

.15506 

.21629 

.22967 

.20673 

.25875 

.24294 

.28907 

.32214 
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III,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 

is  due  to  the  decreasing  nozzle  exit  pressure  from  the  increase  in  nozzle  area 
ratio.  Figure  1 1 1 -12  shows  the  resulting  plume  contours  for  0%  bleed  flow. 

Changing  the  primary  nozzle  size  was  accomplished  by  keeping 
the  same  primary  nozzle  area  ratio  but  changing  the  throat  area.  A  change  in 
throat  area  corresponds  to  a  change  in  the  Mode  i  to  Mode  II  thrust  ratio. 
Table  1 1 1 -V I  lists  the  change  in  bleed  flow  requi reinents  for  change  in  throat 
radius.  As  the  throat  radius  increases,  decreasing  thrust  ratio,  the  bleed 
flow  requirement  increases  slightly.  For  the  0%  bleed  flow  case,  the  plume 
attachment  point  moves  upstream,  and  the  plume  boundary  pressure  increases 
with  increasing  throat  radius.  Figure  1 11-13  shows  the  plume  contours  for  the 
0%  bleed  flow  case. 

The  effect  of  primary  expansion  nozzle  contour  was  evaluated 
by  comparing  aerodynamic  model  predictions  for  the  baseline  Bell  contour  noz¬ 
zle  with  an  exit  angle  of  16°  to  results  calculated  for  a  10°  conical  nozzle. 
The  low  exit  angle  conical  nozzle  is  predicted  to  require  less  bleed  flow 
(4.5%  vs  5.9%)  which  subsequently  results  in  a  lower  base  pressure  ratio 
(0.0717  to  0.05489).  These  decreased  values  are  due  to  the  change  of  nozzle 
exit  angle  which,  when  combined  with  the  Prandtl -Meyer  expansion  angle  at  the 
nozzle  lip,  results  in  a  smaller  plume  expansion  angle.  Having  a  plume  bound¬ 
ary  with  a  shallow  angle  requires  less  bleed  flow  and  a  smaller  base  pressure 
to  turn  the  plume  and  allow  it  to  attach  tangentially  to  the  secondary  nozzle 
wal  1 . 


In  summary,  the  dual  throat  aerodynamic  bleed  flow  model 
developed  during  NAS  8-32666  has  been  updated  to  include  a  procedure  for  pre¬ 
dicting  the  base  pressure  generated  in  Mode  II  operation  with  zero  bleed  flow 
and  provide  for  dissimilar  exhaust  gases  to  be  used  for  the  primary  and  bleed 
flows.  These  changes  have  enhanced  the  analysis  capability  of  the  Mode  II 
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Figure  1 1 1 - 1 2 .  Effect  of  Varying  Primary  Nozzle  Area  Ratio  on  Plume  Attachment  for  0%  W  /W 


TABLE  1 1 1 -VI 
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Ill,  B,  Aerodynamic  Bleed  Flow  Model  (cont.) 


dual  throat  operation  and  provide  a  design  tool  to  permit  tradeoffs  of  various 
geometric  and  operating  parameters. 

The  improved  model  has  been  used  to  evaluate  bleed  flow 
requirements  as  a  function  of  both  bleed  flow  properties  and  geometric  design 
parameters.  These  results  are  summarized  in  Table  1II-VII.  This  model, 
together  with  other  dual-throat  design  and  analysis  improvements,  has  also 
been  utilized  to  assess  changes  in  predicted  performance  of  the  dual-throat 
concept  as  discussed  in  later  sections  of  this  report.  The  bleed  flow  model, 
with  suitable  modifications,  also  provided  the  basis  for  prediction  of  the 
base  flow  region  of  the  dual  expander  concept  as  described  in  Section  V  of 
this  report. 


C.  MODE  II  BOUNDARY  LAYER  LOSS  VERIFICATION 
1.  Previous  Simplified  Approach 

During  Contract  NAS  8-32666,  the  Mode  II  boundary  layer  loss 
was  obtained  by  using  the  turbulent  boundary  layer  chart  procedures  outlined 
in  Appendix  B  of  CP IA  No.  178  with  some  additional  considerations  to  the  aero¬ 
dynamics  associated  with  Mode  II  operation.  The  primary  assumptions  for  this 
calculation  are  that  the  boundary  thrust  decrement  is  additive  and  propor¬ 
tional  to  the  momentum  thickness.  The  boundary  layer  loss  was  obtained  by 
using  the  TBL  chart  procedures  for  a  conventional  Rao  nozzle  having  an  area 
ratio  defined  by  the  point  of  plume  attachment  and  the  primary  throat.  A 
ratio  of  the  shear  layer  momentum  thickness  (at  the  point  of  attachment  given 
by  the  aerodynamic  bleed  flow  program)  to  the  momentum  thickness  given  by  the 
TBL  chart  procedures  was  used  to  proportion  the  conventional  nozzle  boundary 
layer  loss  to  approximate  the  actual  thrust  layer  loss  obtained  for  the  nozzle 
from  the  attachment  point  to  the  secondary  nozzle  exit,  yielding  the  total 
Mode  II  boundary  layer  loss.  This  procedure  is  outlined  in  Figure  1 1 1-14. 
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Ill,  C,  Modi'  II  Boundary  Layer  loss  Verification  (cont.) 

These  calculations  were  made  by  assuming  an  adiabatic  wall  to  simplify  the 
calculation  of  loss  due  to  heat  extraction  from  the  boundary  as  well  as  the 
1 1 1 ’t  due  to  tiie  addition  of  tins  heat  to  the  incoming  regenerative  coolant 


?.  BLIMP  Analysis 

The  BLIMP  computer  program  was  used  to  provide  a  rigorous 
analysis  of  the  boundary  layer  loss  during  Mode  II  operation  for  comparison  to 
the  Amplified  approach  described  above.  BLIMP  was  successfully  run  using  the 
Kendal 1  and  binary  diffusion  options  for  the  Mode  II  dual  throat  baseline  con¬ 
figuration  with  full  bleed  flow.  These  options  are  consistent  with  the  stan¬ 
dardized  JANNAF  performance  methodology.  To  accomplish  this,  it  was  necessary 
to  modify  the  aerodynamic  bleed  flow  model  to  generate  additional  points  used 
in  the  spline  fit  to  describe  the  effective  Mode  II  "nozzle  contour"  and  to 
calculate  the  momentum  thickness  at  the  plume  attachment  point,  including  the 
contribution  due  to  the  bleed  flow.  Also,  the  TDK  program  was  modified  to 
accept  *  he  increased  spline  fit  ,  i  zc  needed  to  calculate  the  pressure  profile 
to  '.!•••  used  m  BLIMP. 


The  BLIMP  run  was  started  at  the  plume  attachment  point  (Sta- 
•  ion  1)  and  continued  throughout  the  secondary  nozzle.  The  BLIMP  input  param¬ 
eter  3(1),  which  represents  the  length  of  a  flat  plate  required  to  generate 
the  (Psi'od  momentum  thickness,  was  varied  until  the  momentum  thickness  at 
Station  1  emdled  the  momentum  thickness  calculated  by  the  bleed  flow  program. 
In  this  manner,  'he  conditions  for  the  first  calculation  station  within  the 
BLIMP  solution  are  representative  (at  least  in  terms  of  momentum  deficit)  of 
the  conditions  calculated  Tor  the  shear  layer  using  the  bleed  flow  model.  In 
'  f  tie  i  words,  t  tie  I’.Ll'T  analysis  was  started  with  initial  boundary  layer  pro¬ 
files  which  are  representat i ve  of  the  boundary  layer  growth  within  the  primary 
nozzle  and  shear  layer,  including  the  effects  of  bleed  flow  addition. 
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Ill,  C,  Mode  II  Boundary  Layer  Loss  Verification  (cont.) 


Starting  with  these  initial  conditions,  the  BLIMP  computer  program  calculated 
the  additional  boundary  layer  growth  in  the  secondary  nozzle  downstream  of  the 
plume  attachment  to  the  nozzle  exit.  Standard  BLIMP  parameters  at  the  final 
station  were  then  used  to  calculate  the  net  boundary  layer  loss,  as  shown  in 
Table  III-VIII. 


Table  III-VIII  also  shows  a  comparison  between  the  boundary 
layer  loss  calculated  from  the  BLIMP  analysis  and  that  determined  by  using  the 
simplified  procedure  from  NAS  8-32666.  As  can  be  seen,  the  net  loss  based  on 
the  BLIMP  results  is  more  than  twice  as  great  as  the  loss  calculated  by  using 
the  simplified  procedure.  Since  the  BLIMP  solution  is  by  far  the  more  rigor¬ 
ous  procedure,  it  is  considered  the  best  available  current  estimate  for  Mode 
II  boundary  layer  losses.  It  should  be  noted,  however,  that  even  the  net 
boundary  layer  loss  from  the  BLIMP  solution  includes  the  performance  loss 
resulting  from  the  addition  of  the  bleed  flow. 

As  a  result  of  the  poor  comparison  between  the  simplified  and 
rigorous  boundary  layer  loss  procedures,  a  new  simplified  procedure  was  devel¬ 
oped  in  order  to  provide  loss  estimates  close  to  those  calculated  from  the 
BLIMP  procedure  without  extensive  and  detailed  analysis.  While  the  revised 
procedure  still  uses  the  TBL  chart  as  a  basis  for  the  loss  calculation,  it 
scales  the  results  to  account  for  shear  layer  effects  in  a  different  manner 
than  the  previous  method. 

Briefly,  the  TBL  chart  method  is  used  to  calculate  the  momen¬ 
tum  thickness  for  a  conventional  nozzle  expansion  from  the  primary  nozzle  to 
the  plume  attachment  point.  The  momentum  thickness  calculated  from  the  shear 
model  is  then  divided  by  this  conventional  momentum  thickness  to  define  a 
parameter,  RTHETA,  which  is  used  in  the  simplified  procedure  to  account  for 
upstream  influences  on  the  boundary  layer.  Thus, 
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TABLE  1 1 ! -VI 1 1 


-CP.M_pARijO!LO-r  boundary  layer  loss  calculat i ons 


Mode  II  Baseline  Configuration  (See  Table  III-III) 

Bleed  Flow  6.3 


Boundary  Layer  Performance 
Pa  rameters 

Total  Flowrate  (lbm/sec) 

Gross  Thrust  Loss  (IbF) 

Gross  Isd  Loss  (sec) 

Total  Heat  Loss  (Btu/sec) 

Isp  Gain  From  Reqen  Heat  (sec) 

Net  RLL  (sec) 

' 1  Rased  on  TW/TS  -1.0  RTHETA  = 

:,'tised  on  TW/TS  --  0.?  RTHFTA  = 

^  nased  on  Difference  Between  Boundary 


BOUNDARY 

LAYER  ANALYSIS  TECHNIQUE 

BLIMP 

NAS  8-32666 
Simpl ified 
Procedure 

Revised 

Procedure 

483.3 

483.3 

483.3 

8095 

3335 

7974 

16.8 

6.9(1^ 

16.5^2) 

37,050 

0 

N/A 

3.3 

0 

3.4(3) 

13.5 

6.9 

13.1 

1.00 

2.22 

Layer  Loss  for  TW/TS  =  0.2  &  TW/TS  =  1.0 


'  w  -  Wall  Temperature 
f  Stagnation  Temperature 

s 


63 


Ill,  C,  Mode  II  Boundary  Layer  Loss  Verification  (cont.) 


RTHETA  = 


0shear/0TBL 
model  chart 


Plume  attach  point 


Previously  defined  trends,  such  as  those  shown  in  Figure  III-15,  are  used  to 
estimate  the  change  in  the  RTHETA  value  for  further  expansion  in  the  secondary 
nozzle  section.  This  final  RTHETA  value  (i.e.,  at  the  secondary  nozzle  exit) 
is  used  to  ratio  the  momentum  thickness  calculated  for  a  conventional  nozzle 
expansion  from  the  primary  nozzle  thioat  to  the  secondary  nozzle  exit,  i.e.. 


eExit 


0Exit,  TBL 


x  RTHETA 


The  boundary  layer  calculations  are  made  by  using  the  actual 
wall  temperature  conditions.  In  the  event  of  regenerative  cooling,  the  new 
approach  also  takes  into  account  the  increase  in  the  IspoDE  due  t0 
increased  propellant  enthalpy  transferred  to  the  regenerative  coolant. 


This  approach  was  incorporated  into  a  computer  program  called 
DT*BLL0SS.  The  inputs  required  for  the  program  are  1)  characteristic  velocity 
(Cp*).  2)  chamber  pressure  (PCp),  3)  throat  radius  (Rtp),  4)  ratio  of 
specific  heats  (yp).  3)  Mode  II  area  ratio  (£eny),  6)  ODE  Isp  at  £eng> 

7)  area  ratio  at  plume  attachment  point  (6-plm ) ,  and  8)  momentum  thickness 
(0Bieed)  the  plume  attachment  point.  The  first  six  parameters  are 
obtained  from  ALRC's  Dual  Throat  Geometry  and  Performance  Prediction  Program 
(FD  0169).  The  last  two  parameters  come  from  the  Aerodynamic  Bleed  Flow 
Program  (FD  0170). 

The  results  of  this  "revised  procedure"  for  the  baseline  case 
are  compared  to  both  the  BLIMP  rigorous  and  previous  simplified  cases  in  Table 
III-VIII.  As  noted,  the  new  procedure  provides  loss  predictions  which  are 
similar  to  the  BLIMP  calculations.  The  impact  of  this  new  procedure  on  the 
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performance  predictions  for  Mode  II  dual  throat  operation  is  identified  quan¬ 
titatively  in  Section  III,F  of  this  report. 

D.  SHOCK  MODEL 

The  objective  of  this  subtask  was  to  develop  a  method  for 
obtaining  two-dimensional  axially  symmetric  flow  field  solutions  that  include 
the  effects  of  shocks  to  estimate  shock  effects  on  performance  and  thermal 
operating  characteri sties.  Mode  II  operation  at  design  conditions  will  neces¬ 
sarily  contain  shock  waves,  the  effect  of  which  may  be  negligible  at  low  area 
ratios  but  which  are  significant  at  high  area  Mode  II  conditions. 

Since  shock  wave  development  in  a  nozzle  is  fundamentally  a  two- 
dimensional  problem,  a  two-dimensional  computer  analysis  must  be  used  to  pre¬ 
dict  the  effect  of  shocks  on  performance.  The  computer  code  selected  for  this 
purpose  is  the  VNAP  code  (Ref.  4).  It  has  been  modified  for  use  in  calcu¬ 
lating  both  the  flow  field  within  the  dual  throat  nozzle  during  Mode  II  opera¬ 
tion  and  the  resulting  performance. 

The  VNAP  program  was  selected  for  performing  the  dual  throat  noz¬ 
zle  calculations  because  it  is  the  only  computer  code  readily  available  that 
is  capable  of  treating  supersonic  nozzle  flow  with  imbedded  shocks  in  a  gen¬ 
eral  manner.  It  is  capable  of  calculating  viscous  as  well  as  inviscid  steady 
and  unsteady  internal  flow.  The  approach  used  is  to  solve  the  Navier-Stokes 
equations  for  two-dimensional ,  time-dependent,  compressible  flow.  The  fluid 
is  assumed  to  be  a  perfect  gas.  An  explicit  artificial  viscosity  is  used  for 
the  shock  calculations.  The  steady-state  solution  is  obtained  as  the  asympto¬ 
tic  solution  for  large  time.  The  nozzle  walls  are  transformed  to  a  rectangu¬ 
lar  grid  so  that  the  resulting  partial  differential  equations  can  be  solved  by 
using  the  second-order  MacCormack  finite-difference  scheme. 
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Ill,  D,  Shock  Model  (cont.) 


Therefore,  it  appears  that  the  VNAP  calculation  results  in  a  performance  which 
is  approximately  2%  greater  than  an  accurate  method  of  characteristic  solution 
(TDK)  and  that  the  .978  correction  factor  should  be  utilized  for  the  two  addi¬ 
tional  VNAP  cases  which  contain  flow  shocks. 

Plots  were  generated  showing  velocity,  Mach  number,  pressure,  tem¬ 
perature,  and  density  profiles  to  help  visualize  the  flow  field  situation. 
Examples  of  these  plots  for  the  conventional  nozzle  analysis  are  shown  in 
Figures  1 1 1 -16  through  1 11-21.  To  help  facilitate  the  computer  storage 
requirement,  the  nozzle's  flow  fields  were  arbitrarily  divided  into  two 
segments,  referred  to  as  upstream  and  downstream  on  the  referenced  figures. 

The  VNAP  simulation  of  the  dual  throat  baseline  Mode  II  configura¬ 
tion  predicted  a  nozzle  performance  of  480.2  seconds.  This  was  corrected  to 
469.6  seconds  by  using  the  correction  factor  obtained  from  the  conventional 
nozzle  simulation.  The  resulting  Mode  II  nozzle  divergence  efficiency  was 
found  to  be  0.988,  compared  to  a  conventional  nozzle  efficiency  of  .996. 

Thus,  the  divergence  efficiency  of  the  dual  throat  Mode  II  expansion  is 
approximately  0.8%  less  than  that  of  a  conventional  nozzle  expanding  to  the 
same  area  ratio.  This  is  only  0.3%  greater  than  the  value  used  during  NAS 
8-32666. 


Plots  of  the  flow  field  (Figures  III -22  to  III-27)  show  a  right 
running  shock  wave  forming  in  the  nozzle  due  to  the  compression  caused  by  the 
highly  contoured  secondary  nozzle  wall.  This  shock  structure,  however,  never 
reaches  the  downstream  secondary  wall  and  thus  has  only  a  slight  influence  on 
performance  and  heat  transfer. 

The  VNAP  program  provided  a  means  for  determining  the  gas-side 
boundary  conditions  used  in  the  thermal  analysis  of  the  secondary  throat 
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Figure  III-17. 
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Figure  I II -21 .  Conventional  Nozzle  Downstream  Pressure  Plot 
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Figure  II 1-22.  Dual  Throat  Baseline  Upstream  Velocity  Vector  Plot 
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-25.  Dual  Throat  Basel ine  Downstream  Mach  Number  Plot 
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Figure  IH-26.  Dual  Throat  Baseline  Upstream  Pressure  Plot 


Figure  111-27.  Dual  Throat  Baseline  Downstream  Pressure  Plot 


Ill,  D,  Shock  Model  (cont.) 


reyion  for  the  case  having  zero  bleed  flow  (the  case  expected  to  be  the  most 
adverse  Mode  II  operating  condition).  Pressures  and  mass  velocities  were 
obtained  at  different  axial  locations  and  used  as  boundary  conditions  for  the 
design  thermal  analysis  (discussed  in  Section  IV). 

The  flow  field  plots  (Figures  111-28  through  III -30 )  show  the  for¬ 
mation  of  the  plume  attachment  shock.  The  plots  are  indicative  of  a  weak 
shock  formation.  This  can  also  be  seen  in  the  wall  pressure  profile  shown  in 
Figure  1 11-31.  The  static  pressure  rise  is  approximately  15%  larger  than  the 
value  calculated  by  the  aerodynamic  bleed  flow  model,  but  it  is  approximately 
45%  of  what  would  be  expected  with  a  discrete  shock.  This  value  appears  to  be 
consistent  with  the  Nash  factor  value  of  0.4  chosen  during  calibration  of  the 
bleed  flow  model,  as  described  in  Section  III,B,2. 

t.  SECONDARY  NOZZLE  CONTOUR  OPTIMIZATION 

During  Mode  I  operation,  the  dual  throat  nozzle  operates  as  a  con¬ 
ventional  nozzle.  The  secondary  nozzle  contour  is  designed  by  using  the  opti¬ 
mization  t'chnique  developed  by  G.  V.  R.  Rao.  This  technique,  which  is 
derived  from  the  calculus  of  variations  and  employs  the  method  of  characteris¬ 
tics,  is  used  in  ALRC's  nozzle  design  program  (FD  0002)  and  will  calculate 
contours  to  provide  maximum  thrust  for  a  specified  length  and/or  area  ratio. 
The  Rao  method  is  applicable  when,  together  with  an  initial  wall  contour,  a 
shock-free  supersonic  data  line  is  given  across  the  nozzle.  In  ALRC's  pro¬ 
gram,  this  initial  data  line  is  automatical ly  provided  by  using  the  modified 
Hall  method  to  approximate  the  solution  of  the  transonic  flow  equations  in  the 
nozzle  throat  region.  In  general,  however,  the  method  can  be  applied  to  any 
supersonic  approach  flow  that  is  shock-free.  This  fact  has  been  utilized  in 
the  study  presented  here  to  provide  the  design  computer  program  with  the 
capability  of  designing  optimum  nozzle  contours  for  dual  throat  thrusters. 

The  method  used  to  desiyn  nozzles  for  Mode  II  operation  is  described  below. 
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Figure  m-28.  Dual  Throat  Baseline  Zero  Bleed  Flow  Velocity  Vector  Plot 
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yure  1 1 1-29.  Dual  Throat  Baseline  Zero  Bleed  Flow  Pressure  Profile 
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Figure  III-30.  Dual  Throat  Baseline  Zero  Bleed  Flow  Mach  Number  Profile 
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Figure  III-31.  Dual  Throat  Baseline  Zero  Bleed  Flow  Wall  Pressure  Profile 


Ill,  E,  Secondary  Nozzle  Contour  Optimization  (cont.) 

1.  Design  Method 

The  aerodynamic  model  presented  in  Section  III,B  can  be  used 
to  design  a  dual-throat  thruster  assembly  which,  by  use  of  bleed  flow,  oper¬ 
ates  with  a  plume  shape  and  attachment  position  such  that  the  flow  field  is 
nearly  shock-free.  Once  the  approach  flow  is  specified  up  to  the  plume 
attachment  position,  the  Rao  procedure  can  be  used  to  calculate  an  optimal 
contour  for  the  exterior  wall  downstream  of  the  attachment  position.  This 
wall  contour  is  typically  quite  long  because  the  thruster  mission  requires  a 
high  expansion  ratio. 

The  design  computer  program  has  been  modified  to  allow  an 
optimum  nozzle  contour  to  be  calculated  where  the  approach  flow  is  determined 
by  an  input  nozzle  contour.  This  feature  of  the  program  is  called  the 
"Extended  Nozzle  Wall  Option."  The  steps  used  to  implement  this  option  are  as 
follows: 


First,  the  wall  geometry  subroutines  of  the  TDK  method  of 
characteristics  program  were  transferred  to  the  ALRC  program  and  checked  out. 
These  are  subroutines  WALL,  CUBIC,  and  XSLP.  The  maximum  number  of  points 
allowed  for  inputting  the  nozzle  wall  using  the  spline  fit  option  (IWALL  =  4) 
was  increased  from  20  to  40.  This  same  modification  was  also  made  to  the 
aerodynamic  model  computer  program. 

Next,  subroutines  WALL1  and  WALL2  were  revised  to  allow  the 
circular  arc  portion  of  the  wall  contour  to  be  centered  at  a  predetermined 
position.  Previously  it  had  been  necessary  to  position  the  circular  arc  at 
the  nozzle  throat  and  center  it  at  (1,  0).  A  new  wall  point  MOC  subroutine, 
called  WALLPT,  was  also  written  and  checked  out. 
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Ill,  E,  Secondary  Nozzle  Contour  Optimization  (cont.) 

A  logic  subroutine,  called  RRC,  was  written  to  construct  a 
nozzle  flow  field  for  a  given  nozzle  wall  geometry  (i.e.,  as  constructed  by 
subroutine  WALL)  using  Right  Running  Characteristics  (RRC * s ) -  The  flow  field 
constructed  (shown  in  Figure  III -32 ) ,  is  bounded  as  follows: 

(a)  The  upstream  boundary  is  the  initial  data  line  RRC. 

(b)  The  upper  boundary  is  the  nozzle  wall  contour. 

(c)  The  lower  boundary  is  the  flow  axis. 

(d)  The  downstream  boundary  is  the  RRC  extending  from  the 
wall  end  point. 

Finally,  the  main  program,  FD002,  and  the  logic  program, 
LOGIC,  were  modified  to  accept  the  above.  The  program  was  checked  out  for  the 
baseline  case  for  the  dual  throat  nozzle.  That  is,  an  optimum  nozzle  contour 
of  fixed  length  and  area  ratio  was  designed  for  the  secondary  nozzle  wall, 
assuming  the  approach  flow  to  be  that  given  by  the  primary  nozzle  and  plume 
(aerodynamic  wall)  with  5.3%  bleed  flow. 

2.  Optimization  Analysis  Results 

A  parametric  analysis  was  performed  to  investigate  the  effect 
of  optimizing  the  secondary  nozzle  contour  for  either  Mode  I  or  Mode  II. 

Table  III-IX  and  Figure  1 1 1 -33  show  the  results  of  this  analysis,  including  a 
comparison  of  the  Mode  II  nozzle  contours  according  to  the  following  param¬ 
eters: 


-32.  Secondary  Nozzle  Optimization  Program  Procedure 


Ill,  E,  Secondary  Nozzle  Contour  Optimization  (cont.) 

Contour  (1)  is  a  conventional  Bell  nozzle  with  an  expansion 
ratio  of  222:1.  The  chamBer  throat  and  exit 
point  of  this  nozzle  are  those  of  the  primary 
nozzle  of  the  dual  throat  baseline  contour. 

Contour  (2)  results  from  optimizing  the  secondary  nozzle  for 
Mode  I  operation.  The  plume  attaches  to  a 
secondary  nozzle  contour  designed  for  an  46:1 
area  ratio  n  Mode  I  operation. 

Contour  (3)  has  been  designed  in  t»rms  of  "Extended  Nozzle 
Wall  Option"  discussed  in  this  section.  The 
approach  flow  passes  through  the  dual  throat 
primary  nozzle  and  the  plume  expansion,  subse¬ 
quently  attaching  at  the  secondary  nozzle 
throat.  Contour  (3)  is  the  optimal  contour 
given  this  approach  flow.  This  nozzle  has  an 
expansion  ratio  of  222:1. 

Contours  (1)  and  (3)  have  the  same  throat  and  exit  point, 
but,  as  can  be  seen  from  Figure  1 1 1-33,  they  are  very  different  in  shape. 
Contours  (2)  and  (3)  have  different  throat  locations  but  the  same  exit  point. 
Contour  (2)  is  designed  to  yield  maximum  thrust  during  Mode  I  operation. 
Contour  (3)  is  designed  to  yield  maximum  thrust  during  Mode  II  operation.  As 
can  be  seen  from  Figure  1 1 1-33 ,  Contour  (3)  has  less  curvature  than  Contour 
(2)  in  order  to  turn  the  high  Mach  number  approach  flow  isentropically.  The 
net  performance  effect  of  the  two  optimization  options  is  shown  in  Table 
III-IX. 


Throat  Secondary  Nozzle  Contour  Comparison 


Ill,  E,  Secondary  Nozzle  Contour  Optimization  (cont.) 


TABLE  III-IX 

COMPARISON  OF  NOZZLE  DIVERGENCE  EFFICIENCIES  FOR 
TWO  TYPES  OF  OPTIMUM  NOZZLE  CONTOURS 


Operating 

Mode: 

Mode  I 
£=  45:1 

Mode  II 
£=  222:1 

Secondary  Nozzle 

MODE  I 

MODE  II 

MODE  I 

MODE  II 

Optimizied  For: 

(Contour  2) 

(Contour  3) 

(Contour  2) 

(Contour  3) 

*^D  IV 

.991 

.988 

.982 

.985 

Source:  TDE  Option  TDK  Program 

The  absolute  values  of  the  divergence  efficiencies  shown  in 
Table  III-IX  are  somewhat  questionable  because  of  the  simplified  treatment  of 
the  flow  recompression  used  in  the  TDE  analysis.  The  relative  efficiency 
values  appear  to  be  valid,  however.  Therefore,  optimizing  the  secondary  noz¬ 
zle  contour  for  Mode  I  or  Mode  II  will  only  slightly  change  the  performance. 

F.  DUAL  THROAT  PERFORMANCE 

Many  of  the  model  improvements  described  previously  have  a  direct 
impact  on  the  performance  prediction  methodology  for  dual  throat  engines. 
Improvements  in  the  aerodynamic  model  have  resulted  in  less  restrictive  calcu- 
lational  assumptions,  tying  it  directly  to  the  shear  model  and  overall  bound¬ 
ary  layer  model  via  the  evaluation  of  the  momentum  thickness  at  the  plume 
attachment  point.  Results  of  detailed  two-dimensional  flow  analyses  (VNAP  and 
TDE)  have  further  verified  the  fluid  conditions  within  the  dual  throat 
thruster  and  the  resulting  performance  efficiencies. 


Ill,  F,  Dual  Throat  Performance  (cont.) 


Presented  within  this  section  are  the  results  of  an  assessment  in 
the  changes  in  the  predicted  dual  throat  performance  based  on  the  use  of  the 
improved  aerodynamic  model  and  updated  simplified  performance  model.  First, 
an  analysis  of  the  baseline  design  condition  for  Mode  II  is  presented  and  com¬ 
pared  to  a  conventional  gas-generator  cycle  engine.  Secondly,  the  performance 
of  a  limited  set  of  parametric  cases  calculated  during  NAS  8-32967  is  compared 
to  values  calculated  by  using  the  updated  simplified  performance  model. 

1.  Baseline  Design  Performance 

The  baseline  operating  point,  shown  in  Table  III-X,  is  the 
same  as  that  recommended  in  Contract  NAS  8-32967.  For  this  design  condition, 
a  bleed  flowrate  of  5.3%  is  predicted  to  be  required  for  shock-free  (zero 
pressure  loss)  plume  attachment.  However,  the  bleed  flow  can  be  reduced  to 
3.5%  with  only  a  small  loss  in  total  pressure  (vl%),  and  thus  this  test  condi¬ 
tion  was  selected  as  a  reasonable  design  tradeoff  for  this  comparison. 

The  calculated  specific  impulse  for  Mode  II  operation  is 
shown  in  Table  III -X I  and  compared  to  a  conventional  gas-generator  cycle 
engine  operating  at  the  same  condition.  The  results  of  this  comparison  show 
that  the  performance  of  the  dual  throat  engine  in  Mode  II  is  approximately  6 
seconds  (1-1/4%)  lower  than  that  of  the  conventional  engine.  As  noted,  most 
of  this  difference  in  performance  is  attributed  to  the  less  efficient  expan¬ 
sion  process  of  the  dual  throat  engine  in  Mode  II  which  results  in  a  greater 
divergence  efficiency.  This  performance  penalty  must,  of  course,  be  traded 
off  with  other  dual  throat  advantages  as  a  result  of  its  unique  capabilities. 

2.  Comparison  of  Simplified  Performance  Parametrics 

The  effect  of  the  previously  described  model  improvements  on 
performance  parametric  trends  was  briefly  investigated  by  using  the  updated 


m 


92 


TABLE  III-X 


DUAL  THROAT  BASELINE  OPERATING  POINT 


MODE  I  SEA  LEVEL  THRUST  2.669  MN  (600K  IbF) 
MODE  I  TO  MODE  II  THRUST  RATIO  3.2 


PRIMARY  NOZZLE 

CHAMBER  PRESSURE 
MIXTURE  RATIO 
PROPELLANT 

SECONDARY  NOZZLE 

CHAMBER  PRESSURE 
MIXTURE  RATIO 
PROPELLANT 


27.58  MPa  (4000  psia) 
7.0 

lox/lh2 


19.31  MPa  (2800  psia) 
2.8 

LOX/RP-1 
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TABLE  III-XI 


COMPARISON  OF  THE  DUAL  THROAT  AND  CONVENTIONAL 
GAS-GENERATOR  ENGINE  CONCEPTS 


PROPELLANTS: 

lo2/lh2 

MIXTURE  RATIO: 

6.0  (Overall) 

CHAMBER  PRESSURE: 

27.58  MPa  (4000  psia) 

AREA  RATIO: 

222:1 

THRUST  (VAC): 

956  KN  ( 21 5K  lbF) 

TOTAL  FLOWRATE: 

215.6  Kg/sec  (475  lbm/sec) 

DUAL  THROAT 

-  MODE  II 

CONVENTIONAL 

GG  CYCLE 

BLEED  FLOW 

=  3.5% 

GG  FLOW  = 

3.5% 

PRIMARY  O/F 

=  7.0 

TCA  O/F  = 

7.0 

ISP  ODE 

=  480.4  sec. 

lSP  ODE  = 

48D.4  sec 

nDIV 

=  .988 

nDIV 

.996 

nERE 

=  .990 

nERE 

.990 

nKIN 

=  .999 

nK  IN 

.999 

AISP  BL 

=  13.4  sec. 

A  *SP  BL 

5.7  sec. 

*SP  TCA 

=  456.2  sec. 

ISP  TCA  = 

467.5  sec 

!SP  ENG 

=  456.2  sec. 

!SP  ENG 

462.2  sec 
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III,  F,  Dual  Throat  Performance  (cont.) 

dual  throat  simplified  model.  The  results  for  five  parametric  cases  examined 
during  NAS  8-3/196 /  were  compared  with  the  present  model  predictions  for  both 
Mode  I  and  II  operation.  These  cases  were  used  to  investigate  the  effects  of 
chamber  pressure  and  thrust  ratio  (FVdCj  Mode  i/Fvac  Mode  2 )  on  en9lne 
performance. 


In  Mode  I,  the  vacuum  specific  impulse  predictions  and  trends 
are  essentially  unchanged  from  the  previously  reported  results.  However, 
because  the  calculational  technique  was  changed  to  require  equal  streamtube 
static  pressure  at  the  secondary  throat  during  Mode  1  operation,  a  significant 
difference  in  the  sea  level  specific  impulse  trend  with  thrust  ratio  was 
noted.  This  is  illustrated  in  Figure  111-34.  The  current  model  predicts  less 
variation  in  the  sea  level  specific  impulse  at  the  extremes  of  the  thrust 
ratio.  Note  that  at  the  baseline  thrust  ratio  of  3.2,  the  present  model  pre¬ 
dicts  a  sea  level  specific  impulse  which  is  approximately  6  sec  greater  then 
the  previous  model  predictions.  Engine  geometry  trends  with  thrust  ratio  were 
also  impacted,  as  noted  in  Fiyure  111-35.  The  secondary  throat  size,  rather 
than  increase  as  previously  predicted,  is  now  expected  to  actually  decrease 
slightly  with  thrust  ratio  (while  holding  sea  level  thrust  constant). 

In  all  cases,  t  fie  Mode  1!  specific  impulse  is  lower  with  use 
of  the  updated  performance  model  because  of  larger  boundary  layer  and  diver¬ 
gence  losses.  The  magnitude  of  this  difference  varies  from  case  to  case,  but 
is  generally  1  to  2%  less.  While  the  magnitude  of  the  Mode  II  performance  is 
less  than  previously  predicted,  the  certainty  of  the  prediction  is  much 
greater  as  a  result  of  the  work  completed  during  this  study.  In  addition,  as 
noted,  the  Mode  II  performance  predicted  for  the  recommended  baseline  is  only 
1-1/4%  less  than  that  of  a  conventional  gas-generator  engine,  and  thus  the 
previous  simplified  model  predictions  appear  optimistic  ii  this  light. 
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Ill,  Dual  Throat  Aerodynamic  Model  and  Performance  Prediction 
Improvement  (cont.) 

G.  RECOMMENDATIONS 


Design 


Droach 


The  previous  methodology  determines  primary  nozzle  area  ratio 
by  the  following  correlation: 


where  EPPM  is  an  empirical  factor.  A  value  of  EPPM  =  1.8  was  chosen  on  the 
basis  of  the  cold-flow  test  results  from  NAS  8-32666.  Nozzle  spacing  was 
chosen  to  be  0.8,  the  length  of  the  conical  portion  of  the  secondary  conver¬ 
gence  section. 


The  primary  nozzle  area  ratio  and  axial  spacing  between  the 
primary  and  secondary  throats  should  be  chosen  to  ensure  that  the  secondary 
flow  is  choked  at  the  secondary  throat,  rather  than  at  the  primary  lip,  by 
allowing  enough  flow  area  on  the  secondary  side  of  the  primary  nozzle  lip. 

This  design  approach  requires  small  primary  nozzle  area  ratios  and  large  noz¬ 
zle  spacing.  Secondly,  the  nozzle  spacing  chosen  must  result  in  the  ambient 
pressure  at  the  primary  nozzle  lip  being  low  enough  to  cause  the  primary  flow 
to  be  choked.  This  implies  that  the  nozzle  spacing  should  be  small.  Analyses 
of  bleed  flowrate  trends  for  Mode  II  indicate  that  bleed  flow  requirements 
decrease  with  increasing  primary  nozzle  area  ratio  and  decreasing  nozzle 
spacing.  These  trends  are  shown  in  Figure  1 1 1-36-  The  primary  nozzle  area 
ratio  and  nozzle  spacing  chosen  must  be  a  compromise  of  all  these  require¬ 
ments. 
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2.  Test  Programs 

The  aerodynamic  bleed  flow  model  has  three  empirical  param¬ 
eters  (identified  in  Section  II1,B).  Two  of  these  parameters  have  been  cali¬ 
brated:  the  plume  scaling  factor  and  the  NASH  factor,  using  cold-flow  and 
NASA/MSFC  test  data.  The  third  parameter,  the  spreading  factor,  still  needs 
to  be  calibrated.  This  parameter  will  influence  the  value  of  the  mass  frac¬ 
tion  of  primary  flow  entrained  in  the  recirculation  region  and  the  resulting 
recirculation  temperature  calculated.  This  may  impact  the  cooling  require¬ 
ments  in  the  plume  attachment  region.  A  test  design  is  required  that  allows 
thermal  data  acquisition  in  these  regions.  This  may  be  done  either  through  a 
calorimetric  convergence  section  or  placement  of  thermocouples. 

3.  Model  Improvements 

The  aerodynamic  bleed  flow  model  computer  program  currently 
determines  the  plume  boundary  shape  by  using  a  correlation  of  the  form  sug¬ 
gested  by  Herron.  The  intersection  of  the  plume  boundary  and  the  secondary 
nozzle  wall  is  determined  and,  together  with  the  slope  of  these  two  surfaces, 
used  to  calculate  recompression  criteria  for  the  shear  layer.  Because  a 
correlation  is  used,  the  plume  boundary  is  limited  in  range  and  is  necessarily 
somewhat  inaccurate.  The  effect  of  specific  heat  ratio  is  grossly  simplified. 
The  point  of  intersection  and  the  slope  of  the  plume  and  wall  at  the  inter¬ 
section  are  also  in  error.  The  error  in  streamline  slope  change,  i.e.,  the 
deflection  angle  of  the  plume  streamline,  can  be  more  than  20%. 

Figure  III-37a  shows  isobars  calculated  by  the  MOC  for  a 
cold-flow  nozzle  configuration.  The  streamline  shape  was  calculated  by  the 
correlation  discussed  above.  The  plume  streamline  should  be  constant  pres¬ 
sure,  but  is  actually  correct  to  within  only  about  15%.  (The  isobar  labeled 
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Ill,  G,  Recommendations  (cont.) 

X  approximately  follows  the  plume  streamline.)  Figure  III -37b  shows  charac¬ 
teristic  lines  calculated  by  the  MOC  interior  to  the  nozzle  and  plume  boundary 
streamlines.  In  Figure  1 1 1 -37c ,  the  flow  field  has  been  divided  into  two 
regions  that  are  separated  by  the  RRC  that  comes  down  from  the  primary  nozzle 
lip.  Region  I,  upstream  of  this  RRC,  is  independent  of  the  plume  shape. 

Region  I  can  be  calculated  by  using  the  TDK  Program  and  then  input  to  the 
Bleed  Flow  program. 

Next,  Region  II  can  be  calculated  and  used  to  replace  the 
plume  correlation.  It  is  defined  as  follows: 

Upstream  boundary:  RCC  ending  Region  I 

Left-hand  corner:  Prandtl -Meyer  fan  expanded  to  the 

base  pressure 

Upper  boundary:  Constant  pressure  boundary  (=  base 

pressure) 

Downstream  boundary:  LRC  from  the  upstream  boundary  to 

the  plume  attachment  point  on  the 
secondary  nozzle 

It  would  be  necessary  to  calculate  Region  II  using  the  MOD 
within  the  inner  search  procedure  of  the  Bleed  Flow  program,  i.e.,  directly 
replace  the  plume  correlation  with  a  MOC  calculation.  As  can  be  seen  from 
Figure  1 1 1-37b ,  however,  it  is  not  necessary  to  calculate  many  points  (only 
perhaps  300).  The  advantage  would  be  that  the  MOC  solution  will  be  highly 
accurate.  Thus,  it  is  recommended  that  the  aerodynamic  model  be  modified  to 
incorporate  a  MOC  calcualtion  for  the  plume  boundary. 
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IV. 


DUAL  EXPANDER  ANALYSIS 


A.  PRELIMINARY  GEOMETRIC  ANALYSIS 

This  study  was  the  first  effort  in  technically  analyzing  the  dual 
expander  nozzle  concept.  The  analysis  determined  the  dependence  of  perform¬ 
ance,  throat  area,  nozzle  and  steamtube  area  ratios,  mass  flow,  and  nozzle 
length  with  the  variation  in  Mode  I  sea  level  thrust,  thrust  ratio  (Mode  1/ 
Mode  II),  primary  and  secondary  chamber  pressures,  chamber  pressure  ratio 
(Pcs/Pcp),  and  nozzle  exit  pressure.  The  parametric  range  analyzed  was  as 
follows: 

Thrust  2.22  -  8.90  MN  (5  x  105  -  2  x  10&  lbF) 

Thrust  Ratio  1.0  -  5.0 

Chamber  6.89  -  55.16  MPa  (1000  -  8000  psia) 

Pressure 

Pc  Ratio  0.25  -  1.0 

Exit  Pressure  34.5  -  101  KPa  (5.0  -  14.7  psia) 

To  ensure  a  valid  comparison  with  the  dual  throat  analysis,  L0X/ 
RP-1  was  chosen  as  the  primary  nozzle  propellant,  and  LOX/LH^  was  chosen  for 
the  secondary  nozzle.  Mixture  ratios  used  were  2.8  and  7.0  for  the  primary 
and  secondary  chambers,  respectively. 

Table  IV-I  lists  the  parametric  values  used  for  each  case.  The 
resulting  performance,  area  ratios,  throat  areas,  flowrates,  and  efficiencies 
are  listed  in  Table  IV-I  I  - 

Mode  I  performance  calculations  were  based  on  the  simplified 
JANNAF  methodology  and  were  similar  to  the  streamtube  approach  used  during  the 
Dual  Throat  Geometry  Program  (FD  0169).  Mode  II  performance  was  determined  by 
applying  the  simplified  JANNAF  methodology  to  a  Bell  nozzle  having  the  same 
throat  area,  area  ratio,  and  length  as  the  secondary  nozzle  during  Mode  II 
operation.  The  procedures  were  programmed  into  a  Dual  Expander  Preliminary 
Geometry  and  Performance  Prediction  Program. 


TABLE  IV- I 


DUAL  EXPANDER  CASES  ANALYZED 


CASE 

SEA  LEVEL 

THRUST 

MN  (lbF) 

THRUST 

RATIO 

p 

^CP 

MPa  (psia) 

pcs 

MPa  (psia) 

p 

tXIT 

KPa  (psia) 

PCS^ 

1 

2.67  (6 

X 

105) 

2.37 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

2 

2.67  (6 

X 

105) 

2.37 

41.37  (6000) 

10.34  (1500) 

34.5  (5.0) 

.25 

3 

2.67  (6 

X 

105) 

2.37 

41.37  (6000) 

31.03  (4500) 

34.5  (5.0) 

.75 

4 

2.67  (6 

X 

105) 

2.37 

41.37  (6000) 

41.37  (6000) 

34.5  (5.0) 

1.0 

5 

2.67  (6 

X 

105) 

1 .25 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

6 

2.67  (6 

X 

105) 

3.50 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

7 

2.67  (6 

X 

io5) 

4.50 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

8 

2.67  (6 

X 

105) 

5.0 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

9 

2.67  (6 

X 

!06) 

2.37 

13.79  (2000) 

6.89  (1000) 

34.5  (5.0) 

.5 

10 

2.67  (6 

X 

105) 

2.37 

27.58  (4000) 

13.79  (2000) 

34.5  (5.0) 

.5 

11 

2.67  (6 

X 

io5) 

2.37 

55.16  (8000) 

27.58  (4000) 

34.5  (5.0) 

.5 

12 

2.22  (5 

X 

io5) 

2.37 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

13 

4.44  (1 

X 

106) 

2.37 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

14 

6.67  (1. 

.5 

x  106) 

2.37 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

15 

8.90  (2 

.0 

6x 

x  10  ) 

2.37 

41.37  (6000) 

20.68  (3000) 

34.5  (5.0) 

.5 

16 

2.67  (6 

X 

IO5) 

2.37 

41.37  (6000) 

20.68  (3000) 

68.9  (10.0) 

.5 

17 

2.67  (6 

X 

io5) 

2.39 

41.37  (6000) 

20.68  (3000) 

101  (14.7) 

.5 

18 

2.67  (6 

X 

105) 

1.0 

- 

20.68  (3000) 

34.5  (5.0) 

- 
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IV,  A,  Preliminary  Geometric  Analysis  (cont.) 


Figures  IV-1  through  IV-3  show  the  dependence  of  engine  perform¬ 
ance  and  area  ratio  on  thrust  split.  Decreasing  Mode  I  performance  as  thrust 
ratio  increases  is  due  to  the  increased  contribution  of  the  lower- performing 
primary  flow  to  the  overall  performance.  The  Mode  II  delivered  Isp  increases 
with  thrust  ratio  because  the  Mode  II  area  ratio  is  increased.  This  increase 
in  area  ratio  occurs  because  the  secondary  throat  decreases  with  respect  to 
the  nozzle  exit  area  as  the  thrust  ratio  increases. 


The  Mode  I  area  ratio  increases  with  thrust  split  because  of 
increases  in  primary  throat  area  with  thrust  split.  From  geometry,  the  Mode  I 
area  ratio  can  be  defined  as: 


£ 


I 


A.  /A+_  +  e 
p  tp  ts 

+ir 


s 


The  primary  (£p)  and  secondary  (€$)  streamtube  area  ratios  are  constant 
for  a  given  ratio  of  specific  heats  (r)  and  pressure  ratio  (Pe/Pc) . 

Therefore,  increasing  Atp  results  in  a  larger  £i* 

Figures  IV-4  and  IV-5  show  the  effect  of  variation  in  thrust  on 
specific  impulse  and  throat  area.  Throat  area  will  increase  with  increasing 
thrust,  but  the  specific  impulse  will  essentially  be  unchanged. 

Figures  IV-6  through  IV-8  show  the  effects  of  chamber  pressure  on 
throat  area,  area  ratio,  and  specific  impulse.  Increasing  chamber  pressure  or 
chamber  pressure  ratio  results  in  decreasing  throat  areas  and  increasing  area 
ratio,  in  turn  resulting  in  an  increase  in  specific  impulse. 

Figures  IV-9  and  IV-10  show  how  decreasing  nozzle  exit  pressure 
will  result  in  increasing  nozzle  area  ratio  and  specific  impulse. 
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Figure  IV-2.  Engine  I  versus  Thrust  Ratio 


THRUST-  ,  =  6  x  10  LBF 


THRUST  RATIO 


6000  PSIA 


Throat  Area  versus  Thrust 


versus  Thrust 


Figure  iv-7.  Throat  Area  versus  Chamber  Pressure 


versus  Chamber  Pressure 


6000  PSIA 


Nozzle  Area  Ratio  versus  Exit  Pressure 
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IV,  A,  Preliminary  Geometric  Analysis  (cont.) 

The  previously  listed  geometric  parameters  were  used  in  deter¬ 
mining  the  dual-expander  nozzle  design.  The  primary  nozzle  contour  was  found 
by  truncating  a  bell  nozzle  optimized  for  the  primary  flow  field  area  ratio  to 
some  specified  length  or  area  ratio.  An  annular  nozzle  with  a  conical  inter¬ 
nal  expansion  section  was  positioned  such  that  the  flow  angle  and  static  pres¬ 
sure  at  the  exit  of  the  internal  expansion  section  matched  the  flow  angle  and 
static  pressure  at  the  primary  nozzle  exit. 

These  steps  are  part  of  the  procedure  developed  for  designing  the 
dual-expander  nozzle.  The  nomenclature  for  this  procedure  is  shown  in  Figure 
IV— 1 1.  Derivation  of  the  equations  is  given  in  Appendix  B. 

The  design  procedure  is  as  follows: 

1.  Calculate  the  overall  engine  area  ratio,  primary  and  secon¬ 
dary  flows  streamtube  area  ratios,  throat  areas,  and  equivalent  nozzle  lengths 
for  a  given  set  of  operating  conditions. 

2.  Generate  a  Rao  nozzle  contour  for  the  primary  nozzle  by  using 
the  calculated  streamtube  area  ratio  and  nozzle  length.  Truncate  this  contour 
to  the  desired  area  ratio  value.  The  primary  nozzle  exit  radius  (REp)» 

flow  angle  (£>l)»  anc*  2-0  pressure  are  now  known. 

3.  Once  the  secondary  nozzle  chamber  pressure  and  the  desired 
exit  pressure  are  known,  calculate  the  required  conical  section  area  ratio.  A 
conical  section  is  used  because  the  2-D  to  1-D  Mach  number  ratio  will  be  close 
to  one.  The  conical  section  exit  area  (A£c )  can  now  be  calculated  for  the 
secondary  flow. 


i 
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IV,  A,  Preliminary  Geometric  Analysis  (cont.) 

4.  Pick  a  value  for  either  the  gap  {&£c)  or  the  radius 
(R£c)  at  the  conical  section  exit  and  solve  for  t he  other  value  by  using 
Equation  1. 


2tt  6^ 

c  c 


O) 


5.  Choose  an  inner  lip  thickness  (Dl)  and  calculate  the  secon¬ 
dary  flow  inclination  angle  (0S)  from  Equation  2. 


Cos  8 


2  ( Rr  -  Rf  -  D  Cos  6  ) 

_ Lc  bp _ _ _ _ 


(2) 


aE. 


6.  The  exit  geometry  for  the  conical  section  is  now  known.  The 
conical  half-angle  can  also  be  found  by 


9c  =  9s  '  0L 


The  next  step  is  to  find  the  proper  secondary  nozzle  throat 
radius  (Rts)  and  gap  (<&ts)  that  occur  for  the  inclination  angle,  0S, 
and  secondary  nozzle  throat  area  (At s )  calculated  in  Step  1. 

7.  The  secondary  nozzle  throat  radius  (Rts)  is  calculated 
from  Equation  4. 


Rt 


2  sec  8, 


(’ 


sin  9L 


sin  0$  cos  8 


4 


Rt. 


s  =  0 


2  sec  6.  I  RE_  +  D,  cos  8,  -  REc  sin  0L 
5  1  p  L  L  sin  0s  cos  8c 


(4) 


TT 


Equation  (4)  is  of  the  form 


ax  +  bx  +  c  =  0 
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Rt  =  -b  +  |  b  -  4ac 


2a 


/  sin  0. 

a  =  2  sec  e$  M  sTn  F'  coTT 

\  S  c 


b  =  -2  sec  0  /  Rr  r  D.  Cos  0. 

s  1  Ep  L  L 


sin  0p 

_ _ 

sin  0  cos  0 
s  c 


2- 

8.  Once  Rt  is  known,  the  throat  gap  (<s  )  can  be  found 


by  using  Equat ion  5: 


2'S  t  Rt 
s  Ls 


(5) 


The  length  of  the  conical  section  (Xc)  is  found  by  the  following  equation: 


R 


F 


R 


t 


c  s 

Sin  9 

s 


(6) 


9.  The  remainder  of  the  secondary  nozzle  contour  is  approximated 
with  a  third-order  parabola: 


R 


A  +  BX  + 


2 

CX  + 


DX~ 


The  constants  A,  B,  C,  D  are  determined  by  specifying  the  coordinates  and 
angle  of  the  outer  conical  section  wall  and  specifying  the  coordinates  and 
angle  at  the  overall  nozzle  exit.  The  exit  radius  is  found  from  Step  1,  and 


IV,  A,  Preliminary  Geometric  Analysis  (cont.) 


the  exit  angle  is  assumed  to  be  the  same  as  calculated  for  the  complete  pri¬ 
mary  nozzle  contour. 

Equations  (1)  through  (6)  show  that  once  the  exit  radius 

(Rff.)  and  gap  )  are  picked,  there  is  only  one  set  of  values, 
c 

es,  Rts>  6ts.  xc,  that  will  satisfy  these  equations  and  result  in 
a  symmetrical  conical  nozzle  section.  Table  I V- III  lists  different  values  of 
rEc»  5Ec»  and  °L  an<*  their  influence  on  throat  gap  and  location, 
throat  orientation,  and  length  of  the  conical  section.  For  a  constant  lip 
thickness,  increasing  the  exit  radius  results  in  decreasing  the  exit  gap,  the 
throat  orientation  angle,  and  the  distance  of  the  secondary  throat  from  the 
central  axis.  The  secondary  throat  gap  and  the  length  of  the  conical  section 
both  increase  with  increasing  RE(;.  Increasing  the  lip  thickness  for  a 
given  REc  and  <SEc  causes  the  secondary  throat  inclination  angle  to 
increase  and  moves  the  throat  further  from  the  central  axis.  Both  the  throat 
gap  and  conical  section  length  get  smaller. 

B.  DUAL-EXPANDER  AERODYNAMIC  BLEED  FLOW  MODEL 

1.  Model  Development 

An  aerodynamic  model  has  been  developed  for  Mode  II  operation 
and  is  presented  in  this  section.  This  model  differs  from  the  dual-throat 
thruster  Mode  II  aerodynamic  model  in  several  respects.  First,  the  input 
geometry  for  the  dual  expander  aerodynamic  model  is  very  simple.  It  is  only 
necessary  to  specify  the  throat  area,  AE$,  and  lip  angle,  The  plume 
model  is  also  very  simple.  The  flow  at  the  exit  of  the  annular  nozzle  is 
assumed  to  be  one-dimensional  in  a  cross  section  normal  to  the  nozzle  lip. 

The  flow  is  assumed  to  follow  a  Prandtl -Meyer  expansion  around  the  nozzle  lip 
until  an  angle  is  reached  such  that  the  pressure  of  the  base  region  is 
achieved.  The  plume  boundary  is  taken  to  be  a  straight  line  inclined  at  this 
angle. 


121 


DUAL  EXPANDER  GEOMETRIC  INELUENCES 


IV,  B,  Dual -Expander  Aerodynamic  Bleed  Flow  Model  (cont.) 


The  shear  layer  is  superimposed  on  the  plume  boundary.  With 
one  exception,  the  assumptions  made  in  formulating  the  shear  layer  model  are 
the  same  as  those  used  for  the  dual  throat  (described  in  Section  III).  The 
exception  is  that  the  shear  layer  is  assumed  to  have  axial  symmetry.  Assuming 
axial  rather  than  planer  symmetry  significantly  complicates  the  formulation; 
however,  as  discussed  in  Section  III,  this  assumption  is  necessary  for  the 
dual  expander  because  of  the  proximity  of  the  shear  layer  to  the  flow  axis. 

In  Reference  18,  Mueller  has  presented  an  analysis  for  an  axisymmetric  shear 
layer  for  a  gas  with  constant  properties.  The  analysis  presented  here  allows 
for  the  addition  of  bleed  flow  into  the  base  recirculation  region.  As  was  the 
case  with  the  dual  throat  aerodynamic  bleed  flow  model,  the  bleed  flow  gas  can 
have  properties  different  from  those  of  the  primary  stream.  The  dual  expander 
shear  model  is  illustrated  in  Figures  I V- 1 2  and  I V- 1 3.  Figure  I V- 1 2  shows  the 
shear  layer  developing  along  the  jet  boundary  until  a  position  is  reached 
where  the  shear  layer  is  turned  and  transformed  into  a  wake.  The  position  at 
which  the  shear  layer  is  turned  is  called  the  "neck."  Turning  the  flow 
abruptly  at  the  neck  results  in  a  recompression  shock  (shown  in  Figure  I V- 1 2 ) . 
The  distance  from  the  exit  of  the  interior  nozzle  to  the  neck,  Z3»  is  shown 
in  Figure  I V- 1 3.  The  streamline  deflection  angle  at  the  neck,  83,  and  the 
trailing  wake  radius,  rw,  are  also  shown  in  Figure  IV-13.  The  cal culational 
procedure  is  illustrated  in  Figure  IV-14.  As  was  the  case  for  the  dual  throat 
aerodynamic  bleed  flow  model,  the  problem  requires  solving  three  nonlinear 
simulateous  equations  for  three  unknowns.  The  three  unknown  quantities  are 
Pb/^op*  Hr/Hop  and  Ypr.  The  three  equations  are  (A),  (B),  and  (C) 
in  Figure  IV-14.  These  equations  are: 


7T  Z  , 


VWp 


cos 
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igure  IV-13.  Dual  Expander  Aerodynamic  Model  Geometry 


INPUT: 


Gas  Propertiesj 
Cp,  Mw,  AND  H0 


FOR  BOTH  STREAMS 
FOR  OUTER  STREAM 


Mach  NO. 

Recompression  Factor,  N 
%  Bleed. 


Geometryj 


R  '  9B''  rw*  rbj  ap 


Estimates 


METHOD: 


PB/Po  ,  '  VHo  '  Ypr 


Newton's  Method,  given  above  estimates  find 
converged  values  that  satisfy  condition  of 

(A)  continuity  for  total  flow 

(B)  energy  conservation 

(C)  continuity  of  base  flow  compostion 


Special  Case:  0%  BLEED. 


VP0  '  VH0  '  Y 

CONSTANT  PRESSURE  BOUNDARY  SHAPE 


Figure  IV-14.  Calculational  Procedure  for  the  Dual  Expander  Aerodynamic  Model 


OUTPUT: 
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In  the  above  equations,  the  integrals  I i ,  I 3 ,  Ji,  and  J3  are 
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IV,  B,  Dual -Expander  Aerodynamic  Bleed  Flow  Model  (cont.) 

An  input  manual  for  the  dual  expander  aerodynamic  bleed  flow 
model  is  presented  in  Appendix  C. 

2.  Model  Calibration 


The  dual  expander  base  flow  model  contains  the  following 
three  parameters  that  must  be  determined  from  experimental  data: 

o  =  The  jet  spreading  parameter 

N  =  The  Nash  recompression  factor 

k  =  A  proportionality  factor  relating  the  radius  of  the 
trailing  wake  to  the  shear  layer  velocity  thickness 

The  approach  taken  in  selecting  each  of  these  parameters  is 
discussed  below.  Note  that  all  of  the  data  used  is  for  air  {y  =  1.4)  without 
bleed  flow. 


o,  the  Jet  Spreading  Parameter 

The  classical  jet  spreading  expression  measured  by  Korst 
(Ref.  6)  for  air  was  used,  i.e., 


o=12+  2.758Mb 

For  incompressible  flow  (M^  =  0),  this  expression  gives  the  well-known 
result  that 


a 


0 


12. 


Many  alternate  expressions  for  o  are  available  in  the  literature  (see  Ref.  7, 
11,  18). 


IV,  B,  Dual -Expander  Aerodynamic  Bleed  Flow  Model  (cont.) 

N,  the  Nash  decompression  Factor 

The  data  presented  by  Mueller  (Ref.  18)  was  used  to  determine 
N.  This  data  is  for  the  flow  of  air  past  a  cylinder.  Measured  values  of  base 
pressure  ratio  vs  the  Mach  number  of  the  approach  flow  are  presented  in  Figure 
IV-15.*  The  theory  predicts  that  base  pressure  depends  only  weakly  on  the 
trailing  wake  radius  ratio.  Measured  values  of  the  trailing  wake  radius  ratio 
are  shown  in  Figure  IV-16. 

The  computer  model  was  calibrated  by  inputting  the  flow 
conditions  that  correspond  to  the  experimental  conditions  and  by  inputting  the 
trailing  wake  radius  from  Figure  IV-16.  It  was  found  that 

N  =  .4 

gave  the  best  fit  to  the  data.  This  is  consistent  with  the  treatment  used  in 
the  dual  throat  aerodynamic  model  computer  program.  It  should  be  noted  that 
Mueller  (Ref.  18)  used  a  value  of  N  =  1  (i.e.,  the  full  value  from  the  shock) 
and  a  different  expression  for  a  o  (i.e.,  an  expression  that  gives  values  that 
are  greater  than  those  given  by  the  Korst  expression). 

K,  Shear  Layer/Wake  Proportionality  Factor 

For  the  case  of  air  without  bleed  flow,  the  trailing  wake 
radius  ratio  can  be  obtained  from  Figure  IV-16.  In  order  to  obtain  a  similar 
result  for  situations  with  bleed  flow  and  with  gases  other  than  air,  the 
following  assumption  was  made. 


^Figures  IV-15  and  IV-16  are  from  Reference  18. 
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IV,  B,  Dual-Expander  Aerodynamic  Bleed  Flow  Model  (cont.) 

It  was  assumed  that  the  radius  of  the  trailing  wake  is 
proportional  to  the  velocity  thickness  of  the  shear  layer  at  the  neck,  6. 
Since  this  value  is  proportional  to 

(x/r b)/o 

where  x/r^  is  the  nondimensional  shear  layer  length,  it  follows  that 

Vrb  =  k  (x/rb)o 

where 

K  is  the  proportional i ty  constant. 

It  was  found  that 


K  «  4 

i.e.,  this  approximate  value  was  consistently  calculated  when 

a  =  12  +  2.758Mb 

N  =  .4 

and  the  cases  shown  in  Figure  IV-16  were  calculated. 

C.  DUAL  EXPANDER  FLOW  FIELD  ANALYSIS 

Perfonnance  calculations  were  carried  out  for  both  Mode  I  and  Mode 
II  operation.  In  order  to  carry  out  these  calculations,  it  was  convenient  to 
divide  the  dual  expander  flow  field  into  regions  that  could  be  analyzed 
separately.  These  regions,  numbered  I  through  V,  are  shown  in  Figures  IV-17a 
and  IV-17b. 
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REGION  II 


REGION  III 


DUAL  EXPANDER  MODE  I  OPERATION 
BOTH  ENGINES  ON 


b.  DUAL  EXPANDER  MODE  II  OPERATION 
OUTER  ENGINE  ON,  INNER  ENGINE  OFF 


Figure  IV-17.  Dual  Expander  Modes  I  and  II  Operation 


IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 

Region  I  is  the  flow  field  inside  the  interior  thruster  during 
Mode  I  operation.  The  downstream  boundary  of  Region  I  is  the  exit  plane  of 
the  primary  thruster.  Region  I  can  be  conveniently,  accurately,  and  inexpen¬ 
sively  calculated  by  using  the  TDK  computer  program  (Ref.  3). 

Region  II  is  the  flow  field  inside  the  annular  thruster.  This 
flow  field  is  the  same  during  both  Mode  I  and  Mode  II  operation.  The  down¬ 
stream  boundary  of  Region  II  is  the  exit  plane  of  the  annular  thruster. 

Region  II  can  be  calculated  by  using  the  VNAP  computer  program  (Ref.  4)  that 
was  discussed  in  Section  III.  VNAP  is  particularly  well  suited  to  calculating 
the  transonic  flow  field  in  an  annular  nozzle.  The  VNAP  calculations  are 
somewhat  lengthy  and  require  dividing  Region  II  into  subregions  in  order  to 
obtain  sufficient  accuracy. 

Region  III  is  the  dual  expander  flow  field  during  Mode  I  operation 
(both  engines  on)  that  is  downstream  from  Regions  1  and  II.  In  order  to  mini¬ 
mize  shock  losses,  the  engine  is  designed  to  operate  so  that  the  static  pres¬ 
sures  in  Regions  I  and  II  are  equal  at  the  point  where  these  two  flow  streams 
come  together.  This  point  is  the  origin  of  a  contact  discontinuity,  or  slip 
line,  between  the  two  streams.  The  TDK  program  has  been  modified  to  calculate 
the  flow  field  in  Region  III.  In  order  to  use  the  option,  it  is  necessary  to 
input  flow  conditions  along  the  exit  planes  of  Regions  I  and  II.  The  two  flow 
streams  that  comprise  Region  III  are  allowed  to  have  different  values  of  stag¬ 
nation  pressure,  specific  heat  ratio,  and  molecular  weight. 

Region  IV  is  the  inviscid  portion  of  the  flow  field  that  occurs 
downstream  of  Region  II  during  Mode  II  operation  (inner  thruster  off,  outer 
thruster  on).  The  interior  boundary  of  Region  IV  is  determined  by  the  shape 
of  the  viscous-inviscid  boundary  of  the  base  recirculation  flow  and  wake. 
Region  IV  can  contain  shock  waves.  The  VNAP  computer  program  (Ref.  4)  has 
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IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 

been  used  to  calculate  the  flow  field  in  this  region  so  that  the  effect  of 
shock  waves  can  be  included.  In  order  to  maintain  accuracy,  it  is  necessary 
to  divide  Region  IV  into  subregions. 

Region  V  is  the  base  flow  reci rcul at  ion  and  wake  region  that  is 
calculated  by  the  dual  expander  aerodynamic  model  described  in  the  previous 
section. 


A  prototype  design  for  the  dual -expander  thruster  was  generated  by 
using  the  procedure  defined  in  Section  IV, A,  and  performance  calculations  were 
carried  out  by  analyzing  the  five  flow  regions  described  above.  These  calcu¬ 
lations  are  described  in  the  following  subsections. 

1.  Calculations  for  Mode  II  Operation 

Calculating  the  dual  expander  performance  for  Mode  II  opera¬ 
tion  requires  calculating  Regions  V,  II,  and  IV  that  were  shown  in  Figure 
I V - 1 7b.  Results  obtained  for  each  region  are  discussed  below.  The  purpose  of 
these  calculations  was  to  determine  the  influence  of  primary  nozzle  truncation 
and/or  the  introduction  of  bleed  flow.  In  order  to  isolate  these  influences, 
it  was  necessary  to  keep  the  same  secondary  nozzle  contour  and  vary  the  pri¬ 
mary  nozzle  exit  area  only. 

Region  V 

Region  V  was  calculated  for  the  prototype  engine  for  three 
bleed  flow  values  and  for  three  truncation  positions  for  the  dividing  wall. 
Results  from  these  nine  cases  are  presented  in  Table  I V- 1 V ,  labeled  as  Cases  2 
through  10.*  The  bleed  flow  values  used  were  zero  bleed,  "moderate"  bleed, 
and  "maximum"  bleed. 


*Case  1  was  used  as  a  sample  case  for  VNAP  for  program  checkout  prior  to  com¬ 
pleting  Cases  2  through  10. 


IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 


The  aerodynamic  model  is  a  small  bleed  flow  model.  If  the 
bleed  flowrate  becomes  too  large,  it  will  lead  to  the  following:  the  base 
pressure  region  will  not  close;  there  will  be  no  recirculation,  the  model  will 
no  longer  be  applicable.  "Maximum"  bleed  was  defined  as  an  amount  just  short 
of  causing  this  condition.  "Moderate"  bleed  flow  was  taken  as  half  of  the 
"maximum"  bleed  flow  value. 

"Minimum"  truncation  corresponds  to  arbitrarily  choosing  to 
truncate  the  primary  nozzle  to  a  geometric  area  ratio  of  8:1.  "Maximum"  wall 
truncation  was  obtained  by  truncating  the  dividing  wall  just  downstream  of  the 
Mach  one  position  on  the  annular  nozzle  side  of  the  wall.  The  "moderate"  wall 
truncation  position  was  taken  as  halfway  between  "maximum"  and  "minimum"  trun¬ 
cation. 


The  parameters  tabulated  in  Table  I V - 1 V  were  illustrated  in 
Figure  IV-13.  Inputs  to  the  dual  expander  aerodynamic  model  that  differ  in 
each  case  are  1)  the  exit  radius  of  the  interior  nozzle,  r^,  2)  the  Mach 
number  of  the  approach  flow  at  the  end  and  upper  side  of  the  dividing  wall, 
Mp,  and  3)  %  bleed  for  inner  engine.  The  approach  flow  angle  is  31°.  The 
operating  conditions  for  these  cases  are  listed  in  Table  IV-V. 

TABLE  IV-V 


MODE  II  OPERATING  CONDITIONS 

Outer  Nozzle 

Bleed  to 

Inner  Nozzle 

Flow 

LOX/LH2 

GG 

Mixture  Ratio 

7.0 

.72 

Y 

1.18 

1.389 

Cp 

.81371  cal / g m  °K 

2.04454  cal/gm 

Molecular  Weight 

16. 

3.467 

Total  Temperature 

3778°K  (6800°R) 

712°K  (1281°R) 

Total  Pressure 

20.63  MPa 

N/A 

Throat  Area 

(3000  psia) 

231  cm2  (43.5  in. 2) 

N/A 

IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 


The  geometry  of  the  shear  layer  and  wake  as  calculated  by  the 
dual  expander  aerodynamic  model  is  defined  by  the  parameters  0^ 
z3/rb>  ancl  rw/rb*  Calculated  values  for  these  parameters  are  listed 
in  Table  IV- 1 V.  Also  listed  is  the  value  of  the  shock  angle,  qs,  corres¬ 
ponding  to  a  deflection  angle  of  @^. 

The  approach  Mach  number  at  the  three  truncation  positions 
was  input  from  the  VNAP  results  for  Region  II.  The  shear  layer  and  wake  shape 
calculated  for  Region  V  was  used  to  define  the  lower  boundary  for  Region  IV. 
The  Region  II  and  IV  calculations  are  discussed  below. 

Region  II 

Figure  IV-18  shows  an  enlargement  of  Regions  I  and  II.  The 
upper  wall  of  Region  II  is  steep,  having  a  maximum  wall  angle  of  67°.  To 
obtain  accuracy,  Region  II  was  divided  into  an  upstream  subregion  and  a  down¬ 
stream  subregion,  and  these  were  calculated  separately  by  using  VNAP. 

Mach  number  contours  calculated  for  the  upstream  subregion 
are  shown  in  Figure  IV-19.  Inlet  conditions  for  this  subregion  were  pre¬ 
scribed  to  be  one-dimensional  parallel  flow  along  a  vertical  cross  section  at 
the  extreme  left  of  Figure  IV-19.  The  calculations  show  that,  along  the  lower 
wall,  the  flow  is  turned  inward  abruptly  while  subsonic  (M  =  .4).  Along  the 
upper  wall,  the  nozzle  throat,  a  circular  arc  with  a  small  radius  of  curva¬ 
ture,  causes  the  flow  to  be  choked.  The  mass  flowrate  through  the  nozzle  is 
found  to  be  only  about  91%  of  the  one-dimensional  flow  value.  At  the  exit  of 
this  subregion,  the  flow  is  found  to  be  supersonic,  and,  therefore,  indepen¬ 
dent  of  the  conditions  downstream. 
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IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 


Mach  number  contours  calculated  for  the  downstream  subregion 
of  Region  II  are  shown  in  Figure  IV-20.  Inlet  conditions  for  this  subregion 
are  the  exit  conditions  of  the  upstream  subregion.*  An  interesting  feature  of 
this  plot  is  that  the  flow  near  the  lower  wall  is  found  to  expand  like  a  one¬ 
dimensional  flow.  This  tends  to  validate  the  design  approach  of  using  a  coni¬ 
cal  section  in  the  annular  portion  of  the  secondary  nozzle  so  that  the  pres¬ 
sure  may  be  determined  by  using  one-dimensional  calculations.  As  was  dis¬ 
cussed  previously,  it  is  important  to  be  able  to  estimate  this  pressure  value 
when  designing  the  nozzle.  It  can  also  be  seen  from  Figure  IV-20  that  the 
flow  is  being  compressed  by  the  upper  wall.  Downstream,  this  compression 
develops  into  a  shock. 

The  Region  II  calculations  are  used  to  supply  inlet  condi¬ 
tions  to  Region  IV  and  also  to  the  Region  III,  Mode  I,  calculations. 

Region  IV 

Inlet  line  conditions  for  Region  IV  were  supplied  by  selecing 
the  appropriate  flow  cross  section  from  the  VNAP  results  for  Region  II.  The 
upper  wall  boundary  was  determined  by  using  the  design  approach  described  in 
Section  IV, A.  The  lower  boundary  was  defined  by  the  parameters  63) 

Z3/ rb,  and  rw/rt>- 


Region  IV  was  divided  into  an  upstream  subregion  and  a  down¬ 
stream  subregion.  The  following  5  plots  were  prepared  for  each  case: 


*Each  region  is  actually  overlapped  to  increase  the  accuracy.  The  last  line 
or  two  calculated  by  VNAP  is  not  accurate  because  of  the  outflow  boundary 
condition  that  must  be  imposed  by  the  method  to  close  the  problem.  (See 
Ref.  4.)  Thus,  the  VNAP  calculations  have  been  overlapped  and  the  last 
two  lines  were  discarded. 


IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 


(1)  velocity  field,  (2)  pressure,  (3)  density,  (4)  temperature,  and  (5)  Mach 
number  contours.  Only  the  Mach  number  contours  are  discussed  in  this  report. 
Mach  number  contours  for  these  subregions  are  presented  in  Figures  IV-21  and 
IV-22  for  case  2,  i.e.,  zero  bleed  and  no  truncation.  Thus,  Figures  IV-19, 
-20,  -21,  and  -22  show  the  entire  Mode  II  flow  field  predicted  for  the  base¬ 
line  nozzle  contour  and  zero  bleed  flow. 

The  lower  boundary  in  Figure  IV-21  is  the  shear  layer.  The 
lower  left-hana  corner  in  the  figure  shows  the  Prandtl-Meyer  expansion  fan  at 
the  lower  lip  of  the  annular  nozzle  where  the  flow  turns  through  45.64°.  The 
finite-difference  method  used  by  VNAP  cannot  accurately  calculate  this  expan¬ 
sion  fan  because  the  Mach  number  gradients  are  too  great  for  the  relatively 
coarse  mesh  spacing  required  by  VNAP.  Thus,  this  expansion  fan  is  a  source  of 
error  in  the  VNAP  calculations. 

The  lower  boundary  of  Figure  IV-22  is  the  shear  layer  and 
wake.  A  compression  shock  is  induced  at  the  corner  where  the  flow  is 
deflected  through  an  angle  of  14.64°.  This  deflection  angle  corresponds  to  a 
shock  angle  of  23.4°.  Figure  IV-22  clearly  shows  this  "neck"  shock;  however, 
the  shock  angle  appears  to  be  less  than  23.4°.  VNAP  predicts  that  the  neck 
shock  will  exit  the  nozzle  far  from  the  upper  wall. 

A  shock  can  also  be  seen  in  the  upper  half  of  Figure  IV-22. 
This  shock  is  induced  by  compression  of  the  flow  by  the  upper  wall. 

Resul ts 


Results  of  the  Mode  II  performance  predictions  are  summarized 

in  Table  IV-VI . 
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Downstream  Mach  No.  Contours,  Case  2,  Zero  Truncation,  0%  Bleed 


IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 

For  cases  2,  3,  and  4,  where  the  wall  dividing  the  upper  and 
lower  nozzles  were  as  shown  in  Figure  iV-18,  the  approach  Mach  number  was 
found  to  be  2.79.  With  zero  bleed  flow,  the  predicted  Prandtl-Meyer  turning 
angle  at  the  nozzle  lip  was  31  +  14.64  =  46.64°.  This  turn  gives  a  Mach  num¬ 
ber  of  4.62  for  the  shear  layer  boundary.  For  a  Mach  number  this  high,  very 
little  bleed  flow  is  required  to  open  the  wake.  The  aerodynamic  model  pre¬ 
dicted  a  trailing  wake  radius  ratio  approaching  one  for  a  bleed  flow  of  only 
.07%.  The  values  of  base  pressure  predicted  for  cases  2,  3,  and  4  are  shown 
in  column  one  of  Table  1 V- V I ,  and  are  found  to  be  quite  low.  For  cases  5,  6, 
and  7,  which  have  moderate  wall  truncation,  and  cases  8,  9,  and  10,  which  have 
maximum  wall  truncation,  the  approach  Mach  numbers  are  2.19  and  1.016,  respec¬ 
tively.  Consequently,  the  predicted  values  of  base  pressure  and  bleed  flow 
required  to  open  the  wake  are  larger  than  for  the  zero  truncation  cases. 

The  remainder  of  Table  1V-VI  summarizes  the  value  of  mass 
flow,  thrust,  and  specific  impulse  that  were  predicted  by  VNAP.  Mass  flowrate 
and  thrust  were  calculated  by  using  two  different  methods: 

1 )  Initial  Line  Method 

"'initial  =  Integral  across  the  VNAP  start  line  of 
line  Figure  IV-21. 

Fi  =  Integral  of  axial  moment  'm  +  pressure 

across  the  VNAP  stjrt  line  of  Figure 
IV-21,  +  integral  of  axial  pressure  along 
both  the  upper  and  lower  boundaries 
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IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 


2)  Exit  Line  Method 

Integral  across  the  VNAP  line  at  the 
end  of  the  nozzle 

Integral  of  axial  momentum  +  pressure 
across  the  VNAP  line  at  the  end  of  the 
nozzle 

Two  different  values  of  specific  impulse  were  calculated  by  using  the  above 
quantities: 

(ISP)l  _  Fl/^initial 
1  i  ne 

(ISP)2  =  F2/mexit 
1  i  ne 

Referring  to  Table  IV-VI,  it  can  be  seen  that  there  is  substantial  error 
growth  in  both  the  mass  flow  and  thrust  calculations,  especially  for  cases  8, 
9,  and  10.  This  is  primarily  because  of  the  inability  of  VNAP  to  properly 
calculate  the  large  gradient  region  in  the  vicinity  of  the  Prandtl-Meyer 
expansion  fan.  These  gradients  are  especially  large  in  the  last  three  cases 
because  of  the  low  Mach  number  (1.016)  of  the  flow  upstream  of  the  fan. 

Nearly  all  of  the  observed  error  growth  occurs  in  this  region,  i.e.,  in  cal¬ 
culating  the  first  several  VNAP  solution  lines.  The  Initial  Line  Method 
cannot  be  used  to  predict  the  specific  impulse,  because  a  thrust  calculation 
with  error  growth  should  not  be  ratioed  to  a  mass  flow  value  witli  no  error 
growth.  The  Exit  Line  Method,  however,  seems  to  give  a  consistent  prediction 
of  specific  impulse  for  all  9  cases.  This  is  because  the  mass  flow  and  thrust 
integrals  are  very  similar  and  experience  nearly  the  same  error  growth.  All 
nine  cases  are  found  to  yield  a  specific  impulse  of  434  +_  1.5  seconds.  Cor¬ 
recting  these  values  by  using  the  VNAP  correction  factor  (.978)  obtained 
during  the  dual  throat  analysis,  the  Mode  II  divergence  efficiency  was  found 
to  be  100%  +  .2%. 
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IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 

Although  the  results  for  Mode  II  operation  obtained  with 
VNAP  appear  to  be  qualitatively  correct,  it  is  clear  that  a  more  accurate 
method  must  be  used  to  predict  the  nozzle  performance  and  that  a  cold-flow 
test  program  is  needed  to  verify  these  results. 

2.  Performance  Calculations  For  Mode  I  Operation 

Calculating  the  dual  expander  performance  for  Mode  I  opera¬ 
tion  required  calculating  Regions  I,  II,  and  III  shown  in  Figure  I V- 17a . 
Results  obtained  for  each  region  are  discussed  below. 

Regions  I  and  II 

Results  of  the  Region  II  calculations  were  available  from  the 
Mode  II  analysis.  Region  II  results  were  discussed  in  the  preceding  subsec¬ 
tion  and  are  applicable  to  both  modes  of  operation.  Next,  Region  I  was  cal¬ 
culated  by  using  the  TDK  computer  program,  with  operating  conditions  for  the 
inner  nozzle  as  shown  below  in  Table  IV- VI I . 


Flow 

Mixture  Ratio 
Y 

CP 

Molecular  Weight 
Total  Temperature 
Total  Pressure 
Throat  Area 


TABLE  IV-VII 

MODE  I  OPERATING  CONDITIONS 
Outer  Nozzle  Inner  Nozzle 


L0X/LH2 

7.0 

1.18 

.8137  cal/gm°K 
16. 

3778°K  (6800°R) 

20.68  MPa  (3000  psia) 

280.6  cm^  (43.5  in. 2) 


LOX/RP-1 

2.8 

1.199 

.44985  cal/gm°K 

26.6 

3972°K  (7150°R) 

41.37  MPa  (6000  pcia) 

220.6  cm2  (34.2  in. 2) 
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IV,  C,  Dual  Expander  Flow  Field  Analysis  (cont.) 


It  was  intended  to  adjust  the  chamber  pressure  of  the  inner  nozzle  so  that  the 
calculated  exit  pressure  would  match  the  VNAP  exit  pressure.  However,  a 
satisfactory  match  was  obtained  on  the  first  attempt.  Pressures  calculated 
along  the  dividing  wall  for  Regions  I  and  II  are  shown  in  Figure  IV-23.  Pres¬ 
sures  at  the  downstream  end  of  the  wall  are  nearly  equal,  about  105  psia  on 
both  sides.  It  is  important  to  match  the  static  pressures  at  the  end  of  the 
wall  to  avoid  forming  a  strong  shock  wave  at  the  position  where  the  two  flow 
streams  come  together. 

Figure  IV-23  demonstrates  that  a  large,  negative  pressure 
gradient  occurs  along  the  upper  wall  and  that  a  much  lower  gradient  occurs 
along  the  lower  wall.  Thus,  if  the  wall  is  truncated,  a  large  difference  in 
static  pressure  will  result  at  the  end  of  the  wall.  For  example,  at  the 
"moderate"  truncation  position,  the  upper  and  lower  wall  pressures  are  found 
to  be  1970  (286)  and  979  (142)  KPa  (psia),  respectively.  For  "maximum"  trun¬ 
cation,  the  values  are  11606  (1683)  and  1248  (181)  KPa  (psia).  In  order  to 
match  the  pressures  at  the  lip  for  these  truncations  and  still  operate  at  the 
desired  chamber  pressure  requires  a  complete  redesign  of  the  nozzle. 

Region  III 

Region  III,  shown  in  Figure  I V- 1 7a ,  is  the  dual  expander  flow 
field  during  Mode  I  operation  (both  engines  on)  that  is  downstream  from 
Regions  I  and  II.  The  TDK  program  has  been  modified  to  calculate  Region  III. 
Previously,  TDK  would  calculate  the  flow  field  only  if  the  difference  between 
the  streamtube  and  stagnation  pressure  was  within  10%  of  each  other.  For  the 
dual  expander  this  difference  is  20%. 

Region  III  is  divided  into  an  inner  region  and  an  outer 
region.  The  inner  reyion  is  composed  of  exhaust  gases  from  the  L0X/RP-1 
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propel lants  at  a  mixture  ratio  of  2.8.  The  outer  region  is  composed  of 
exhaust  gases  from  the  LOX/LH2  propellants  at  a  mixture  ratio  of  7.0. 

Operating  conditions  for  the  two  engines  are  as  shown  in  Table  I V - V 1 1 . 

Inputs  along  the  start  line  for  the  inner  region  were  input 
from  the  exit  plant  of  the  TDK  calculation  for  Region  I.  Inputs  along  the 
start  line  for  the  outer  region  were  input  from  the  last  vertical  line  of  the 
VNAP  calculation  of  Region  II.  Values  input  are  pressure,  density,  flow 
angle,  and  velocity  at  each  point.  The  two  regions  are  divided  by  a  slip 
line,  i.e.,  a  double  valued  streamline  with  equal  values  of  pressure  and  flow 
angle,  but  different  values  of  velocity,  density,  Mach  number,  and  temperature 
on  each  side. 


The  TDK  calculations  were  made  by  removing  the  thickness  of 
the  dividing  wall  and  then  correcting  the  TDK  calculations  for  mass  flow  and 
thrust.  This  calculation  showed  the  divergence  efficiency  to  be  .995, 
matching  the  simplified  prediction. 

0.  RECOMMENDATIONS 

1 .  Mode  II  Performance  Prediction  Program 

After  careful  analysis  of  the  existing  Mode  II  performance 
prediction  methods  for  the  flow  field  of  the  dual  expander  propulsion  concept, 
it  is  now  felt  that  a  method  of  characteristics  analysis  and  computer  program 
capable  of  more  precise  prediction  techniques  is  needed.  A  sketch  illus¬ 
trating  the  geometry  and  the  major  features  of  the  flow  field  is  presented  in 
Figure  IV-24.  In  view  of  the  complexity  of  this  flow  field,  a  finite- 
difference  method  that  smears  Prandtl -Meyer  expansions  and  oblique  shock 
waves,  such  as  the  VNAP  program,  is  considered  inadequate  to  fully  determine 
the  potential  performance  of  the  dual  expander  concept. 


IV,  D,  Recommendations  (cont.) 


The  general  features  of  the  proposed  flow  model  are  pre¬ 
sented,  along  with  a  discussion  of  the  overall  numerical  algorithm.  The  pro¬ 
gram  would  be  less  complex,  faster  and  cheaper  to  run  than  VNAP.  Such  advan¬ 
tages  allow  the  program  to  be  run  parametr ical ly  to  determine  an  optimum  Mode 
II  design.  The  possibility  also  exists  to  expand  the  program  to  analyze  dual 
throat  Mode  II  operation. 

2.  Flow  Model 


The  flow  model  should  be  developed  for  steady  axisymmetric 
flew.  The  flow  should  be  allowed  to  be  rotational  to  account  for  the  entropy 
gradient  created  by  the  oblique  shock  wave  emanating  from  the  inner  wall.  The 
gas  would  be  assumed  to  have  constant  specific  heats  and  molecular  weight. 

The  supersonic  initial  value  line  (IVL)  would  be  input  to  the  analysis  in 
tabular  fonn.  Both  the  outer  wall  and  inner  wall  would  be  input  in  tabular 
form.  A  Prandtl -Meyer  expansion  would  emanate  from  point  L  on  the  inner  wall. 
A  discrete  left-turning  obliquie  shock  wave  would  emanate  from  point  N  on  the 
inner  wall  and  be  reflected  from  the  outer  wall  if  it  intersects  that  wall. 

All  unit  processes  for  calculating  the  flow  field  would  be 
based  on  the  method  of  characteristics.  The  interior  point  unit  process  would 
be  constructed  by  following  the  two  Mach  lines.  Outer  and  inner  wall  unit 
processes  would  follow  the  appropriate  Mach  line  and  the  streamlinfe  ( i . e . ,  the 
wall  itself). 


3.  Overall  Numerical  Algorithm 

The  overall  logic  would  follow  left-running  Mach  lines.  The 
first  left-running  line  would  emanate  from  the  second  point  on  the  initial- 
valve  line  (IVL)  below  the  outer  wall,  point  12.  That  Mach  line  would  be 
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extended  through  the  flow  field  until  it  intersects  the  outer  wall.  Subse¬ 
quent  left-running  Mach  lines  would  emanate  from  the  remaining  points  on  the 
IVL,  in  decreasing  order,  until  the  Mach  line  from  point  L  has  been  extended 
to  the  outer  wall.  Left-running  Mach  lines  would  then  be  generated  through 
the  Prandtl -Meyer  expansion  fan  at  point  L  until  the  required  amount  of  flow 
turning  at  point  L  has  been  achieved.  Each  of  these  left-running  Mach  lines 
would  be  extended  across  the  flow  field  to  the  outer  wall.  A  right-running 
Mach  line  would  then  be  extended  from  the  second  point  on  the  last  left- 
running  Mach  line  of  the  expansion  fan,  point  A,  to  the  inner  wall.  A  left¬ 
running  Mach  line  would  be  extended  from  that  wall  point  across  the  flow  field 
until  it  intersects  the  outer  wall.  This  procedure  would  be  repeated  until 
the  left-running  Mach  line  emanating  from  point  N  on  the  inner  wall  has  been 
extended  across  the  flow  field.  A  left-running  oblique  shock  wave  would  then 
be  generated  at  point  N  and  extended  across  the  flow  field.  This  shock  wave 
would  either  intersect  the  outer  wall  or  pass  beyond  the  nozzle  exit.  If  it 
intersects  the  outer  wall,  a  shock  reflection  unit  process  would  be  employed 
to  determine  the  flow  properties  across  the  reflected  shock  wave.  Left¬ 
running  Mach  lines  would  be  extended  from  points  on  the  inner  wall  downstream 
of  point  N  until  a  left-running  Mach  line  passes  beyond  point  K  on  the  outer 
wall  (the  exit  lip  point).  If  the  inner  wall  is  a  constant  radius  section, 
the  solution  is  complete.  If  not,  the  solution  would  be  continued  until  the 
flow  field  is  determined  on  the  inner  wall  at  point  W.  The  nozzle  thrust 
would  then  be  determined  by  calculating  the  momentum  flux  and  pressure  thrust 
across  the  tabular  initial-value  line,  line  IL,  and  adding  to  it  the  integral 
of  the  pressure  forces  along  the  outer  wall. 


PRELIMINARY  HEAT  TRANSFER  ANALYSIS 


The  objective  of  Task  IV  was  to  conduct  a  preliminary  thermal  analysis 
for  both  the  dual  throat  and  dual  expander  configurations.  The  dual  throat 
studies  represent  a  continuation  and  refinement  of  the  work  initiated  in  Ref. 
2.  Vehicle  studies  following  that  effort  resulted  in  a  change  in  the  base¬ 
line  design  point:  primary  and  secondary  chamber  pressures,  respectively, 
were  increased  from  20.68/14.48  MPa  (3000/2100  psia),  to  27.58/19.31  MPa 
(4000/2800  psia),  and  the  Mode  I  secondary  (LOX/RP-1)  thrust  was  increased 
from  60  to  70  percent  of  the  total.  Both  changes  make  chamber  cooling  more 
difficult  and  result  in  a  baseline  design  point  not  explicitly  considered  in 
the  parametric  studies  of  Ref.  2.  A  new  baseline  design  was  developed  which 
included  special  emphasis  on  the  lip  region  formed  by  the  termination  of  the 
primary  chamber  and  on  the  effect  of  primary  chamber  plume  impingement  in  the 
secondary  throat  region  during  Mode  II  operation.  New  parametric  design 
studies  investigated  the  effects  of  independent  changes  in  chamber  pressure 
and  Mode  I  thrust  split. 

Previous  thermal  analysis  of  the  dual  expander  concept  was  limited  to  a 
very  preliminary  I R&D  study  (Ref.  21)  which  indicated  problems  in  apportioning 
the  limited  hydrogen  coolant  flow  among  the  three  heated  surfaces.  Three  can¬ 
didate  baseline  designs  were  investigated  in  the  present  study:  (1)  regenera¬ 
tive  cooling  with  hydrogen  only;  (2)  regenerative  cooling  with  hydrogen  and 
oxygen;  (3)  a  transpiration-regenerative  cooling  combination  with  hydrogen.  A 
special  study  of  the  supersonic  lip  region  formed  by  the  termination  of  the 
primary  chamber  was  also  included  in  the  baseline  design  effort,  conducted  for 
a  primary  chamber  pressure  of  41.37  MPa  (6000  psia)  (L0X/RP-1),  secondary 
chamber  pressure  of  20.68  MPa  (3000  psia)  (LOX/H?).  and  a  Mode  I  primary 
thrust  equal  to  59  percent  of  the  698  MN  (600K  lbF)  total.  Parametric  studies 
investigated  the  effects  of  changes  in  thrust,  thrust  split,  and  both  chamber 
pressures  for  the  first  cooling  concept  noted  above. 
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V,  Preliminary  Heat  Transfer  Anlaysis  (cont.) 

A.  SPECIAL  STUDIES 
1.  Lip  Regions 

One  of  the  heat  transfer  problem  areas  unique  to  the  dual 
throat  and  dual  expander  concepts  is  the  lip  formed  by  the  termination  of  the 
primary  chamber.  During  Mode  I  operation,  the  primary  and  secondary  chamber 
flows  create  a  wake  region  behind  this  lip.  Therefore,  a  special  study  of  the 
lip  regions  was  conducted  to  identify  the  boundary  conditions  established  by 
the  wake  flow  and  the  design  problems  associated  with  regenerative  cooling  of 
a  corner  heated  on  both  sides. 

Wake  region  heat  transfer  data  in  the  literature  are  primar¬ 
ily  for  supersonic  flow  over  axisymmetric  bodies,  with  the  sidewall  boundary 
layer  thickness  small  in  relationship  to  the  base  diameter.  The  concept  of 
interest  to  us  incorporates  what  is,  essentially,  a  plane  surface  with  a  rela¬ 
tively  thick  boundary  layer  and  subsonic  flows  for  the  dual  throat  and  super¬ 
sonic  flows  for  the  dual  expander.  Because  of  these  differences,  a  conserva¬ 
tive  design  approach  has  been  adopted,  using  the  literature  data  as  a  guide. 
The  latter  indicate  that  heat  fluxes  at  the  center  of  the  wake  region  are 
typically  80-90  percent  of  the  sidewall  heat  flux  immediately  upstream  of  the 
wake  region  (see  Refs.  22  and  23),  but  that  a  significant  reduction  in  heat 
flux  is  observed  as  it  moves  away  from  the  center  of  the  wake  toward  the 
sidewall  (Refs.  22-24).  Figure  V-l  shows  the  radial  heat  flux  variation 
observed  in  Ref.  22;  Ref.  24  indicates  a  greater  reduction,  with  the  local 
base  heat  flux  being  proportional  to  (1  -  r/Rg)l*2  for  r/Rg  _<  0.8. 

Preliminary  mechanical  design  studies  indicate  that  a  conven¬ 
tional  rectangular  channel  construction,  as  illustrated  in  Figures  V- 2  and  V-3 
for  the  dual  throat  and  dual  expander,  respectively,  can  be  utilized  in  the 
lip  regions.  Therefore,  the  emphasis  in  the  present  study  was  on  the  lip 
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region  corners,  where  a  severe  cooling  problem  occurs  since,  at  a  land  center- 
line,  heat  input  from  two  surfaces  must  be  conducted  in  three  directions  to 
reach  cooled  surfaces.  A  three-dimensional  SINDA  model  of  such  a  corner  was 
developed  and  is  shown  in  Figure  V- 4.  The  model  does  not  include  any  part  of 
the  el ectroformed  nickel  closure  and  thus  is  limited  to  that  part  of  the 
corner  adjacent  to  the  coolant  channels.  Since  the  electroforined  nickel  can 
be  relatively  thin  in  the  lip  region,  the  corner  model  covers  a  significant 
fraction  of  the  lip  half-thickness.  Therefore,  the  procedure  used  in  the  pre¬ 
sent  study  to  provide  a  conservative  design  relative  to  the  boundary  condi¬ 
tions  of  Figure  V-l  was  as  follows: 

(1)  Apply  a  unifonn  heat  transfer  coefficient  over  the  end 
wall  equal  to  80  percent  of  the  sidewall  coefficient. 

(2)  Set  the  lip  channel  depth  equal  to  the  sidewall  channel 
depth  (symmetrical  corner  model),  thereby  accommodating 
an  even  higher  heat  flux  at  the  center  of  the  lip. 

It  should  be  noted  that  the  sidewall  boundary  layer  was  assumed  to  be 
unaffected  by  the  throat  region  shock  structure.  In  addition,  the  coolant 
heat  transfer  coefficients  were  considered  to  be  unaffected  by  the  turn  at  the 
lip  corner;  the  actual  enhancement  due  to  curvature  was  omitted  to  compensate 
for  the  lack  of  a  channel  corner  radius  and  increased  wall  thickness  in  the 
wall  conduction  model. 

The  corner  model  was  used  to  define  the  lip  region  channel 
depth  required  for  each  baseline  design  to  limit  the  corner  temperature  to 
approximately  756°K  (900°F).  Figure  V-5  shows  a  typical  gas-side  wall  temper¬ 
ature  distribution  over  a  land  centerline  in  the  vicinity  of  a  lip  corner  for 
the  dual  throat  concept;  temperatures  over  the  channel  centerline  are  only 
slightly  lower. 
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Figure  V-4.  Lip  Region  3-D  Thermal  Model 
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V-5.  Dual  Throat  Lip  Region  Wall  Temperatures 
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In  order  to  provide  lip  region  channel  depths  for  the  para¬ 
metric  design  studies,  we  used  a  simple  one-dimensional  conduction  model  with 
an  effective  wall  conductance  and  coolant  surface  area  defined  from  results  of 
the  three-dimens ional  SINDA  model.  The  validity  of  this  model  was  demon¬ 
strated  by  successfully  predicting  a  dual  throat  baseline  design,  using  a  dual 
expander  design  to  define  the  one-dimensional  model  parameters.  In  all  cases, 
the  total  pressure  drop  for  the  two  lip  region  corners  was  taken  as  three 
velocity  heads.  Results  from  the  3D  baseline  and  ID  parametric  models  are 
presented  in  subsequent  sections  as  part  of  the  overall  results  for  each 
concept. 


2.  Secondary  Throat  Analysis  for  Dual  Throat  Mode  II 

During  normal  Mode  II  operation  of  the  dual  throat  concept, 
the  secondary  chamber  bleed  flow  causes  the  primary  chamber  plume  to  attach  to 
the  secondary  wall  downstream  of  the  throat  such  that  the  plume  flow  is  nearly 
parallel  to  the  wall.  In  this  case,  throat  heat  fluxes  are  expected  to  be 
well  below  those  during  Mode  I  operation.  Of  interest  herein  is  the  possi¬ 
bility  of  bleed  flow  failure  and  impingement  of  the  primary  plume  upstream  of 
the  secondary  throat,  with  resultant  heat  fluxes  that  might  be  more  severe 
than  Mode  I  val ues. 

A  number  of  studies  of  rocket  exhaust  plume  impingement  on 
adjacent  vehicle  structures  were  found  in  the  literature  (Refs.  25-27). 
Analyses  and  correlation  of  these  data  split  the  problem  into  strong  and  weak 
shock  regions,  with  the  former  treated  as  a  stagnation  point  region  as  far  as 
heat  transfer  is  concerned.  However,  because  of  the  effect  on  the  plume 
boundary  of  the  back  pressure  in  the  secondary  chamber,  formation  of  a  strong 
shock  region  was  not  considered  likely  in  the  present  problem.  Therefore,  the 
VNAP  program  was  used  to  define  the  inviscid  flow  field  for  zero  bleed  flow. 
This  program  provides  an  exact  inviscid  solution,  and  its  use  for  the  nominal 
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bleed  flow  case  was  already  planned  as  part  of  Task  I.  A  boundary  layer 
analysis  based  on  the  edge  conditions  determined  by  VNAP  was  used  to  define 
the  heat  transfer  to  the  secondary  wall. 

Figure  V- 6  shows  the  Mode  II  pressure  and  mass  velocity  pro¬ 
files  at  the  secondary  wall  from  VNAP  for  the  baseline  design  with  zero  bleed 
flow.  The  maximum  pressure  and  mass  velocity  due  to  the  plume  impingement 
occurs  4.3  an  (1.7  in.)  upstream  of  the  throat.  The  freestream  velocity, 
shown  in  Figure  V - 7 ,  is  essentially  constant  in  the  region  of  interest,  i.e., 
from  the  attachment  point  through  the  zone  of  maximum  heating.  This  feature  of 
the  inviscid  flow  was  used  to  advantage  herein,  since  it  greatly  simplifies 
the  integral  momentum  equation  for  the  boundary  layer.  For  constant  velocity, 
this  equation  becomes 

!?.  +  1  0  =  C,/2 

dx  oe  dx  T 

This  form  of  the  momentum  equation  was  solved  by  using  the  Van  Driest  skin 
friction  law  for  compressible  flow  based  on  the  Von  Karman  turbulent  mixing 
length  (Ref.  28): 

ReQ  =  0.191  — - 

"e 

in  which  is  a  function  of  Mach  number  and  wall  temperature  ratio.  The  Mach 
number  and  freestream  static  temperature  distributions  from  VNAP  are  nearly 
uniform.  A  design  wall  temperature  of  811°K  (1460°R),  consistent  with  the 
Mode  I  design  studies  described  later,  was  used. 

The  momentum  thickness  at  the  attachment  point  must  be  spe¬ 
cified  to  start  the  above  boundary  layer  analysis.  Two  cases  were  considered. 
In  one  case,  the  momentum  thickness  of  the  plume  shear  layer  defined  by  the 
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Task  I  bleed  flow  aerodynamic  model  was  used.  This  approach  assumes  that  the 
shear  layer  is  not  disturbed  by  its  interaction  with  the  wall.  In  the  second 
case,  the  shear  layer  was  assumed  to  be  nearly  destroyed  by  its  interaction 
with  the  wall,  such  that  the  initial  momentum  thickness  for  the  wall  boundary 
layer  is  very  small.  Results  reported  for  this  case  are  for  an  initial  thick¬ 
ness  of  .0254  mm  (0.001  in.),  compared  with  a  shear  layer  momentum  thickness 
of  .861  mm  (0.0339  in.).  Further  reductions  in  initial  momentum  thickness 
have  little  effect  on  the  Stanton  number  at  the  location  of  the  maximum  free- 
stream  mass  velocity.  Therefore,  these  two  cases  are  considered  to  bound  the 
actual  interaction  of  the  shear  layer  with  the  wall. 

Stanton  numbers  were  defined  from  the  skin  friction  coeffi¬ 
cient,  using  a  Reynolds  analogy  factor  of  1.25.  This  value  is  a  compromise 
between  the  value  of  1.19,  calculated  from  the  theory  of  Ref.  28  for  a  turbu¬ 
lent  Prandtl  number  of  0.9,  and  the  value  of  1.32  defined  by  the  Colburn 
factor  Pr~-/3.  Results  from  the  boundary  layer  analysis  are  summarized 
below  for  the  location  of  maximum  freestream  mass  velocity: 

Initial  Momentum  Thickness,  mm  (in.) 

.861  (0.0339)*  .0254  (.001) 

0,  .024 7  .00244 

Refl  x  10-4  3.35  .331 

St  .00179  .00263 

hg,  KW/m2-°K  30.32  (.0103)  44.45  (.0151) 

(Btu/in.2-sec-°F) 

Heat  Flux,  MW/m2  86.1  (52.7)  126.3  (77.3) 

(Btu/i n. 2-sec) 

♦Defined  by  plume  shear  layer 

These  results  compare  with  a  maximum  Mode  I  heat  transfer 
coefficient  of  28.38  KW/m2-°K  (0.00964  Btu/i n.2-sec-°F)  and  heat  flux 
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of  85.46  MW/m2  (52.3  Btu/ i n. 2-sec )  in  the  secondary  throat  region.  Thus, 
when  the  initial  momentum  thickness  for  the  zero  bleed  flow  boundary  layer  is 
equal  to  that  from  the  plume  shear  layer,  the  maximum  Mode  1  and  Mode  II  heat 
fluxes  are  essentially  equal.  A  similar  result  for  Mode  II  is  obtained  by 
scaling  the  maximum  heat  flux  from  the  Task  I  BLIMP  run  for  nominal  bleed  flow 
according  to  the  local  VNAP  wall  pressures  (a  procedure  suggested  by  various 
flow  reattachment  data): 

0  ,  =  A  /  Pno  bleed  ] 

no  bleed  *  nom.  bleed  P  .  I 

nom.  bleed / 

BLIMP  \  VNAP  / 

_397\  7 

98. 9j  57  Btu/ in.  -sec  (93.14  MW/m2) 

The  nominal  uleed  flow  BLIMP  was  started  with  an  initial  momentum  thickness 
from  the  corresponding  plume  shear  layer  (D  streamline)  since  the  plume  flow 
is  very  nearly  parallel  to  the  wall.  When  the  initial  momentum  thickness  for 
the  zero  bleed  flow  boundary  layer  is  very  small,  the  maximum  Mode  II  heat 
flux  is  almost  50%  higher  than  that  for  Mode  I.  The  ability  of  the  secondary 
chamber  cooling  system  to  accommodate  this  higher  heat  flux  is  discussed  in 
Section  V,C,1. 

B.  CHANNEL  DESIGN  PROCEDURES 

Two  computer  programs  for  regenerative  cooling  were  used  for  the 
[resent  studies:  one  with  the  dual  throat  geometric  details  built  in,  the 
other  for  the  dual  expander.  These  programs  were  derived  from  the  SCALER 
program  for  conventional  geometry  which  was  developed  specifically  for  para¬ 
metric  design  studies.  With  these  programs,  it  is  economically  feasible  to 
generate  a  relatively  large  number  of  parametric  design  points  and  still 
obtain  a  detailed,  multi-station  analysis  of  a  rectangular  channel  at  each 
design  point.  The  chamber  geometry  and  the  local  gas-side  heat  transfer 
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coefficients  and  coolant  heat  loads  are  scaled  from  reference  input  to  other 
thrusts  and  chamber  pressures.  At  each  station,  the  programs  iterate  to 
determine  the  channel  depth  required  to  satisfy  a  gas-side  wall  temperature 
limit,  which  can  be  specified  as  a  function  of  closeout  wall  temperature 
consistent  with  cycle  life  and  creep  criteria  (Figure  V-8).  The  only  simpli¬ 
fying  assumption  is  that  gas-side  wall  temperature  differences  between  the 
reference  input  and  scaled  cases  have  a  negligible  effect  on  gas-side  heat 
transfer  coefficients  and  heat  loads.  Normally,  gas-side  wall  temperature 
limits  are  well  known  in  advance  so  that  local  reference  gas-side  heat 
transfer  analyses  can  be  run  at  appropriate  wall  temperatures. 

Two-dimensional  conduction  effects  in  rectangular  channels,  as 
well  as  the  spatial  variations  of  the  coolant  heat  transfer  coefficient,  are 
approximated  as  shown  in  Figure  V-9.  The  hot  wall  adjacent  to  the  channel, 
the  entire  land,  and  that  part  of  the  external  wall  adjacent  to  the  channel 
are  represented  as  fins;  that  part  of  the  external  wall  adjacent  to  the  land 
is  assumed  to  be  isothermal.  A  simple  resistance  network  represents  the  hot 
wall  adjacent  to  the  land  and  connects  the  hot-wall  fin  to  the  land  fin.  Dif¬ 
ferent  coolant  heat  transfer  coefficients  are  associated  with  each  fin.  That 
for  the  hot  wall  is  based  on  the  average  of  the  centerline  ( TWL2 )  and  corner 
(TCORN)  wall  temperatures,  while  the  coefficient  for  the  external  wall  is 
based  on  its  centerline  wall  temperature  (TBS).  The  land  coefficient  is 
merely  a  weighted  average  of  the  other  two.  In  this  way,  it  has  been  possible 
to  simulate  the  wall  temperatures  predicted  by  the  SINDA  network  analyses  very 
closely. 


Channel  geometry  constraints  imposed  in  the  present  study  due  to 
fabrication  and  flow  distribution  considerations  were  as  follows: 


168 


a# 


V,  B,  Channel  Design  Procedures  (cont.) 


Minimum  wall  thickness,  mm  (in.) 
Minimum  land  width,  mm  (in.) 
Minimum  channel  width,  mm  (in.) 
Minimum  channel  depth,  mm  (in.) 
Maximum  channel  depth/width 


.635  (.025) 
1.016  (.040) 
1.016  (.040) 
.889  (.035) 
5:1 


Wall  thicknesses  were  also  required  to  satisfy  the  strength  criteria  of  Figure 
V-10.  The  number  of  coolant  channels  was  determined  by  utilizing  the  above 
minimum  dimensions  for  both  the  land  and  channel  widths  at  either  the  throat 
or  the  location  in  a  circuit  of  the  minimum  radius. 


Channel  layouts  were  limited  to  combinations  of  regions  with  con¬ 
stant  channel  width  or  constant  land  width  (straddle  milling).  Increases  in 
flow  area  relative  to  the  throat  beyond  those  possible  with  channel  depth 
increases  were  accomplished  by  constant  land  width  regions  or  by  step  changes 
in  channel  width.  At  the  pressure  levels  associated  with  the  present  baseline 
designs,  only  very  small  increases  in  channel  width  .102-. 127  mmg  (.004  -  .005 
in.)  above  the  1.016  mm  (0.040  in.)  minimum  were  possible  without  increasing 
the  wall  thickness  above  .635  mm  (0.025  in.).  Details  of  the  channel  layouts 
for  each  concept  are  presented  in  subsequent  sections. 

C.  DUAL  THROAT 

1.  Baseline  Design 

As  shown  schematically  in  Figure  V-ll,  three  parallel  cooling 
circuits  are  used  in  all  dual  throat  designs  presented  herein.  A  separate 
circuit  is  included  for  the  inner  annulus,  lip  region,  and  part  of  the  primary 
nozzle;  in  previous  dual  throat  analyses  (Ref.  20),  these  regions  were 
included  in  the  primary  chamber  circuit.  Since  the  primary  chamber  dictates 
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Dual  Throat  Coolant  Flow  Schematic 


V,  C,  Dual  Throat  (cont.) 

the  overall  system  pressure  drop,  use  of  a  separate  circuit  to  accommodate  the 
pressure  drop  required  for  lip  cooling  significantly  reduces  the  overall  pres¬ 
sure  drop.  The  primary  nozzle  area  ratio,  £j,  at  the  interface  between 
these  circuits  was  a  parameter  in  the  baseline  design  study;  a  maximum  value 
of  1.92  is  possible  from  a  mechanical  design  standpoint.  Note  that  the 
secondary  circuit  extends  to  an  area  ratio  of  7.8,  after  which  an  oxygen- 
cooled  tube  bundle  would  be  used  for  the  remainder  of  the  secondary  nozzle.  A 
parametric  design  study  of  such  a  tube  bundle  was  presented  in  Ref.  2. 

Channel  layout  details  for  the  baseline  design  are  given  in 
Table  V-I.  The  primary  circuit  utilizes  straddle  milling  in  part  of  the  con¬ 
vergent  section  in  order  to  increase  the  channel  width  from  the  throat  values. 
In  order  to  avoid  the  channel  aspect  ratio  limit  in  the  barrel,  it  was  neces¬ 
sary  to  increase  the  channel  width  above  that  consistent  with  the  strength 
criteria  for  the  minimum  wall  thickness.  Therefore,  the  barrel  wall  thickness 
was  increased  to  .762  mm  (0.030  in.).  In  the  secondary  circuit,  the  minimum 
wall  thickness  of  .635  mm  (0.025  in.)  was  used  in  the  barrel,  which  limited 
the  channel  width  to  1.143  mm  (0.045  in.),  i.e.,  only  .127  rmi  (0.005  in.) 
greater  than  the  throat  channel  width.  In  order  to  increase  the  flow  area 
away  from  the  throat  region  more  rapidly  than  provided  for  by  straddle 
milling,  a  step  change  in  width  occurs  in  the  convergent  section.  This  barrel 
design  was  not  optimized  for  Mode  I,  since  throat  cooling  in  Mode  II  dictates 
most  of  the  circuit  pressure  drop  and  since  the  total  pressure  drop  does  not 
appear  to  be  a  problem.  The  throat  land  width  is  maintained  in  the  first  part 
of  the  secondary  nozzle,  with  the  wall  thickness  increasing  to  1.270  mm  (0.050 
in.)  to  accommodate  the  increased  channel  widths.  A  channel  width  consistent 
with  the  strength  criteria  for  the  1.270  mm  (0.050  in.)  wall  thickness  is 
maintained  in  the  second  part  of  the  secondary  nozzle.  A  uniform  channel 
width  of  1.016  mi..  (0.040  in.)  was  maintained  in  the  lip  circuit.  This  is 
necessary  in  the  primary  nozzle  segment  of  the  circuit  in  order  to  provide  the 
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most  effective  channel  geometry  in  the  lip  region.  However,  the  inner  annulus 
segment  could  be  straddle-milled,  which  would  provide  a  small  reduction  in 
circuit  pressure  drop. 

Figure  V— 1 2  shows  the  pressure  drops  required  for  the  primary 
circuit  for  two  inlet  area  ratios  as  a  function  of  the  fraction  of  the  hydro¬ 
gen  flow  used  in  the  circuit.  Lach  point  on  these  curves  represents  a  dif¬ 
ferent  channel  depth  profile,  as  defined  by  the  cycle  life  criterion  of  Figure 
V-8  and  the  maximum  depth/width  ratio  of  5.  As  the  hydrogen  flow  fraction 
initially  increases,  the  local  channel  depth  increases  for  a  fixed  coolant 
velocity;  however,  the  increased  depth  results  in  a  larger  effective  coolant 
surface  area,  thereby  reducing  the  required  coolant  velocity  and  allowing  an 
even  greater  channel  depth.  The  reduced  velocity  and  increased  hydraulic 
diameter  result  in  a  significant  decrease  in  required  pressure  drop.  However, 
as  the  coolant  flow  continues  to  increase,  substantial  portions  of  the  chamber 
reach  the  channel  aspect  ratio  limit  resulting  in  a  pressure  drop  increase. 
Figure  V - 1 2  indicates  that  reducing  the  inlet  area  ratio  of  the  primary  cir¬ 
cuit  reduces  its  pressure  drop;  the  coolant  Mach  number,  which  exceeds  0.30 
for  higher  area  ratios,  is  also  reduced.  In  addition,  the  hydrogen  flow  frac¬ 
tion  at  the  minimum  pressure  drop  decreases  slightly  as  the  inlet  area  ratio 
is  reduced,  thereby  providing  more  flow  for  the  other  circuits. 

Figure  V-l?  also  includes  the  required  pressure  drop  charac¬ 
teristics  for  the  secondary  circuit  for  Mode  I  operation.  A  hydrogen  flow 
fraction  of  0.30  was  selected  for  the  secondary  circuit  in  order  to  provide  a 
large  pressure  drop  margin  for  redesign  of  the  throat  region  for  Mode  II  oper¬ 
ation  with  a  bleed  flow  malfunction.  The  Mode  II  secondary  throat  analysis 
discussed  previously  (Section  V,A,2)  indicated  that  local  heat  fluxes  in  the 
primary  plume  impingement  region  may  be  as  much  as  50%  higher  than  for  Mode  I. 
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Figure  V-12.  Dual  Throat  Baseline  Cooling  Circuit  Pressure  Drop  Characterisi tcs 
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A  composite  design  with  revised  channel  depths  in  the  Mode  II  impingement 
region  was  developed  for  the  worst-case  boundary  condition  of  Section  V ,  A  ,2 ; 
Figure  V - 1 3  compares  the  Mode  I  and  II  heat  flux  profiles.  The  pressure  drop 
for  the  composite  design  during  Mode  I  operation  is  4.45  MPa  (645  psi),  which 
is  well  below  the  primary  circuit  pressure  drop  of  Figure  V- 1 2 .  A  detailed 
design  of  the  lip  region  circuit  was  completed  for  an  interface  area  ratio  of 
1.58,  assuming  a  primary  circuit  flow  fraction  equal  to  the  optimum  value  of 
0.40.  In  all  cases,  the  total  hydrogen  flow  fraction  available  for  cooling 
was  assumed  to  be  0.94,  based  on  the  gas  yenerator/staged  combustion  mixed 
cycle  studies  of  Ref.  2.  Therefore,  the  flow  fraction  available  for  the  lip 
region  circuit  was  0.24.  Table  V - 1 1  summarizes  the  baseline  designs  for  all 
three  cooling  circuits.  Since  the  lip  circuit  pressure  drop  is  well  below 
that  of  the  primary  circuit,  it  is  apparent  that  further  reduction  in  t fie 
interface  area  ratio  is  desirable  in  order  to  minimize  the  cooling  system 
pressure  drop. 

2.  Parametric  Study 

Chamber  pressure  and  thrust  split  were  varied  independently 
in  the  dual  throat  parametric  study.  Primary  chamber  pressure*  was  varied 
from  20.68  to  34.47  MPa  (3000  to  5000  psia)  at  the  baseline  thrust  split. 
However,  at  34.47  MPa  (5000  psia),  the  coolant  Mach  number  was  impractical; 
therefore,  the  highest  pressure  for  which  results  are  reported  is  31.06  MPa 
(4500  psia).  Mode  I  thrust  fror  the  secondary  propellants  was  varied  from  60 
to  80%  at  the  baseline  chamber  iressures.  The  design  procedure  in  each  case 
was  similar  to  that  for  the  baseline,  except  that  composite  secondary  designs 
including  the  effects  of  Mode  11  were  not  generated.  Primary  and  secondary 
circuit  designs  for  Mode  I  were  defined  as  a  function  of  hydrogen  flow 

♦Secondary  chamber  pressure  was  always  70%  of  the  primary,  as  defined  by 
aerodynamic  considerations. 
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Figure  V-13.  Comparison  of  Secondary  Throat  Region  Heat  Fluxes 
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fraction,  as  shown  previously  In  Figure  V-12  for  the  baseline.  The  secondary 
circuit  flow  fraction  was  selected  to  provide  a  large  pressure  drop  margin 
relative  to  the  minimum  primary  circuit  pressure  drop;  this  margin  is  avail¬ 
able  to  provide  smaller  channel  depths  in  the  Mode  II  impingement  region.  The 
remaining  coolant  flow  was  divided  between  the  primary  and  lip  region  cir¬ 
cuits,  with  the  primary  flow  selected  to  give  the  minimum  pressure  drop 
whenever  possible.  For  the  higher  thrust  split  and  chamber  pressure  cases,  it 
was  necessary  to  select  a  primary  flow  less  than  optimum,  resulting  in  a  pres¬ 
sure  drop  up  to  11$  higher  than  minimum.  Due  to  the  time  required  to  design 
the  lip  region  circuit,  the  flow  split  between  the  primary  and  lip  circuits 
was  not  optimized  by  balancing  the  required  pressure  drops. 

The  parametric  designs  are  summarized  in  Table  V-III. 

Figures  V-14  and  V-15  show  the  coolant  pressure  drop  and  Mach  number  as  a 
function  of  primary  chamber  pressure  and  thrust  split,  respectively.  At  the 
highest  pressure  and  thrust  split,  the  coolant  Mach  number  has  reached  0.37; 
therefore,  further  increases  in  these  parameters  are  not  practical  with  the 
present  channel  layout. 

D.  DUAL  EXPANDER 

1.  Baseline  Designs 

As  shown  in  Figure  V-16,  the  initial  design  considered  was  a 
hydrogen  regeneratlvely  cooled  system  with  three  parallel  cooling  circuits. 
Note  that  these  circuits  are  conceptually  the  same  as  for  the  dual  throat 
designs.  Although  the  lip  region  cooling  problem  is  less  severe  for  the  dual 
expander,  the  inner  annulus  circuit  is  much  longer  and  includes  a  throat 
region. 


TABLE  V-III 

DUAL  THROAT  PARAMETRIC  DESIGN  SUMMARY 
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Channel  layout  details  for  the  baseline  design  are  included 
in  Table  V-I.  The  primary  circuit  layout  is  similar  to  the  dual  throat  secon¬ 
dary  circuit.  Channel  aspect  ratios  in  the  barrel  are  within  the  design 
limit,  so  there  was  no  motivation  for  increasing  the  wall  thickness  to  provide 
a  wider  channel.  A  uniform  channel  width  of  1.016  mm  (0.040  in.)  was  main¬ 
tained  in  the  barrel  and  convergent  sections  of  the  secondary  circuit  because 
of  the  small  channel  depths  required.  A  uniform  channel  width  was  also  used 
over  the  entire  inner  annulus  surface.  It  was  found  that  a  1.118  mm  (0.044 
in.)  width  provided  a  lower  pressure  drop  than  the  minimum  allowable  width  of 
1.016  mm  (0.040  in.).  The  latter  value  was  maintained  in  the  primary  nozzle 
section  of  the  inner  annulus  circuit.  This  section  defined  the  number  of 
channels  in  the  circuit  and  is  shown  separately  in  Table  V-I. 

Two  area  ratios  were  considered  for  the  interface  between 
the  primary  circuit  inlet  and  the  inner  annulus  circuit  outlet:  2.34,  which 
represents  the  maximum  from  a  mechanical  design  standpoint,  and  1.66.  Figure 
V-17  shows  the  primary  circuit  pressure  drop  characteristics  for  these  two 
interface  area  ratios,  along  with  the  secondary  circuit  characteristics.  The 
coolant  flow  fraction  for  the  latter  is  set  at  0.24  by  the  minimum  channel 
depth  constraint.  Inner  annulus  flow  requirements  for  Mach  number  control, 
the  remaining  hydrogen  flow  fractions  available  for  the  primary  circuit,  and 
the  corresponding  primary  circuit  pressure  drops  from  Figure  V-17  are  as 
follows: 


Interface 
Area  Ratio 


Inner  Annulus 
Flow  Fraction 


Primary 
Flow  Fraction 


Primar; 
Ap.  MPa 


'U) 


11.86  (1720) 

12.86  (1860) 


The  reduced  flow  available  for  the  primary  circuit  with  the  lower  interface 
area  ratio  offsets  the  reduced  length  and  increases  the  primary  circuit 
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pressure  drop.  In  both  cases,  the  inner  annulus  flow  required  for  coolant 
Mach  number  control  forces  the  primary  flow  to  be  less  than  optimum.  Table 
V-IV  summarizes  the  circuit  designs  for  an  interface  area  ratio  of  2.34.  In 
view  of  the  Mach  numbers  in  the  primary  and  inner  annulus  circuits  and  the 
channel  depth  in  the  secondary  circuit,  it  is  clear  that  cooling  limitations 
have  been  reached  in  both  the  primary  and  secondary  chambers.  The  secondary/ 
primary  chamber  pressure  ratio  was  selected  in  anticipation  of  such  a  result. 
However,  a  slightly  lower  ratio  would  have  allowed  the  primary  circuit  pres¬ 
sure  drop  to  be  optimized. 

Figure  V-17  indicates  that  the  i\ydrogen  pressure  drop  can  be 
reduced  to  10.69  MPa  (1550  psi)  by  increasing  the  primary  circuit  flow  frac¬ 
tion  to  0.40  with  an  interface  area  ratio  of  1.66.  This  can  be  accomplished 
by  using  oxygen  to  cool  the  secondary  circuit.  Table  V-V  summarizes  such  a 
design;  it  should  be  noted  that  the  secondary  circuit  channel  design  has  not 
been  optimized  for  oxygen  cooling. 

Figure  V-18  shows  the  effect  of  transpiration  cooling  the 
throat  region  of  the  primary  chamber  for  an  interface  area  ratio  of  1.66.  The 
upstream  end  of  the  transpiration  region  was  varied,  with  the  downstream  end 
fixed  at  area  ratio  1.44,  which  is  3.74  cm  (1.47  in.)  aft  of  the  throat. 

Since  the  secondary  chamber  was  assumed  to  be  hydrogen  cooled,  only  30%  of  the 
total  hydrogen  flow  was  available  for  the  two  primary  cooling  circuits.  Con¬ 
sequently,  the  significant  coolant  flow  fractions  required  for  the  longer 
transpiration  sections  severely  limit  the  available  regenerative  cooling  flow 
and  result  in  increased  pressure  drops.  The  minimum  regenerative  pressure 
drop  in  this  case  is  about  11.03  MPa  (1600  psi),  which  is  higher  than  the 
optimum  value  of  Figure  V-17  that  can  be  achieved  by  oxygen  cooling  of  the 
secondary. 


UAL  EXPANDER  BASELINE  COOLINS  CIRCUITS 
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DUAL  EXPANDER  BASELINE  WITH  OXYGEN-COOLED  SECONDARY 
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Figure  V-18.  Dual  Expander  Transpiration  Cooling 
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2.  Parametric  Study 

The  dual  expander  parametric  study  was  based  on  the  hydrogen 
regeneratively  cooled  concept  of  Table  V- IV.  Independent  variations  of  thrust 
and  Mode  I  thrust  split  at  the  baseline  chamber  pressures  were  investigated, 
along  with  various  chamber  pressure  combinations  at  the  baseline  thrust  and 
thrust  split.  Table  V-VI  indicates  the  design  points  attempted  and  identifies 
those  for  which  chamber  designs  could  not  be  obtained;  characteristics  of  the 
successful  designs  are  summarized  in  Table  V-VI  I . 

In  all  cases  with  increased  chamber  pressures  or  Mode  I 
thrust  split  relative  to  the  baseline,  designs  could  not  be  obtained  due  to 
insufficient  hydrogen  flow  for  the  cooling  of  three  circuits.  A  reduction  in 
secondary  chamber  pressure  to  10.34  MPa  (1500  psia)  with  the  primary  chamber 
pressure  unchanged  resulted  in  convergence  failure  due  to  large  pressure  drops 
in  the  primary  circuit.  The  latter  results  from  the  reduction  in  coolant 
inlet  pressure  associated  with  the  reduction  in  secondary  chamber  pressure. 
From  the  baseline  and  parametric  results,  it  appears  that  the  present  baseline 
design  point  is  near  the  optimum  chamber  pressure  ratio  and  is  very  close  to 
the  chamber  pressure  vs  thrust  split  characteristic  which  represents  a  prac¬ 
tical  cooling  limit.  Future  parametric  studies  should  define  additional 
points  on  such  a  characteri st ic  and  consider  smaller  chamber  pressure  incre¬ 
ments  to  optimize  the  pressure  ratio. 

E.  SUMMARY  OF  REC  ".TS 

A  satisfactory  dual  throat  baseline  design  was  obtained  by  using 
three  parallel  cooling  circuits.  The  primary  circuit  coolant  flows  from  area 
ratio  1.58  in  the  primary  nozzle  to  the  primary  injector;  this  circuit  defines 
the  system  pressure  drop,  so  its  coolant  flow  fraction  is  selected  to  provide 
the  minimum  pressure  drop.  A  separate  circuit  cools  the  inner  annulus 
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surface,  the  lip  region,  and  the  rest  of  the  primary  nozzle.  The  majority  of 
the  channel  pressure  drop  In  this  circuit  is  required  for  the  lip  region,  for 
which  a  three-dimensional  SINDA  model  was  developed.  A  counterflow  circuit 
cools  the  outer  surface  of  the  secondary  chamber.  This  circuit  is  a  composite 
design  considering  both  Mode  I  and  Mode  II  operation;  the  channel  design  over 
a  short  region  upstream  of  the  secondary  throat  is  defined  by  the  primary 
plume  impingement  associated  with  a  bleed  flow  malfunction  in  Mode  II.  An 
exact  solution  for  the  invlscid  plume  impingement  problem  was  obtained  with 
the  VMAP  program.  This  solution  provided  edge  conditions  for  an  integral 
boundary  layer  analysis  which  defined  the  wall  heat  transfer.  The  initial 
momentum  thickness  for  the  boundary  layer  analysis  was  varied  in  order  to 
bound  the  problem,  with  the  channel  design  based  on  the  worst  case;  the 
maximum  heat  flux  ranged  from  essentially  the  Mode  I  value  to  almost  501 
higher. 


Table  V-I I  summarizes  individual  cooling  circuit  characteristics 
for  the  dual  throat  baseline  design.  Although  the  fuel  system  pressure  drop 
of  8.00  MPa  (1160  psl)  is  set  by  the  primary  circuit,  the  coolant  Mach  number 
is  slightly  higher  in  the  secondary  circuit  when  the  latter  is  designed  to 
accommodate  the  maximum  postulated  Mode  II  heat  flux  with  no  bleed  flow. 
Parametric  study  results  indicate  that  it  is  the  coolant  Mach  number  rather 
than  pressure  drop  which  limits  increases  in  chamber  pressure  and  Mode  I 
thrust  split.  Increasing  the  chamber  pressures  to  31.02/21.72  MPa  (4500/3150 
psl a)  or  the  thrust  split  to  80/20  results  in  a  primary  circuit  Mach  number  of 
0.37. 


The  dual  expander  designs  are  cooled  with  three  parallel  circuits 
similar  to  those  of  the  dual  throat.  However,  the  inner  annulus  circuit  is 
now  a  major  circuit  in  terms  of  length  and  heat  load,  and  it  includes  a  throat 
region.  Regenerative  cooling  of  the  baseline  design  point  with  hydrogen  only 
is  difficult.  Coolant  flows  in  the  secondary  and  inner  annulus  circuits  must 
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be  selected  to  avoid  minimum  channel  depth  and  Mach  number  limits,  respec¬ 
tively,  and  the  remaining  flow  available  for  the  primary  is  less  than  optimum. 
As  a  result,  the  system  pressure  drop  of  11.86  MPa  (1720  psi)  is  .896  MPa  (130 
psl)  above  the  minimum  primary  circuit  value.  Using  oxygen  to  cool  the  secon¬ 
dary  circuit  and  reducing  the  interface  area  ratio  between  the  primary  circuit 
inlet  and  the  inner  annulus  circuit  outlet  allows  the  hydrogen  pressure  drop 
to  be  reduced  to  10.69  MPa  (1550  psi).  These  designs  are  sumamrized  in  Tables 
V-IV  and  V-V.  Transpiration  cooling  the  throat  region  of  the  primary  chamber 
was  not  effective  since  the  hydrogen  flow  fractions  required  for  transpiration 
cooling,  coupled  with  the  flow  requirements  of  the  other  circuits,  severely 
limit  the  flow  available  for  regenerati vely  cooling  the  remainder  of  the 
primary  circuit  and  result  in  relatively  high  pressure  drops. 

Parametric  study  designs  could  not  be  obtained  with  Increased 
chamber  pressures  or  increased  Mode  I  thrust  split,  or  with  the  secondary 
chamber  pressure  reduced  to  10.34  MPa  (1500  psia).  The  latter  results  from 
the  high  primary  chamber  pressure  drop  associated  with  reduced  coolant  Inlet 
pressure.  These  results,  along  with  the  baseline  design  details,  indicate  the 
baseline  design  point  is  near  the  optimum  chamber  pressure  ratio  and  is  very 
close  to  the  chamber  pressure  vs  thrust  split  characteristic  which  represents 
a  practical  cooling  limit. 


VI. 


HOT-FIRE  DATA  ANALYSIS 


A.  PRE-TEST  ANALYSIS 

The  overall  objective  of  this  task  was  to  provide  support  for  the 
hot-fire  testing  of  the  dual  throat  hardware  by  NASA/MSFC.  ALRC's  first 
objective  was  to  review  the  MSFC's  test  plan  and  provide  whatever  consultation 
might  be  necessary  during  testing.  Secondly,  ALRC  would  review  the  hot-fire 
data  and  interpret  test  results  by  using  the  aerodynamic  model  and  previous 
cold-flow  data.  Test  data,  would  dso  be  used  to  correlate  dual  throat  nozzle 
Mode  II  performance.  Thirdly,  corrections  to  the  aerodynamic  model  believed 
advisable  were  to  be  incorporated  into  the  model,  and  further  model  improve¬ 
ments  were  to  be  recommended. 

The  hardware  tested  was  the  same  as  that  used  during  cold-flow 
testing  (NAS  8-32666),  (Ref.  1).  For  these  tests,  the  small  primary  nozzle 
(ep  =  1.72)  and  the  large  secondary  nozzle  (e$  =  3.0)  were  used.  A  layout 
of  the  test  hardware  is  shown  in  Figure  VI-1.  The  proposed  test  plan 
called  for  testing  to  be  conducted  at  three  different  nozzle  spacings,  using 
GH2  as  the  bleed  gas.  Recommendations  were  made  by  ALRC  for  testing  to  be 
done  with  various  bleed  flow  gas  compositions  and  temperatures  to  verify  the 
improvements  made  to  the  aerodynamic  bleed  flow  model  to  simulate  such  condi¬ 
tions.  Funding  and  time  limitations  did  not  permit  testing  a  survey  of  bleed 
flow  gases,  although  one  series  of  tests,  at  the  larger  nozzle  spacing,  was  per¬ 
formed  with  GN^  as  the  bleed  gas.  The  tests  conducted  are  summarized  in  Table  VI-I. 

The  improved  aerodynamic  bleed  flow  model  was  used  to  define  the 
"blowoff"  condition  to  provide  a  safe  first  test  operation.  The  predicted 
trend  in  bleed  flowrate  versus  nozzle  spacing  is  shown  in  Figure  VI-2,  for 
both  a  GH2  and  a  GN2  bleed.  A  significant  influence  of  bleed  flow 
properties  on  operating  characteristics  was  predicted.  These  predictions 
indicated  that  one  and  one-half  to  two  times  the  bleed  flow  would  be 
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Figure  VI-2.  Predicted  Trend  in  Bleed  Flowrates  versus  Nozzle  Spacing 
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required  when  using  GN2  as  opposed  to  using  GH2.  This  difference  is  due 
to  the  influence  that  GN2  has  on  the  shear  layer  mixing  equations  with  its 
larger  molecular  weight  and  smaller  specific  heat  (Cp).  One  item  of  concern 
was  the  large  bleed  flowrates  predicted  using  GN2  and  for  the  large  nozzle 
spacing  using  GH2.  The  assumptions  used  in  developing  the  bleed  flow  model 
were  expected  to  limit  the  model's  validity  and  applicability  to  small  bleed 
flows  of  approximately  10%  or  less.  As  noted  in  Section  III,  significant 
reductions  in  bleed  flow  from  the  blowoff  condition  can  be  made  without 
measurably  inducing  a  loss  in  total  pressure. 

B.  TEST  RESULTS 

1.  Data  Evaluation 


Wall  pressure  profiles  for  the  various  test  configurations 
are  shown  in  Figures  VI-3  through  VI-5  for  various  nozzle  spacings  and  in 
Figure  VI-6  for  N2  bleed  flow.  In  general,  the  pressure  data  are  reasonable 
and  similar  to  the  previous  cold-flow  test  results  (Ref.  1  ).  A  comparison  of 
the  hot-fire  and  cold-flow  pressure  profiles  is  shown  in  Figure  VI-7  for  the 
large  nozzle  spacing  (Le/Rtp  »  5.06)  with  GN2  bleed.  The  plume 
impingement  point  is  shown  to  move  toward  the  secondary  throat  with  increasing 
bleed  flow,  as  was  seen  in  cold-flow  testing.  Also,  the  plume  impingement 
shock  strength  is  decreased  with  bleed  flow,  which  is  consistert  with  the 
cold-flow  data  and  the  aerodynamic  bleed  flow  model  predictions.  Test  results 
also  show  that  increased  nozzle  spacing  requires  increased  bleed  flow  to 
attach  the  plume  to  the  secondary  throat  which  is  also  consistent  with  both 
the  cold-flow  data  and  the  model  predictions. 

Although  the  hot-fire  data,  overall,  were  consistent  with  the 
cold-flow  data,  some  anomalies  did  exist.  All  the  test  configurations  had 
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some  questionable  pressure  measurements  at  the  secondary  throat.  The  wall 
pressure  would  increase  due  to  the  plume  impingement,  decrease,  and  then 
increase  again  at  the  secondary  throat  before  decreasing  as  a  result  of  being 
expanded  through  the  secondary  nozzle.  Cold-flow  data  show  the  secondary  wall 
pressure  to  increase  due  to  the  plume  impingement  and  then  continually 
decrease  as  the  gas  exits  through  the  secondary  nozle. 

Another  anomaly  is  seen  in  Tests  15,  18,  23,  and  34,  shown  in 
Figures  VI-6,  VI-5,  and  VI-4,  respectively.  The  pressure  in  the  secondary 
chamber  increases  and  then  decreases  prior  to  plume  impingement.  This  trend 
was  not  seen  during  cold-flow  testing  and  is  probably  the  result  of  small 
pressure  measurement  errors. 

Test  8,  shown  in  Figure  VI-3,  is  also  inconsistent  with  all 
other  test  results.  Its  measured  pressure  profile  is  consistent  with  a  mini¬ 
mum  bleed  flow,  yet  the  measured  bleed  flow  was  the  largest  with  that  con- 
figuation.  This  may  be  due  to  a  leak  in  the  bleed  circuit  after  the  flowrate 
measurement  or  an  erroneous  flowrate  measurement. 

2.  Model  Comparison 

Although  some  uncertainties  do  exist  with  the  test  data,  the 
experimental  results  have  verified  the  bleed  flow  model's  predictions.  An 
example  of  this  is  shown  in  Figure  VI-8.  The  bleed  flow  model  predicts  nearly 
exactly  the  recirculation  pressure  for  a  given  bleed  flow.  Also,  the  model's 
prediction  that  more  GN2  bleed  would  be  needed  to  achieve  the  same  recircu¬ 
lation  pressure  as  with  GH2  was  substantiated  with  the  hot-fire  test 
results. 


The  hot-fire  results  also  verified  the  correlation  of  the 
plume  scaling  factor  (FCTR)  with  nozzle  spacing.  Figure  VI-9  shows  that. 
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Figure  VI-8.  Experimental  Results  Have  Verified  Aerodynamic  Model  Predictions 
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for  zero  bleed  flow,  the  bleed  flow  model  predicts  exactly  the  recirculation 
pressure.  Some  inaccuracy  does  exist  between  the  data  and  predictions  at  the 
larger  bleed  flows  for  small  nozzle  spacing.  This  may  have  resulted  due  to 
inaccuracies  in  the  plume  correlation.  Another  important  point  verified  is 
that  a  noticeable  decrease  in  bleed  flow  from  that  required  for  blowoff  may 
result  in  a  very  small  stagnation  pressure  loss  from  plume  impingement. 

3.  Engine  Performance 

Two  means  of  measurement  were  used  to  determine  engine  per¬ 
formance,  characteristic  velocity  (C*),  and  specific  impulse  (Isp).  Table 
V I - 1 1  presents  the  method  used  to  determine  these  parameters.  The  format  for 
determining  specific  impulse  is  the  same  as  that  used  during  previous  cold- 
fl ow  testi ng  (Ref.  1 ). 

Tests  1  through  7  used  only  the  primary  nozzle  to  check  out 
the  facility  and  determine  combustion  performance.  The  resulting  C*  effi¬ 
ciency  proved  to  be  very  low  (70-75%)  probably  due  to  the  addition  of  film 
cooling  holes  in  the  outer  ring  of  the  primary  injector.  Testing  with  the 
dual  throat  configuration  hao  similar  low  primary  combustion  performance, 
resulting  in  a  correspondingly  low  Isp  efficiency.  As  a  result,  this 
unusually  large  combustion  inefficiency  tended  to  mask  the  relatively  smaller 
inefficiencies  attributable  to  the  dual  throat  configuration  (i.e.,  boundary 
layer  and  divergence  losses). 

Nevertheless,  the  performance  determined  from  the  hot-fire 
test  was  correlated  using  the  simplified  prediction  methodology.  Estimates  of  the 
divergence,  boundary  layer,  and  kinetic  efficiencies  were  used  to  define  a  perfect 
injector  performance  for  the  test  configuration.  This  value  was  compared  to  the 
measured  Isp  for  a  typical  test  to  "back  out"  an  energy  release  efficiency  as  shown  on 
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Table  VI-III.  This  method  resulted  in  an  indicated  energy  release  efficiency 
of  76%,  which  is  of  the  same  magnitude  as  the  characteristic  velocity  effi¬ 
ciency  (73%).  Therefore,  while  the  hot-fire  results  cannot  be  used  to  provide 
Vernier  calibration  of  the  dual  throat  performance  model,  it  appears  that  this 
measured  result  is  consistent  with  the  model  and  that  the  unusually  low  spe¬ 
cific  impulse  is  primarily  the  result  of  a  low  energy  release  (combustion) 
efficiency. 


PERFORMANCE  COMPARISON  OF  HOT-FIRE  TESTS/AERODYNAMIC  MODEL 
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PROGRAM  INPUT,  AERODYNAMIC  MODEL 


CARD  1: 


The  1st  card  to  be  input  is  used  to  provide  a  printed  heading  for 
the  computer  output.  Columns  2  throuqh  72  are  available  for  text.  This 
card  must  always  precede  each  case.  The  remainder  of  the  input  data  is 
read  using  NAMELIST  with  the  name  $DATA,  as  described  below. 


SPAT  A 

Geometry;  Option  1,  Cylindrical  Secondary  Wall,  see  Figure  1. 


INPUT  ITEM 

I  WALL.  =  1 , 

RS 

RSP 

THETAP 

EPSP 


DEFINITION 


The  secondary  nozzle  wall  is  a  cylinder 
r*,  radius  of  the  secondary  nozzle  wall 
r*p,  radius  of  the  primary  nozzle 
6p,  exit  cone  half  angle,  primary  nozzle 
Cp ,  expansion  ratio  of  the  primary  nozzle 


The  items  r*  and  r*p,  above,  can  be  input  with  units  of  inches,  feet,  cm, 
etc.  The  input  item  cp,  above,  is  used  only  to  determine  the  Mach  number 
at  the  exit  of  the  primary  nozzle  (not  required  if  MEPI  is  input). 


Figure  i  .  Nozzle  Geometry  for  Option  1 
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Geometry;  Option  2,  Full  Geometry  Option,  see  Figure  ?■ 


I WALL  =  2, 

Complete  wall  geometry  specified  for  both 
primary  and  secondary  nozzles 

Primary  Nozzle 

ECP 

>.c  ,  primary  nozzle  chamber  contraction  ratio 

RIP 

K 

Rip.  wall  radius  connecting  chamber  and  inlet 
of  primary  nozzle,  divided  by  r* 

•Jjp,  inlet  angle  for  primary  nozzle 

THIP 

RWTUP 

Ru  ,  wall  radius  on  upstream  side  of  the 
throat,  primary  nozzle  ,  divided  by  r* 

RSP 

P 

r*D,  throat  radius  of  primary  nozzle  ,  inches 

RWTDP 

^dD*  wall  radius  on  downstream  side  of  the 
thfoat,  primary  nozzle  ,  divided  hy  r* 

n 

THETAP 

K 

!>p,  exit  cone  half  angle,  primary  nozzle 

EPSP 

c  ,  expansion  ratio  of  the  primary  nozzle 

EC 

RI 


THETA1 

RWTU 


RS 

RWTD 


THETA 

EPS 


Secondary  Nozzle 

cc,  secondary  nozzle  contraction  ratio 


Rj ,  wall  radius  connecting  chamber  and 
inlet  of  primary  nozzle  ,  divided  by  r* 


e°j,  inlet  angle  for  secondary  nozzle 

Ry,  wall  radius  on  the  upstream  side  of  the 
throat,  secondary  nozzle  ,  divided  by  r* 

r*,  throat  radius  of  secondary  nozzle,  inches  or  e 

Rrt,  wall  radius  on  downstream  side  of  the 
throat,  secondary  nozzle  ,  divided  by  r* 

e°,  exit  cone  half  angle,  secondary  nozzle 

e,  expansion  ratio  of  the  secondary  nozzle 


Geometry;  Options  3,  4  or  5.  Full  Geometry  Option  with  Spline  Fit 
Exhaust  Contours 


For  the  spline  fit  wall  contour  options,  contour  points  are  input  to 
replace  the  conical  exit  cone  of  either  the  primary  nozzle,  secondary  nozzle, 
or  both.  Other  inputs  are  the  same  as  in  geometry  Option  2,  described  above, 
except  that  the  inputs  for  ep,  and/or  e  are  not  required.  The  exit  cone 
half  angle  is  interpreted  as  the  spline  contour  attachment  angle.  The  spline 
options  are  as  follows: 


1WALL  =  3, 
I WALL  «  4, 
I WALL  =  5, 


Primary  Wall 

Spline 

Cone 

Spline 


The  spline  coordinates  are  input  as  follows: 


Secondary  Wall 

Cone 

Spline 

Spline 


X4P(2)* 

Y4P(2) 

N4P 

TH4P 

X4(2) 

Y4(2) 

N4 

TH4 


Xpj,  axial  coordinates  for  primary  wall 
spline  contour 

yp j ,  radial  coordinates  for  secondary  wall 
spline  contour 

np,  1  ■  1,  2,  ...  np  above.  np  <_  20 

8°ep»  spline  contour  exit  angle  for  pri¬ 
mary  nozzle,  degrees 

xi,  axial  coordinates  for  secondary 
wall  spline  contour 

y^,  radial  coordinates  for  primary  wall 
spline  contour 

n,  1 ,  ■  1 ,  2  . . .  n  above.  N  <_  20. 

e#e,  spline  contour  exit  angle  for  secondary 
nozzle,  degrees 


♦The' first  point  on  a’  spline  contour  is  automatically  calculated  by  the 
program  from  the  exit  contour  half  angle,  6p°,  or  e°,  thus  each  of  these 
input  vectors  begins  with  item  2. 

All  \,y  coord  i  rid  It-  v.ilucs  are  normal i/cd  (i.e.  divided)  by  the 
corresponding  norz’e  threat  radius. 
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Re  la  t  i  vePos  itioning  of  Pri  inary  and  Secondary 
NozzTes~TlWALL  >  2  Cases  Only) 


Lp/r*p,  axial  distance  from  the  exit  plane 
of  the  primary  nozzle  to  the  throat  plane 
the  secondary  nozzle,  divided  by 
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FLOW  PROPERTIES  AND  OPTIONS 

ASSUMED 

INPUT  ITEM  DEFINITION  VALUE 


FCTR  Factor  for  scaling  plume  shape.  Used  as  a  1. 

multiplier  for  the  parameter  M./y  in  the  Herron 
method. 


IPNCH  If  IPNCH  =  1,  geometry  inputs  for  TDK  will  be  0 

punched  in  the  NAMELIST  format  . 

If  IPNCH  -  0,  no  punched  output,  only  printed 
output. 

MEPI  M  .  (REAL  type)  exit  Mach  No.,  primary  nozzle.  0 

P 

If  Mr  <1  then  values  will  be  calculated  from  e 
EP  P 

assuming  one-dimensional  flow 

NW  nw,  number  of  entries  in  WSWP  13 

n.j.50 


PBPO  Initial  guess  for  the  pressure  ratio  for  the  plume 

boundary.  .  If  not  input,  the  program  will  calcu¬ 
late  a  value  corresponding  to  the  area  of  the 
secondary  nozzle  throat  divided  by  the  area  of  the 
primary  nozzle  throat  using  the  relationship  for 
one-dimensional  flow. 


WLPRN  If  a  value  other  than  zero  is  Input,  print  of  the 

secondary  and  primary  nozzle  wall  tables  will  be 
suppressed  (to  be  used  on  successive  cases  when 
the  wall  geometry  does  not  change.) 

WSWP  (1)  (Ws/Wp)i  *  100,  percent  bl  .ed  flow  given  for 

these  values,  i  =  1,  — ,  n  for  which  base 
pressure  and  recirculation  temperature  and 
primary  mass  fraction  in  the  recirculation 
region  ,  are  solved. 


0. 


0,1 ,2, 3, 4, 5, 

6.7,8.9.10, 

12,15 
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INPUT  ITEM 


DEFINITION 


SUGGESTED 

VALUE 


TRTO 

VTo 

Recirculation  temperature  ratio  initial  guess 

NOTE:  TRTO  is  between  TOBASE  &  1.0 

1.0 

TOBASE 

Tot/To 

Stagnation  temperature  ratio  of  bleed  flow  (secondary 
flow)  to  primary  flow 

DELV 

Perturbation  fraction  for  Newton  solution  procedure 

.0001 

DELVMX 

Maximum  fractional  change  of  P. /P  and  T  /T  during 
solution  00  r  0 

.3 

EPSGGT 

Convergence  criteria  for  bleed  flow  Ws/Wp 

.0004 

EPSTOB 

Convergence  criteria  for  To^/Tq 

.0004 

ITERPR 

Iteration  print  option,  set  equal  to  1  for  print 

0 

ITERMX 

Max  #  of  iterations 

20 

CKN 

Nash  Factor  N 

.4 

LPBPO 

Logical 

Variable 

Find  Flow-off  Base  Pressure  (PBPOBF)  if  true 

True 

LPL0T 

Logical 

Variable 

Draw  Plots  if  true 

False 

XMINPL 

Minimum  axial  coordinate  for  plot  axis 

XMAXPL 

Maximum  axial  corrdinate  for  plot  axis 

XDISPL 

Length  of  axis  (for  axial  coordinate) 

YMINPL 

Minimum  radial  coordinate  for  plot  axis 

YMAXPL 

Maximum  radial  coordinate  for  plot  axis 

YDISPL 

Length  of  axis  (for  radial  coordinate) 

PBPOBF 

Blow-off  Base  pressure  to  provide  a  limit  during 
iteration  or  initial  guess  for  finding  blow-off 
pressure  if  LPBPO  =  .TRUE. 

.01 

(Note:  This  is  then  stored  in  PBPO). 


INPUT  ITEM 
66TB  F 

CPF 

CP0 

MWF 

MW0 

ITYFR 

YFRX 


DEFINITION 


SU66ESTED 

VALUE 


Blow-off  Secondary  Mass  flow  (percent)  to  provide  a 
limit  during  iteration.  (This  is  automatically  set 
if  LPBPO  =  .TRUE.) 

Heat  capacity  (constant  pressure)  of  primary  mass  flow 

Heat  capacity  (constant  pressure)  of  secondary  mass  flow 

Molecular  weight  primary  mass  flow 

Molecular  weight  secondary  mass  flow 

#  of  iterations  to  be  taken  on  mass  flow  fraction 
during  each  base  pass. 

Primary  mass  flow  fraction  in  recirculation  region 
first  estimate  for  YFR. 


2 

.5 
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EXTENDED  NOZZLE  WALL  OPTION 

Given  an  input,  nozzle  wall  contour,  this  option  will  calculate  an 
extended  nozzle  wall  contour  that  gives  maximum  thrust.  This  is 
accomplished  by  first  calculating  ttie  flow  field  that  is  dependent  on  the 
given  input  nozzle  wall  contour,  i.e.  the  downstream  boundary  of  the  flow 
field  is  the  right  running  characteristic  extending  from  the  nozzle  lip 
to  the  flow  axis.  Shocks  are  ignored.  Next,  a  circular  arc  of  radius  R^ 
is  attached  tangent  to  the  nozzle  lip.  The  optimization  options  existing 
in  the  program  are  then  used  to  calculate  the  maximum  thrust  wall  contour. 
The  following  additional  input  items  are  required. 

IWL0PT  =  1  for  the  extended  wall  option.  The  program  assumes 

IWhtil'T  -  0. 

RD  =  wall  radius  ratio  (R/R,  )  of  the  circular  arc,  Rd,  attached 

tangent  to  the  nozzle  lip. 

NOZZLE  GEOMETRY 

The  items  RWTU,  RWTD,  and  RSTAR  are  as  previously  defined,  but  apply 

to  the  input  nozzle  wall  contour.  See  '''igure:  la.b.c.  and  d  for  wall  contour 
wpt ions. 

IWALL  =  wall  option  flag .  The  value 

assigned  this  item  determines 
the  type  of  wall  to  be  specified. 

0  =>  arbitrary  contour 

1  =>  cone;  9,  e 


2.  =>  parabolic  contour 

3  =>  circular  ore  contour 

4  =>  spline  fit  contour 

5  =>  cone;  r  ,  x 

n  n 

If  IWALL  =  0  (arbitrary  contour) 

PW(1)  =  The  wall  coordinate  points 

,  Xj ,  are  to  be  input 

sequentially: 


PW(1)  =  rj,  Xj.  r2>  x2,  .  .  . 

*  *  "  rn '  xn '  ®  ®  *  * 
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'  a  pair  of  zero  coordinates 
marks  the  end  of  the  table. 

A  value  of  n  «  300  is  required. 

This  table  is  normalized  to 
the  throat  radius  and  the 
origin  is  at  the  throat  so  that 
necessarily  r^  =  1,  =  0. 

If  1WAL?  ■■■;  1  (cone;  6,  () 

RWTD  =  as  above. 

THJI  =  exit  cone  half  angle  (decrees). 

EPS  =  expansion  ratio 

2  or  3  (parabolic  contour  or  circular 

~ ~  arc  contour) 

=  as  above 

=  wall  inflection  angle  (degrees) 

RWMAX  =  nozzle  end  point  radial 

coordinate,  r 

n 

ZWMAX  =  nozzle  end  point  axial 

coordinate,  xn 

If  IWALL  =  4  (spline  fit  contour) 

As  in  IWALL  options  1,2,  and  3  the  program  assumes  the  throat  to  begin 
with  a  circular  arc  downstream  of  the  minimum  point.  The  final  point  on  this  cir¬ 
cular  arc  is  located  at  an  Inclination  of  THJI  degrees.  The  contour  exit  angle,  THE 
degrees,  must  also  be  input.  The  contour  to  be  spline  fit  between  THJI  and  THE 
must  be  input  in  tabular  form.  The  first  point  on  the  spline  fit  is  automatically 
located  by  the  program  on  the  circular  arc  at  the  point  of  inclination  THJI. 


RWTD  =  As  above 

THJI  =  RWTD  and  the  angle  THJI  define 

the  end  point  of  the  circular  arc 
and  the  beginning  of  the  spline  fit 
(degrees ) 

NRZS  =  n,  the  number  of  points  to  be  spline 

fit.  n  «  l|0  is  required. 

PWRS(2)  =  radial  wall  coordinated  to  be  spline 

fit. 

PWRS(2)  -  r2,  r3 . rn 


If  IWALL  » 

RWTD 

THJI 
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PWZS(2) 


axial  wall  coordinates  to  be  spline 
fit. 

I’WZS(2)  x„ ,  x  ,  .  .  x 
A  o  n 

THE  =  nozzle  exit  angle 

If  IWAL1.  =  5  (cone:  r  .  x  ) 

- '  n  n 

RWTD  =  As  above 

RWMAX  =  As  above 

ZWMAX  =  As  above 


APPENDIX  B 

DERIVATION  OF  DUAL  EXPANDER  EQUATIONS 


DERIVATION  OF  AREA  EQUATION 


A  =  2it  6  R 


=  Outer  Radius 
=  Inner  Radius 
=  Radius  @  Mid  Point 


=  Angle  Gap  Inclined 


5-  bh  -  \ 

COS  0 


(RH  "  Rl_^sec  0 


The  projected  area  of  the  annulus  can  be  found  by 
Ap  =  TT  (Rh2  -  Rl2) 


The  annulus  area  is  then 


A  =  A. 


COS  0 


A  ■  ■  (rh2  •  rl2> 

COS  0 


*  (Rh  -  Rl)  (Rh  +  Rl) 


cos  0 


*  RH  +  RL 


Rh-  Rl  =  6  cos  0 


Substitute: 


*(icos  0)  2  R„ 


cos  0 


2  if  6  Ru 


t 

I 


I 


DERIVATION  OF  INCLINATION  ANGLE  EQUATION 

Cos  0=2  [Rj-  '  re  '  °L  cos  °L] 
c  P 


From  Figure  B-l 


Ap  =  — —  [A  R2  -  A  R  R  -  A  R  D.  COS  6,  ] 

bc  cos  os  C  C  bp  fcc  L  L 


A  TT 

=  - - -  [Rr  -  Rp  -  D.  cos  0.  ] 

cos  os  C  L  L 

cos  0$  =  4  "  REc  [Re  -  Re  ■  °l  cos  0L^ 
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Figure  B-l 
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Substitute: 


2  IT  R^ 

c 


cos  0 


Dl  cos  eL) 


DERIVATION  OF  THROAT  RADIUS  EQUATION 


2  sec  0$  /I  -  sin 


sin  Q  cos  c  ,  . 

S  C  '  -j 


-) 


-  R* 


2  sec  G  / Rf  +  D.  cos  0.  -  Rr  sin  e 
S  LP  L  L  c 


sin  es  cos  - 


0 


Start  with  the  relationship  from  Fiqure  B-l 


rt  '  rt 

_H _ _L 

COS  O, 


and  r.  =  2  R.  -  r. 

n  Ts  \ 


2  (Rt  -  rt  : 
S  L 


cos  n 


From  Fiqure  B-2 


tL  el 


v -  =  cos 

X  c 


X  sin 


=> 


rt  =  rE 


c  sin  o. 


cos  0 
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From  Figure  B-l 
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INPUT  DESCRIPTION 
DUAL  EXPANDER  SHEAR  LAYER  MODEL 


The  purpose  of  this  computer  program  is  to  calculate  the 
base  pressure  that  will  occur  when  a  uniform  supersonic 
flow  expands  around  the  corner  of  a  bluff  body.  An 
example  is  shown  in  Figure  1,  which  corresponds  to  the 
Dual  Expander  engine  during  Mode  II  operation. 


The  method  of  Korst  is  used  to  predict  the  base  pressure. 

The  method  applies  only  to  flows  that  are  not  fully  separated 
The  shear  layer  must  approach  the  flow  axis  and  be  turned  at 
the  so  called  "neck  region".  The  program  allows  a  bleed 
flow  to  be  specified,  provided  this  secondary  flow  is  not 
so  great  that  the  flow  becomes  fully  separated.  The  pro¬ 
gram  allows  the  molecular  weight,  specific  heat,  and  total 
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enthalpy  of  the  secondary  stream  (bleed)  to  be  different 
than  that  of  the  primary  stream.  The  method  of  using  this 
program  is  described  below. 

Solutions  can  be  obtained  for  the  following  situations: 


I) 

zero  bl eed 

flow, 

2) 

bleed  flow 
properties 

with  properties 
of  the  primary 

identical 
stream , 

to  the 

3) 

bleed  f 1 ow 
properties 

with  properti es 
of  the  primary 

different 
s  tream . 

than 

the 

Case  3 ) , 

a bove ,  i s 

the  general  case 

solved  by 

the 

program 

and  requires  the  determination  of  three  unknown  variables 

(P./P  ,  H  /H  , Y  )  of  three  unknown  variables  such  that 

b  op  r  op  pr 

three  functional  constraints  are  satisfied.  Cases  1)  and 

2),  above,  are  special  cases  of  case  3)  and  require  the 

determination  of  only  one  variable  (P./P  )  such  that  one 

b  op 

function  is  satisfied.  In  all  three  cases,  initial  esti¬ 
mates  of  the  unknown  variables  must  be  supplied.  These 
estimated  values  must  lie  within  the  radius  of  conver¬ 
gence  Newton's  method.  Estimates  can  be  obtained  from 
solutions  previously  found.  Each  run  on  the  computer  will 
use  its  final  values  as  the  first  estimate  for  the  suc¬ 
ceeding  case.  This  allows  parametric  studies  to  be  run 
efficiently  by  the  systematically  changing  of  input. 

The  program  also  uses  bounds  for  the  unknown  variables 

so  as  to  enhance  convergence  of  the  numerical  method. 

Variables  H  /H  „  and  Y  are  automatically  bounded  by  the 
r  op  pr 

program  and  have  caused  no  difficulties  thus  far.  The 

variable  P./P  ,  however,  requires  carefull  consideration, 
b  op  _5 

A  satisfactory  lower  bound,  10  ,  is  assumed;  but  the 

upper  bound,  {P^/PQ^)  max,  must  be  input  by  the  uses  as 
described  in  the  following  steps. 
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Step  1)  Using  y  and  the  approach  flow  Mach  No.,  M  , 
look  uppthe  Prandtl-  Meyer  angle,  v  ,  in  tRe 
Gas  Ta bl es  .  p 


Step  2)  Calculate  the  Pra rdtl -Meyer  angle,  v,  that 
will  turn  the  flow  parallel  to  the  axis  by 
adding  the  approach  flow  angle  to  i.  e., 

v*  v  +  0 . 

P  b 


Step  3)  Look  up  the  flow  pressure  ratio  for  parallel 
flow,  P/P  ,  corresponding  to  v.  (P./P  )  max 

must  be  l$ss  than  this  value,  i.e.,  op 

<W  ""  ■  -95  (P'Po> 

A  good  estimate  for  the  base  pressure  ratio  is 


p./p  .. 


.5  (P/P  ) 


One  card  is  read  and  printed  to  label  the  program  output. 
This  card  is  followed  by  the  N  A  M  r  L  I  S  T  input  data  set  des¬ 
cribed  below,  followed  by  a  listing  of  cards  for  a  sample 
case. 


NAMELIST 

NAME 

EXPLANATION  AND  COMMENTS  | 

UNITS 

• 

OE  FAULT 
VALUE 

SCUALEX 

1 

cps  =  ! 

C  ,  Specific  heat  at  constant  pressure 

ca  1  / gm  °K 

i 

cpp  =  ! 

Ps  for  the  secondary  (bleed)  flow. 

C  ,  Specific  heat  at  constant  pressure 

t 

1 

c  a  1  /  gm  0  k! 

i 

AMWS  =  ! 

Pp  for  the  primary  flow. 

M  ,  Molecular  weight  for  the  secondary 

1 

1 

i 

■ 

AMWP  = 

ws  (bleed)  flow. 

M  Molecular  weight  for  the  primary 

w  p’  flow. 

H  ,  Total  enthalpy  for  the  secondary 

i 

1 

H0S  = 

1 

ca  1  /gm  1 

I 

H0P  = 

°s  (bleed)  flow.  H  =  C  T 

°s  ps  s 

H  ,  Total  enthalpy  for  the  primary 

ca 1 /gm 

i 

P0P  = 

!  °P  flow.  H  =  C  T 

°P  V  P 

P  ,  Total  pressure  for  the  primary 

ps  i  a 

WSWP  = 

0  p  flow. 

W  /W  ,  Ratio  of  the  mass  flow  rate  of 
p  the  secondary  (bleed)  flow  to 
the  mass  flow  rate  of  the  pri- 
ma  ry  flow 

i  r,  Radius  of  base 

none 

SRB  = 

inches 

RWRB  = 

r  /r.  ,  Wake  radius  ratio 

(or  cm )  | 

none  i 

.  5 

AREA  = 

w  b 

The  throat  area  of  the  primary  stream 

.  2 
,  1n 

|  1  /  n 

THETAB  = 

0 k ,  Approach  flow  angle,  i.e.  lip  angle 

(or  cm  j  1 

i 

i 

degrees 

1 

i 

I 

AMP  * 

of  the  lower  wall  of  the  annular 
nozzle. 

Mp,  Mach  No.  of  the  primary  flow 

none 

f 

1 

i 
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NAMELIST 

NAME 


EXPLANATION  AND  COMMENTS 


UNITS 


PBP0P1  = 

Pb/P  ,  1st  estimate  for  the  ratio  of 

p  base  pressure  to  the  total  pres¬ 
sure  of  the  primary  flow. 

none 

PMAX  = 

{Pb/PQ  )max,  upper  bound  for  P ^ / P q  . 

^  Use  procedure  previously  ^ 
described  to  calculate  this 
value. 

none 

HRH0P1 

Hr/HQ  ,  1st  estimate  for  the  ratio  of 
p  the  enthalpy  of  the  gas  in  the 
base  recirculation  region  to 
the  total  enthalpy  of  the  pri¬ 
mary  flow. 

none 

YPR1  = 

Y  ,  Mass  fraction  of  bleed  flow  gas 
r  contained  in  the  base  recircula¬ 
tion  region  . 

none 

AN  = 

N,  Nash  factor 

none 

GC0N  = 

Convergence  criteria  for  yd 

none 

EJCNV  = 

Convergence  criteria  for  n • 

none 

WF(1)  = 

Wi ,  weighting  factors  used  in  testing 
the  convergence  of  the  varables  sought 
by  Newton's  method. 

none 

EPSI  = 

Convergence  criteria  for  Newton's 
method.  Convergence  occurs  when 

.S  (WIFI ) 2<  EPSI 

l  *1 

NUN  = 

n,  number  of  varables  to  be  determined. 

If  n  =  3,  then  P./P  ,  H  /H  ,  and  Y 
bo  r  o  p 

P  p  yr 

will  be  determined. 

If  n=l ,  then  Pb/PQ  will  be  determined. 

P 

Always  set  H/H  =  1 ,  Y  =1 
r  o  p 

P  Kr 

PERTV  = 

Perterbatlon  factor  for  calculating  par¬ 
tial  derivitlves  for  Newton's  method 
(see  write-up  of  subroutine  NEWT) 
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PERTV 


NAM  I  l  I 
NAME 

SK  = 


I  DBUG 

MAX  IT 

MAXITP 

MAXITE 

LAST  = 

SEND 


S  T 


EXPLANATION  AND  COMMENTS 


UN  ITS 


DEFAULT 
_V  AJ.  UE_  _ 


proportionality  factor  relating  shear 
layer  velocity  thickness,  6,  to  the 
wake  radius  ratio,  r  /r  .  The  pro¬ 
gram  has  been  calibrated  versus  data 
for  air  without  bleed  flow  and  it  was 
found  that  k  =  4,  i.e.  it  is  assumed 
that 


r  /  r 
w  b 


=  k 


*/rb)/a 


where 


x/r  is  the  non-dimensional 
b  shear  layer  length 

n  i  s  t  h  e 
meter . 


jet  spreadi ng  para- 


a 


12  +  2  .  7  58M, 


none 


4  . 


jFlag  for  debug  ing  print-out 
If  IDBUG  =  0,  no  debug  print 
|  If  IDBUG  =  1,  print  debug  output 


0 


Maximum  number  of  iterations  allowed  for 
Newton's  method. 


1  0 


Maximum  number  of  iterations  allowed  in 
calculating  y^ 

Maximum  number  of  iterations  allowed  in 

calculating  n  . 

J 

Flag  to  avoid  error  message  that  results 
from  reading  end  of  file  for  input.  For 
the  last  case  to  be  run,  set 

LAST  =  1 


25 

25 

0 


I 
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TABLE  1,  CARD  LISTING  FOR  A  SAMPLE  CASE 


**  CALIPWATIUN  CASfc  **  bAPMAxl.a 

SIHiALEX 

CPS=.2396HQ7,CPP=.^49bB97»  AMwSs,?<v 
HilS=  1  .  ,  H(jP=  1  .  ,  Pl)Ps  l  .  , 


SUfls  1  .  , 

1  HE  T AHsO  .  / 

AMP  =  <?.  3, 

AG  =  1  .  , 
w9wP=o . , 

PM/.x  =  .  079,PBP(jP1  =  .0u, 
NUNS  1  , 

ANs.a, 

LA5T  =  1 , 

SEND 


A  I  W 


A  I- '  w  P  s  e  9  .  , 
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I  !NAi  n  im  ■!  i ■  i ■  1  - 1  i 1 1  ;..-j 


DUAL  NO/ /I  I  A|  i*  WP  Mi..';  ,  .V,:.|  .  .p, 

CONI  r-v  i  fs/i:,  ft  -  j  * 


NAML  corns 

National  Aeronautics  ami  Sp.icr  Ad  iim strut  ion 

Georye  C.  Marshall  spa.*-  M  !■:«■«  ( .  n  t  r. . 

Marshall  Spare  fhyht  f* .liter Ai 


Attn:  Administrative  Service-.  /  AS24D  3 

Technoloqy  Util  ij.il  ir.n  Of  t  b  »•  /  ATO*  1 
Associate  Director  for  Manayer./nt  /  l 

FM12-34 

J.  L.  Sanders  /  PD 1 1  1 

F.  W.  Braam  /  LP24  10 


Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  Ohio  44 135 

Attn:  Library  1 

R.  Priem  /  501-6  1 


NASA  Headquarters 
Washington,  L-.C.  20516 

Attn:  F.  W.  Stephenson  /  RP  1 

Goddard  Space  t  >r*.--r 

Greenbelt,  Miry!. mi  Z0 771 

Attn:  Library  1 

John  F.  Kerr, rt.,. j 
Cocoa  Beach ,  f  1  ••  i  •  •  , 

Attn  l  i lira*  ,  ] 

Lyndon  B.  Johns,  n  '  race  Cin’n 
Houston,  Te<-*s  /?{J6H 

Attn:  library  1 

Lar.glc-y  Research  Center 
langely  Station 
Hampton,  Virginia  23365 

Attn:  l ibrary  1 


Jet  Propulsion  laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91101 

Attn:  Library  i 

Defense  Documentation  Center 
Cameron  Station 
Building  5 
5010  Duke  Street 
Alexandria,  VA  22314 

Attn:  TISIA  2 

NASA  Scinnt.  ifi.  Technical  Infru  :i  ion  r<i<ility 
P.0.  Box  H/‘>/ 

Bal  limn re  Wash  i»*-il on  tntorn.i t  j«n  a  I  i  >  pot  t 
Boltiinoti-,  Kb  2  I, '40 

Attn:  Accessioning  D*»pt.  6 

U.S.  Army  Missile  Co*  -and 

Redstone  Scientific  Information  Center- 

Reds  tone  Arsen.*!  ,  AL  358*1. . 

Attn:  Docur.  t  Selection  3 

Air  Forte  Pocl'i’t  prnp.jl  ,  <(.ii  !  -b-.r  it  ;**y 
Fdw-’ rdr. ,  CA  i  i ‘ '  i 

Af  to  ■  libr.-ry  1 

At  no  Id  I  :  n>  •  1 11  ;  *( f  ' 

Air  fr-t  «  e  ,v  i 

1  o  1 1  a  hi  .ii  i ,  j  r-i  \  j 

At  In  L  itn  ,*  /  1 


ful'  t  <•  !  1  Vtll  ■ 

[‘ •  v  i  -n  ,i*  1 1  ! 

6631  •. 

C •* n-u;  .  (Ctr  l.  ,  l,\  ‘i  1  I 

At  f.  r  i :  Library 

Pratt  and  Whitney  Air; 
Govern-ier.t  ProiR,  ts  D» 
L'n  i  t'.'.l  let  hr  . ;  oow  s  C  . 
West  Palm  Ib-vh.  II  3 

Attn:  Librae, 

Bell  Aerosyst ens ,  Inr. 
Box  1 

Buffalo,  NY  14240 

Attn:  Library 

Boeing  Company 
Spate  Division 
P.0.  Bo a  868 
Seattle,  WA  98124 

Attn;  Library 

A o re r. i  >a .;e  Cftrpo t  \ on 
(  .  LI  Segurdo  P 1  vd 
Los  Angeles ,  CA  90045 

At  tn ■  l  ibrary 

lot  M.<  1  Missil..  and  V, 

1111  I  or  k  hoc  1  Vs,iy 
Sunnyvale,  CA  94088 

At  tn .  L ibrary 

KiirmmrGt  Corporation 
1661*6  Saticey  Strict 
Be*  x  2013  South  Anne  a 
Van  Nuys ,  CA  91409 

Attn;  Library 

Martin  Marietta  Aerospa 
P.0.  Pox  179 
Denver ,  CO  80201 

Attn:  Library 

Me  bonne  1 1  Doug  1  ■» «,  As  f  *-. . 

5301  Pol  sa  A , r 

Hunt  irigi <m  [•••K.n,  LA  n 

Attn:  Library 

Sp  it  e  Di  v i  s  ior 
km  Uell  If*  .  n-t  n»*  I 
12, '14  I  .ilrw.  ■>.!  P.ivd. 
y  ,  CA  r--.. 

|  i Sr  1 1 

til  r-r  ■  f  | :  fr,.  ...  , 

Beth,  Vy.  ’Pi  1  1  7 1 4 
Altr.  I 

Or  n«  ■  * !  ['  yr  .,  ■  r  ■  f  <  , ,  ,  i 

I.O.  P.,,.  H  • 

!  i‘  1--.  FA 

At  * »» .  lit*-  i  r  y 

f.  ■  I  .. . 

'  (  .  r  -  - 

C  •  n  J  .  ■  ,  i  *  •  .  -  *  -  ■  • 


I  it  • 


nt  *.  rn  1 1  ionj  1 


•if*  Group 
vision 

■  .'  -r.it  ion 


hi;  f oinu.il>} 


ce 


nautics 
.  04  7 

(  urn., rat  i  on 


